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ABSTRACT  (Continued) 


authority  digital  electronic  controller  must  utilize  demonstrated  multi- 
variable  design  techniques  to  integrate  adequately  these  complex  system 
functions. 

A  preliminary  design  of  a  controller  for  this  variable  cycle  engine 
is  described.  It  is  implementable  on  a  small  digital  computer  (less  than 
16K  words  of  storage),  and  is  modular  in  design  (subroutine  format). 
Specific  controller  functions  of  transient  regulation,  steady  state 
regulation,  trajectory  generation,  signal  processing,  and  fault  detection 
and  accommodation  are  incorporated  in  a  way  which  allows  experimentation 
with  different  techniques  for  each  function  without  affecting  the  overall 
structure.  Promising  techniques  for  implementing  each  function  are  dis¬ 
cussed. 


FOREWORD 


The  work  described  in  this  report  was  performed  under 
Contract  No.  F33615- 77-C-2096  for  the  Air  Force  Aero  Pro¬ 
pulsion  Laboratory  (AFAPL),  Wright-Patterson  Air  Force  Base. 
The  technical  monitor  for  AFAPL  is  Mr.  Charles  A.  Skira. 
Principal  authors  from  Systems  Control,  Inc.  (Vt)  are  Dr. 
Stephen  M.  Rock  and  Dr.  Ronald  L.  De  Hoff.  Program  manager 
for  Variable  Cycle  Engine  Multivariable  Control  Synthesis  is 
Dr.  W.  Earl  Hall,  Jr.  Computer  programming  support  was 
provided  by  Ms.  Connie  G.  Hayes.  Report  preparation  efforts 
were  directed  by  Ms.  Clare  Walker. 


iii 


TABLE  OF  CONTENTS 


Page 


I.  EXECUTIVE  SUMMARY  . .  1 

1.1  Introduction  .  1 

1.2  Multivariable  Control  Design  Program  for  the 

GE-23  Variable  Cycle  Engine  .  3 

1 . 3  Summary .  5 

II.  CONTROL  REQUIREMENTS  FOR  A  VARIABLE  CYCLE  ENGINE  .  7 

2.1  Introduction  .  7 

2.2  The  UE23-JTDE  Variable  Cycle  Engine  .  7 

2.3  Control  Hardware  Characteristics  .  10 

2.4  Mission  Definition  .  16 

2.5  Control  Requirements  .  19 

2.6  Summary . . .  22 

III.  MULTIVARIABLE  CONTROL  DESIGN  METHODOLOGY .  24 

3.1  Introduction .  2** 

3.2  The  Control  Design  Problem .  24 

3.3  Feedforward  and  Transition  Control  .  .  25 

3.4  Classical  Control  Synthesis  ....  .  29 

3.5  Modern  Control  Synthesis  Methods  .  .  .  29 

3.6  Model  Generation . 39 

3.7  Nonlinear  Model  Development  .  52 

3.8  Summary . 57 

IV.  CONTROL  STRUCTURE  DEFINITION  ....  .  58 


4.1  Introduction . . .  . 

4.2  Overview  of  System . . 

4.3  Nonlinear  Model . .  . . 

4.4  Reference  Point  Schedules  .......... 

4.5  Trajectory  Generator  .....  . 

4.6  Multivariable  Regulator  Logic  ........ 

4.7  Fault-Tolerant  Filter . .  .  .  . 

4.8  Engine  Limit  Logic  . . .  .  . 

4.9  Failure  Accommodation  Logic  .  . 

4. 10  Summary  ....  . 

V.  PRELIMINARY  DESIGN  RESULTS  ............ 


$8 

63 

63 

69 

11 

81 

84 

85 


a*. 

a*  » 

94 


5.1  Linear  Model  Analysis  . .  94 

5.2  Nonlinear  Model  Generation  . . .  118 

5.3  Actuator  Servo  Valve  Compensation  .  122 

5. 4  Transition  Generator  Demonstration  .  134 


v 


TABLE  OF  CONTENTS  (Continued) 


Page 


VI.  SUMMARY  AND  PROGRAM  STATUS  .  138 

6.1  Summary . 138 

6.2  Program  Status . 139 

APPENDIX  A:  DERIVATION  OF  141 

APPENDIX  B:  SECOND  STAGE  MODEL  REDUCTION  .  147 

APPENDIX  C:  MODELING  PROCEDURES  FOR  THE  VARIABLE  CYCLE 

ENGINE . 151 

REFERENCES . 160 


LIST  OF  FIGURES 

Page 

1.1  Variable  Cycle  Engine  .  3 

1.2  Multivariable  Control  Design  Program  for  the  GE-23 

Variable  Cycle  Engine  .  4 

2.1  Variable  Cycle  Engine  -  Controls  and  Outputs  .  .  8 

2.2  Grumman  Design  623  VSTOL-B  Fighter  .  10 

2.3  Functional  Hardware  Description  .  10 

2.4  Control  System  Diagram  .  12 

2.5  Fail-Fixed  Servo  Valve  .  15 

2.6  Grumman  623  Envelope . 16 

2.7  Three  Typical  Mission  Profiles  for  the  GE23/623 

Involving  Air  Combat,  VTOl  Intercept  and  Subsonic 
Surveillance  .  17 

2.8  Partition  of  Operating  Envelope  in  Mission  ...  18 

3.1  Block  Diagram  of  Output  Regulator  Design  ....  kl0 

3.2  Generation  of  Design  Models  . 

3.3  Comparison  of  Time  Histories  for  Second  Stage 

Reduction  Example  .....  . 51 

4.1  Controller  Structure  .  59 

4.2  FADEC  Hardware  .  . . 51 

4.3  Control  Signal  Flow  Diagram  .  65 

4.4  Reference  Point  Schedule  Generator  .......  70 

4.5  Control  Criteria  within  Flight  Envelope  .....  72 

4.6  Control  Saturation  Effects  on  Nj?  Schedule  ....  75 

4.7  Point  Calculation  Flowchart  .....  75 

4.8a  trajectory  Generator  Logic  -  Large  Transient 

Control  Loop . 79 

4,8b  Trajectory  Generator  Logic  -  Small  Transient  Loop  79 

4.9a  T.  Limit  Schedule  .  . . 82 

*4 

4.9b  TBj  Limit  Schedule  for  Increased  Engine  Life  .  ,  82 

4.10  Engine  Limit  Logic . 87 

4.11  Overview  of  Engine  Limit  Logic  . . 88 

4.12a  Normal  Operation  of  Actuator  Servo  Loop  .  91 

4.12b  Operation  of  Actuator  with  Failed  LVPT . 91 

S.la  Definition  of  Flight  Points  vs.  Altitude  and 

Mach  Number . 96 

S.lb  Definition  of  Flight  Points  vs,  PT2  and  TT2  .  .  ,  96 

vii 


LIST  OF  FIGURES  (Continued) 


Page 


5.2  Definition  of  Power  Points  1  through  5  and  A 

through  D  . 

5.3  Variable  Cycle  Engine  -  Controls  and  Outputs  .  . 

5.4  Differential  and  Collective  Spool  Response  Time 

Constants  at  Flight  Points  1  through  6  . 

5.5  Temperature  Response  Time  Constants  at  Flight 

Points  1  through  6  . 

5.6  Eigenvalues  Corresponding  to  Engine  Gas  Dynamics 

at  Each  of  Six  Flight  Conditions  . 

5.7  Eigenvalues  (Gas  Dynamics)  vs.  Flight  Condition 

5.8  Eigenvector  Composition  of  the  Five  Princpal 

Modes  at  SLS/Maximum  Power  (Table  5.7)  . 

5.9  Residue  Spectra  of  Selected  Transfer  Functions.  . 

5.10  Sensitivity  of  Mean  Square  Response  to  Elements 

in  F-Matrix  .....  . 

5.11  Nonlinear  Model  Generation  Procedure  ...... 

5.12  Comparison  of  Linear  Eigenvalue  Locations  and  Root 
Locations  Found  by  Linearising  the  Nonlinear  Model 
at  28  Flight  Points  Spanning  the  Envelope  .... 

5.13  Comparison  of  Idle  to  Military  Power  Step 
Simulated  using  Linear  Model  Derived  at  Inter¬ 
mediate  Power,  the  Nonlinear  Model  Derived  from 
Linear  Matrices  and  a  Model  using  Linearisations 
at  Eaeh  Point  along  the  Trajectory  .  .  .  .  .  .  . 

5.14  Failed- Fixed  Servo  Valve  Characteristic  . 

5.15  Revised  Failed-Fixed  Servo  Valve  Characteristic  . 

5.16  Block  Diagram  of  Actuator  Valve  Servoloop  .  .  .  . 

5.17  Modulation  Interval,  Sample  Inteval  and  Computa¬ 
tion  Interval  . . . . 

5.18  Flowchart  of  Actuator  Valve  Servo  Loop  Simulation 

5.19  Nonlinear  Compensation  for  Actuator  Servovalve 

(FADEC  I) . . 

5.20  Step  Response  of  F'ADEC  Valve  #1  to  a  101  of  Full 
Scale  Step  in  Commanded  Stroke  ......... 

5.21  Step  Response  of  FADEC  Valve  *2  to  a  10$  of  Full 

Scale  Step  in  Commanded  Stroke  .  . . .  . 


97 

103 

107 

107 

108 
108 

109 

111 

115 

118 

120 


121 

123 

123 

12S 

12? 

128 


130 


131 

133 


viii 


LIST  OF  FIGURES  (Concluded) 

Page 


5.22  Comparison  of  Step  Response  and  Linear  Servomech¬ 
anism  Response  for  a  101  Acceleration  to  SLS  .  .  135 

5.23  Demonstration  Trajectory  Generator  for  VCE  Linear 

Simulation  Using  a  Fuel  Flow  Actuator  and  Two 
Nonlinear,  Variable  Rate  Feedbacks  .  136 

5.24  Demonstration  of  Temperature  and  Surge  Limiting 
Acceleration  Trajectory  Using  Nonlinear,  Variable 

Rate  Feedback .  137 


ix 


LIST  OF  TABLES 


Page 


2.1  Sensor  Hardware .  13 

2.2  Actuator  Characteristics  .  14 

2.3  Flight  Domain  Related  Control  Criteria  .  19 

2.4  Engine  Limits .  20 

2.5  Error  Tolerances  on  Set  Points .  21 

2.6  Failure  Tolerance  Specification  .  22 

3.1  Classical  Control:  Advantages  and  Disadvantages  ...  30 

3.2  Modern  Control:  Advantages  and  Disadvantages  .....  32 

4.1  Input/Output  Signal  Categories  .  62 

4.2  Preliminary  Functional  Description  of  Multivariable 

Control  . .  69 

4.3  Mission/Aircraft  Definition  (Control  Performance 

Criteria)  .  . . 71 

4.4  Operating  Constraints  .  73 

4.5  Modes  of  Operation  .  > .  71* 

5.1  Definition  of  Power  Points  . . 93 

5.2  Summary  of  Sixty  Linearisation  Points  . 19 

5.3  Engine  State  Variables  .  102 

5.4  Engine  Inputs  . .  102 

5.5  Engine  Output  Variables  .  . . . . 102 

5*6  Eigenvalues  of  VCE  at  SLS  -  Maximum  Power . 106 

5.7  Principal  Dynamic  Modes  of  VCE  at  SLS  ~  Maximum  Power  .  106 

5.8  State  Variables  for  Reduced  Order  System  . .  110 

$.9  Comparison  of  Eigenvalues*  between  F  and  F  . . 117 

$.10  Preliminary  Nonlinear  Model  Structure  .........  US 


x 


SECTION  I 
EXECUTIVE  SUMMARY 


1.1  INTRODUCTION 

Aircraft  turbine  engine  propulsion  control  systems  have  been 
the  focus  of  extensive  development  in  recent  years.  Improvements 
in  all  types  of  engine  components  have  been  realized  by  optimized 
materials,  and  configuration,  and  the  addition  of  variable  geometry. 
Engine  control  hardware  technology  has  accordingly  progressed  to 
utilize  effectively  the  nearly  unlimited  capabilities  of  digital 
processing  machines.  This  progress  has  occurred  so  rapidly  that 
demonstrated  theoretical  synthesis  techniques  using  the  vastly 
increased  performance  potential  have  been  few.  Recently,  programs 
sponsored  by  the  Air  Force  Aero  Propulsion  Laboratory  and  NASA 
lewis  Research  Center  have  begun  to  exploit  this  area  [1,2]. 

The  AFAPL/NASA  LeRC- sponsored  FIDO  multivariable  controls  program 
is  an  example  of  the  successful  demonstration  of  this  emerging 
technology  on  a  current  production  aircraft  turbine  [2,3]. 

Future  engine  cycles  have  been  defined  and  include  a  signifi¬ 
cant  variable  geometry  capability  to  the  extent  that  the  thermo¬ 
dynamic  engine  cycle  varies  over  wide  limits.  This  new  variable 
cycle  engine  (VCE)  definition,  for  example,  includes  the  latest 
engine  technology.  The  number  of  actuated  variables  have  been 
similarly  expanded  and  resulting  control  strategies  are  sufficiently 
complex  to  require  extension  of  the  current  advanced  multivariable 
control  methodologies. 

Military  aircraft  missions  of  the  1980‘s  and  beyond  are  dic¬ 
tating  the  requirement  for  multimode  integration  of  airframe  and 
propulsion  systems  over  a  wide  flight  envelope.  In  most  present 
designs,  the  multimode  requirement  is  being  met  by  varying  aircraft 
configurations  or  relying  on  active  controls  technology  to  accommo¬ 
date  response  changes  to  varying  flight  conditions.  It 
.  is  becoming  evident  that  much  of  the  multimode  capability  must 
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be  met  by  the  propulsion  system  itself.  The  apparent  tactical 
effectiveness  of  vectored  thrust  further  highlights  the  need  for 
multimode  performance  capability  to  be  at  least  partially  satis¬ 
fied  by  the  propulsion  system. 

The  provision  of  multimode  propulsion  response  capability, 
without  sacrificing  efficiency  or  performance,  is  achievable  by 
engines  with  variable  cycles  which  are  controlled  by  commanded  in¬ 
ternal  geometrical  changes  in  the  engine  itself.  Such  control 
capability  is  achieved  at  the  expense  of  a  significant  increase  in 
engine  complexity,  addition  of  actuators,  and  addition  of  sensors. 
The  subsequent  control  system  complexity  to  maintain  strict  tran¬ 
sient  and  steady  state  performance  specifications  forces  attention 
to  more  accurate  and  reliable  controller  implementations.  The 
evaluation  of  engine  control  technology  is  showing  that  such  com¬ 
plex  engine  control  functions  could  be  implemented  with  digital 
electronic  techniques,  relying  on  hydromechanical  control  hardware 
for  backup,  fail/operate  functions  [4].  Control  synthesis 
.  techniques  for  such  digital  systems  are  now  being  developed  and 
^demonstrate  multivariable  design  techniques  can  integrate  ade¬ 
quately  the  complex  coupling  of  the  system,  its  functions,  modes, 
and  performance  requirements. 

The  variable  cycle  turbine  engine,  GE23-JTDE,  represents  a 
prototype  of  future  multimode  propulsion  plants.  It  is  a  sophis¬ 
ticated  design  of  highly  variable  geometry  and  multiple  control 
inputs.  To  control  such  an  engine,  a  large  number  of  engine  vari¬ 
ables  must  be  sensed  including  engine  pressures,  temperatures, 
rotor  -speeds',  and  airframe  and  inlet  commands.  A  controller  for 
this-  engine  must  therefore  be  delineated  as  multivariable  (i.e., 
manipulating  large  numbers  of  input  and  output  variables)  and 
multifunctional  (i.e*.  perform,  in  addition  to  control  logic,  also 
data : conditioning  and  fault  diagnosis).  The  development  of  such  a 
full  authority  digital  electronic  controller  must  therefore  utilise 
demonstrated  multivariable  design  techniques  to  integrate  ade¬ 
quately  these  complex  system  functions. 


1.2  MULTIVARIABLE  CONTROL  DESIGN  PROGRAM  FOR  THE  GE-23  MRTABLE 

CYCLE  ENGINE 

Systems  Control,  Inc.  (Vt)  is  conducting  an  advanced  develop¬ 
ment  program  for  the  Air  Force  Aero  Propulsion  Laboratory  to  de¬ 
sign  a  multivariable  controller  for  an  advanced  variable  cycle 
engine.  An  associated  contractor,  the  General  Electric  Co.,  Air¬ 
craft  Engine  Group,  is  providing  modeling  and  simulation  data  des¬ 
cribing  the  Joint  Technology  Demonstrator  Engine,  the  GE-23. 

The  GE-23  was  selected  as  the  test  bed  for  advanced  VCE  multi- 
variable  control  design  because  cf  the  availability  of  digital 
hardware  capability  for  direct  application  of  digital  control  logic 
An  engine  cross  section  is  shown  in  Figure  1.1  [5].  The 
engine  actuated  variables  include  two  variable  fuel  flows  and  eight 
variable  geometries.  It  is  a  double-bypass,  single  no2zle,  mixed- 
flow,  augmented  variable-cycle  turbofan  with  forward  and  aft 
variable-bypass  injectors.  Both  fan  blocks  and  turbines  are 
single  stage.  There  are  va  C'e  inlet  guide  vanes  and  variable- 
area  turbine  nozzles. 
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The  development  program  (Figure  1.2)  consists  of  a  Phase  I  con¬ 
trol  definition  activity  and  a  Phase  II  design  demonstration.  This 
report  details  the  results  of  the  Phase  I  control  structure  defini¬ 
tion  and  presents  selected  examples  of  the  application  of  the 
multivariable  control  design  methodology.  The  development  of  a 
demonstrator  MVC  for  the  GE-23  engine  will  be  accomplished  in 
Phase  II.  Validation  of  the  logic  will  be  undertaken  on  a  detailed 
nonlinear  digital  simulation  of  the  engine.  The  controller  will 
then  be  implemented  on  a  hybrid  simulation  facility  at  General 
Electric,  Evandale,  Ohio,  for  detailed  evaluation  and  validation 
of  the  performance  throughout  the  flight  envelope. 


The  successful  completion  of  these  activities  will  represent 
a  major  accomplishment  in  the  application  of  modern  control  tech¬ 
nology  to  advanced  cycle  turbine  propulsion  systems.  The  realiza¬ 
tion  of  improved  engine  performance  and  reliability  relies  heavily 


on  the  control  algorithm  which  effectively  operates  the  engine  to 
satisfy  the  stringent  limiting  criteria  and  performance  require¬ 
ments  in  the  presence  of  failures  and  degradation.  It  is  the  goal 
of  this  program  to  demonstrate  the  flexibility  and  effectiveness 
of  modern,  digital  control  design  methodology  in  solving  the  dif¬ 
ficult  problems  associated  with  the  turbine  engine  application  and 
prototype  a  system  for  installation  in  the  future  engine  implemen¬ 
tation. 


1 . 3  SUMMARY 


This  report  details  the  specification  of  a  multivariable  con¬ 
trol  structure  for  the  GE-23  variable  cycle  engine.  This  report 
represents  a  self-contained  document  which  describes  the  design 
methodology  which  has  been  developed  to  address  the  controls  cri¬ 
teria,  mission  specification  and  hardware  definition  for  operation 
of  the  engine  within  a  typical  installed  envelope.  The  report  is 
meant  to  provide  a  reference  to  the  design  process.  The  final 
report  will  describe  the  synthesis  an*  evaluation  of  a  control 
using  the  procedures  presented  below,  for  the  GE-23.  The  anticipated 
computer  storage  requirement  for  this  control  is  between  8K  and 
12K  words  of  memory. 

The  report  is  organized  as  follows: 

•  Section  II  Control  Requirements  for  th 1  Variable 

Cycle  Engine  ~ 

This  chapter  presents  the  definition  of  the  variable  cycle 
engine  (VCB)  control  problem  for  the  GE-25.  Engine  hard¬ 
ware  and  interfaces  are  described.  The  installed  mission 
definition  is  presented  and  the  contrr?  design  criteria 
are  established.  Each  of  these  elements  is  considered 
during  the  formulation  of  the  control  *t,ucture  and  the 
underlying  synthesis  methodology. 

•  Section  III  Multivariable  Control  Design  Methodology 

This  chapter  details  the  synthesis  procedures  which  are 
required  for  development  of  the  multimode  controller.  The 
calculation  of  the  control  law,  generation  of  trajectories 
and  provision  for  failure  tolerance  and  accommodation  are 
treated  from  a  theoretical  viewpoint. 
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•  Section  IV  Control  Structure  Definition 


The  functional  definition  of  the  control  logic  is  presented 
in  Chapter  IV.  The  reference  schedules,  trajectory  gener¬ 
ator,  multivariable  control  law,  engine  protection  and 
fault  tolerant  filter  modules  are  described.  The  synthesis 
of  these  blocks  is  related  to  the  theoretical  procedures 
presented  in  Chapter  III  and  a  detailed  specification  of 
the  prototype  system  and  alternatives  is  delineated. 


Section  V 


Preliminary  Design  Results 


Several  procedures  which  will  be  utilized  in  the  final  con¬ 
trol  design  have  been  exercised  on  preliminary  data  sup¬ 
plied  by  the  engine  manufacturer.  Results  concerning 
linear  model  analysis,  nonlinear  model  development,  control 
design  methods  and  actuator  compensation  are  presented  and 
evaluated  relative  to  the  theoretical  goals  and  ultimate  design 
applications. 


Section  VI  Summary  and  Conclusions 


Chapter  VI  summarizes  the  control  structure  and  describes 
the  design  and  synthesis  activities  being  undertaken  in 
Phase  II. 


This  report  represents  interim  progress  in  a  major  design 
effort.  At  this  point  in  time,  it  appears  that  the  control  design 
procoss  which  will  bo  a  product  of  this  study  can  have  major  and 
lasting  impact  on  the  production  of  the  next  generation  high  per¬ 
formance  aircraft  propulsion  plants. 
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SECTION  II 

CONTROL  REQUIREMENTS  FOR  A  VARIABLE  CYCLE  ENGINE 

2.1  INTRODUCTION 

The  next  generation  aircraft  weapon  system  will  achieve  higher 
levels  of  performance,  versatility  and  effectiveness  than  its 
predecessors  at  the  cost  of  increased  complexity.  The  variable 
cycle  engine  incorporates  larger  numbers  of  actuated  and  sensed 
variables  to  meet  these  advanced  mission  requirements.  This  chapter 
describes  a  prototype  variable  cycle  engine  which  will  constitute 
the  test  bed  for  demonstration  of  the  multivariable  control  design 
methodology. 

The  General  Electric  GE23-JTDE  engine  is  described.  This 
engine  will  be  controlled  by  an  engine-mounted  microelectronic 
digital  computer  being  developed  under  the  Full  Authority 
Digital  Electronic  Control  (FADEC)  program  by  the  General 
Electric  Company. 

Mission  definition  and  installation  requirements  are  speci¬ 
fied  f*r  an  advanced,  VSTOL-B  military  fighter  concept,  the 
Grumman  623.  Control  requirements  are  described  for  performance, 
response,  limit  protection  and  failure  accommodation  within  the 
framework  of  the  above  control  hardware  and  installed  environment. 
This  chapter  describes  the  control  criteria  being  addressed  by 
the  multivariable  control  design.  Further  information  concerning 
the  design  data  can  be  found  in  Ref.  6. 

2.2  THE  GE23-JTDE  VARIABLE  CYCLE  ENGINE 

The  GE23  is  a  twin-spool  turbofan  engine  (see  Figure  2.1). 

The  compressor  and  high  pressure  turbine  form  the  core  assembly. 
Variable  compressor  stators  (STP25)  and  turbine  area  (STP49)  can 
be  used  to  modulate  the  core  operating  characteristics.  The  low 
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Figure  2.1  Variable  Cycle  Engine 
Controls  and  Outputs 


pressure  turbine  drives  a  two- stage  fan.  A  forward  bypass  valve 
(AE94)  provides  the  capability  to  bypass  airflow  around  the 
second  fan  block.  Variable  guide  vanes  (STP22)  are  provided 
for  the  second  fan  block.  The  low  pressure  turbine  area  (STP49) 
can  be  modulated  to  adjust  turbine  work  output.  Variable  area 
elements  are  included  at  the  augmentor  entrance  (AE16)  and 
exhaust  nozzle  exit  (A8)  to  provide  further  operational 
flexibility. 


The  JTDE's  variable  geometry  gives  it  the  ability  to  oper¬ 
ate  over  a  significant  power  range  at  constant  airflow.  Cycle 
pressure  and  temperature  ratio  is  varied  to  produce  this  oper¬ 
ation.  Conventional  turbojets  or  turbofans  modulate  thrust  by 
reducing  airflow  along  with  pressure  and  temperature  ratio. 

At  intermediate  power  setting  the  JTDE  operates  as  a  low 
bypass  ratio  mixed  flow  augmented  turbofan.  The  forward  bypass 
transfer  valve  is  closed,  the  forward  bypass  injector  is  open  and 
the  engine  is  operating  at  its  maximum  pressure  ratio.  All  air 
flow  is  compressed  by  both  blocks  of  the  fan.  As  power  is  reduced, 
pressure  ratio  can  bo  lowered  by  opening  the  nozzle  or  by  bypass¬ 
ing  air  around  the  second  block  of  the  fan.  Using  control  of  the 
forward  injector  and  the  fan  variable  geometry,  the  second  block 
can  be  unloaded  without  changing  airflow  or  the  operating  condi¬ 
tions  of  the  front  block  of  the  fan.  Fan  turbine  nozzle  area  can 
be  increased  to  reduce  its  power  to  match  the  reduced  power  required 
to  drive  the  fan  in  the  double  bypass  mode.  The  engine  now  oper¬ 
ates  like  a  higher  bypass  ratio  turbofan  for  subsonic  cruise  power 
conditions.  The  aft  bypass  injector  adjusts  to  accommodate  the 
increased  ratio  of  bypass  air  to  core  air.  The  engine  also  retains 
the  ability  to  operate  at  high  specific  thrust  where  high  power 
operation  is  necessary,  i.e.,  takeoff  and  high  Mach  number  opera¬ 
tion. 

The  installation  definition  chosen  for  this  study  is  the 
Grumman  design  623  advanced  VSTOL-B  fighter.  The  aircraft 
(Figure  2.2)  uses  two  GE23  engines  to  operate  from  vertical 
take-off  to  altitudes  in  excess  of  SO, 000  feet  and  speeds  in  ex¬ 
cess  of  Mach  2.  The  623  is  an  advanced  multimission  fighter 
aircraft  having  both  air-to-air  and  air-to-ground  missions  requir¬ 
ing  the  multimode  capability  of  the  variable  cycle  engine. 
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Figure  2.2  Grumman  Design  623  VSTOL-B  Fighter 


2.3  CONTROL  HARDWARE  CHARACTERISTICS 


The  on-engine  control  hardware  consists  of  engine  actuators, 
sensors,  a  microcomputer  and  power  supply,  hardware  failure 
overrides  and  a  hydromechanical  backup  control.  The  FADEC 
computer  is  a  fuel-cooled  LSI  processor  with  specially  designed 
interface  hardware  for  the  actuators  and  sensors.  The  functional 
interconnection  of  these  hardware  elements  is  shown  in  Figure  2.3. 


Figure  2.3  Functional  Hardware  Description 


Inputs  to  the  control  computer  are  required  from  engine 
sensors,  actuator  position  feedbacks  and  the  airframe  data  bus 
(see  Figure  2.4).  Data  from  the  airframe  includes  power  level 
command,  PLA,  and  aircraft  sensor  data-including  Mach  number,  ambient 
temperature  and  pressure,  thrust  mode,  wheels  down  and  squat 
switch  logic.  Data  can  be  received  from  the  inlet  control  sys¬ 
tem  on  this  communication  link;  however,  the  current  Grumman 
design  does  not  require  this  feature. 

Actuator  position  feedback  signals  are  provided  by  linear, 
variable  phase  transducers  (LVPT)  which  measure  stroke  on  the 
servo  pistons.  These  instruments  measure  phase  shift  between 
excitation  and  return  signals  passing  through  a  variable  induc¬ 
tance  device.  A  digital  counter/threshold  detector  is  used  to 
time  the  phase  difference  and  calculate  the  measured  stroke.  The 
phase  converters  are  free  running,  and  measured  values  are 
asynchronously  loaded  into  the  control  processor.  Error  sources 
on  these  signals  are  due  primarily  to  load  defloction  on  the 
linkages  and  effective  area  changes  in  the  variable  geometry 
hardware.  Ten  actuator  stroke  positions  are  used. 

The  sensor  inputs  used  by  the  control  are  listed  in  Table  2.1. 
Pressure  measurements  use  vibrating  crystal  transducers.  Frequency 
conversion  of  these  signals  provides  asynchronous  input  to  the 
computer  memory.  Static  probes  are  located  at  the  discharge  sta¬ 
tions  of  the  two  fan  blocks  and  the  high  pressure  compressor.  Two 
differential  probes  are  located  behind  each  fan  block.  These 
probes  produce  high  response  signals.  Error  sources  arise  from 
rotating  blade  passage,  stall  pressure  variations,  hysteresis, 
temperature  effects  and  distortion.  Accuracy  levels  for  those 
devices  are  in  the  0.54  range  (except  near  stall). 

Accurate  rotor  speed  measurements  (<0.14)  are  obtained  from 
frequency  conversion  of  an  alternator  output  (high  spool)  and 
magnetic  pickup  (low  spool).  Converted  measurements  are  asynchro¬ 
nously  input.  Time  delays  on  these  signals  may  exceed  25  ms 
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Table  2.1 
SENSOR  HARDWARE 


PARAMETER 

SENSOR 

TYPE 

PS21  Absolute 

APjj  Differential 
Pj23  Absolute 
^23  Differential 

PS3  Absolute 

Qua 

'U 

T12 

T2.5 

RTD 

{ 

TU 

J-R 

llght-Off  Detector 

IM 

Alternator 

•V 

Magnetic 

Pickup 

at  low  power.  This  effect  is  caused  by  counter  delays  at  the 
low  shaft  frequencies. 

Two  analog  signals  measuring  fan  block  discharge  temperature 
are  sensed  through  A-D  converters.  Platinum  resistance  thermo- 
meters  (RTD)  provide  an  accurate  linear  temperature/resistance 
relationship  which  is  converted  to  an  analog  voltage  and  sampled 
by  the  computer.  These  devices  have  an  airflow  dependent  time 
constant  which  can  exceed  eight  (8)  seconds  at  high  altitude  * 
low  power  conditions.  An  optical  pyrometer  measures  high  pres¬ 
sure  turbine  blade  temperature.  The  device  senses  infrared  emis¬ 
sions  from  the  rotating  elements.  The  response  time  of  this  probe 
is  extremely  fast;  however*  the  blade  temperature  generally  lags 
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the  gas  temperature  at  this  station  by  several  seconds  due  to 
mass  effects  in  the  rotor. 

An  augmentor  light-off  detector  senses  ultraviolet  emis¬ 
sions  from  burning  fuel  in  the  tailpipe.  This  device  produces 
a  nearly  instantaneous  discrete  indication  of  fuel  ignition 
during  augmented  operation. 

Ten  actuators  are  used  to  modulate  areas  and  flows  in  the 
GE23.  These  actuators  are  described  in  more  detail  in  Chapter  V 
and  Ref.  6.  A  spool  valve  mechanism  is  used  to  modulate  a  fuel 
or  hydraulically  driven  servo  piston.  The  actuated  elements 
are  listed  in  Table  2.2.  The  spool  valve  control  input  is  a 
pulse  width  modulated  current  from  a  FADEC  output  amplifier. 

The  spool  valve  is  designed  to  provide  a  fail/fixed  operating 
characteristic  as  shown  in  Figure  2.5. 

Table  2.2 

Actuator  Characteristics 


VMHAftlE 
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Figure  2.5  Fail -Fixed  Servo  Valve 

If  a  modulator  failure  or  internal  mechanical  failure 
drives  the  spool  to  one  extreme  or  the  other,  the  actuator  is 
fixed  in  its  last  position.  Leakage  around  the  spool  lands 
will  cause  the  actuator  to  drift  to  a  predetermined  setting 
after  a  failure  has  occurred.  This  control  actuator  inter¬ 
face  device  provides  an  extremely  useful  tool  for  incorporat¬ 
ing  control  logic  to  accommodate  actuator  failures  since  the 
necessity  of  compensating  for  fast,  hardover  response  has  been 
eliminated. 

The  hardware  configuration  of  the  Oh 2 5  incorporating  control 
electronics,  sensor,  hydromechanics!  actuators  and  backup 
control  provide  the  control  designer  with  a  flexible  and  power¬ 
ful  capability  to  control  the  engine  reliably  and  efficiently 
throughout  its  envelope.  This  system  represents  a  thoroughly 
integrated  concept  for  failure  compensation  and  fault  tolerance 
which  can,  with  proper  control  software,  produce  an  overall 
design  with  hydromechanical  reliability  and  digital  electronic 
capability. 


IS 


2.4  MISSION  DEFINITION 


The  aircraft  mission  requirements  are  an  important  factor 
in  determining  the  functional  structure  of  t>.-c  controller  and 
the  complexity  of  the  implementation.  The  Grumman  623  is  an 
advanced  VSTOL-B  fighter.  The  aircraft  flight  envelope  is 
shown  in  Figure  2.6.  The  aircraft  has  both  air-to-air  and  air- 
to-ground  missions,  three  of  which  have  been  chosen  for  considera¬ 
tion  during  the  study.  The  primary  mission  is  fighter  escort-air 
combat.  The  flight  profile  is  shown  in  Figure  2.7.  In  addition, 
a  surface  surveillance  and  VTOL  intercept  mission  were  considered. 
These  profiles  are  also  shown  in  Figure  2.7. 


Figure  2.6  Grumman  623  Envelope 

The  mission  definitions  determine  the  type  of  engine  per¬ 
formance  which  is  important  and  the  region  of  the  power  and 
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flight  domain  in  which  this  performance  is  required.  The  flight 
envelope  can  bo  partitioned  into  regions  corresponding  to  each 
of  the  mission  modes.  These  regions  are  shown  in  Figure  2.8. 

By  considering  the  important  operational  characteristics  of 
each  domain,  a  definition  of  control  design  criteria  can  be 
developed.  Table  2.3  lists  the  control  criteria  which  are  associ¬ 
ated  with  each  flight  domain.  The  translation  of  these  per¬ 
formance  criteria  into  operating  point  specifications,  reference 
point  values  and  mode  definitions  is  presented  in  Chapter  IV. 


mach  mm* 

Figure  2.8  Partition  of  Operating  Envelope  into  Mission  Domains 
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Table  2.3 

Flight  Domain  Related  Control  Criteria 


DOMAIN 

CONTROL  CRITERIA 

Maximum  Thrust 

Maximize  augmented  and  dry 
thrust  at  design  points 
with  minimum  SFC  between 
points 

Cruise/Loiter 

Minimum  SFC  at  design 
thrust 

Maximum  Stability 

Maximum  augmented  and  dry 
thrust  at  design  points 
with  Increased  component 
stability  margins 

Maximum  Response 

Maximum  thrust  response 
at  part  power 

Maximum  Life 

Rated  thrust  with  minimum 

temperature/speed 

variation 

2.5  CONTROL  REQUIREMENTS 

In  addition  to  the  performance  definition  discussed  in 
Section  2.4,  the  controller  must  protect  the  engine  from 
limit  exceedances,  accurately  specify  the  set  points  and 
accommodate  failures.  These  design  specifications  are  presented 
below. 

2.S.1  Engine  Limit  Protection 

Engine  limits  arise  in  components  operating  in  a  high  pres¬ 
sure,  high  temperature  environment  (see  Table  2.4).  Structural 
considerations  of  fatigue,  creep  and  strength  dictate  temperature 
and  rotational  speed  limits  for  the  high  pressure  turbine  entrance 
and  the  compressor  discharge  stations.  In  addition,  life  limits 
are  imposed  on  turbine  blade  temperature  rates  during  power  changes. 
Structural  limits  are  approached  for  combustor  pressures  also. 
Combustion  stability  must  be  maintained  by  limiting  fuel-to-air 
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Table  2.4 
Engine  Limits 


•  FANS  AND  COMPRESSOR  LIMITS 

Maximum  Discharge  Pressure 
Maximum  RPM 
Maximum  Corrected  RPM 
Minimum  Surge  Margin 

Maximum  Closure  of  Stators. from  Basic  Schedule 
Maximum  Opening  of  Stators  front  Basic  Schedule 

•  COMBUSTOR  AND  AFTERBURNER  LIMITS 

Rich  Burner  Blowout  Limit 
Lean  Burner  Blowout  Limit 
Minimum  Pressure 
Llghtoff  Limits 

e  TURBINE  LIMITS 
Maximum  RPM 

Maximum  Inlet  Qas  Temperature 

Maximum  Turbine  Blade  Temperature 

Maxlmua  Rate  of  Change  of  Turbine  Blade  Temperature 

e  BYPASS  DUCT 

Maximum  Airflow  or  Velocity 
Minim*  Airflow  or  Velocity 


ratios  between  lean  and  rich  limits  and  maintaining  main  and  aug- 
mentor  combustion  pressures  above  a  minimum  value.  Plow  stability 
in  the  fan  blocks  and  compressor  is  maintained  by  assuring  surge 
margin  limits  [6]  are  maintained.  Duct  airflow  must  be  main¬ 
tained  between  Mach  number  limits  to  ensure  stable,  unchoked  flow. 

Each  of  the  engine  limits  described  above  is  approached 
during  transient  or  steady  state  operation.  For  example,  the 
engine  will  most  efficiently  operate  at  many  intermediate 
power  conditions  at  the  maximum  turbine  temperature  limits.  The 
limiting  values  of  the  variables  represent  desired  goals;  however, 
to  ensure  maximum  performance  and  response,  these  limits  may  be 
exceeded  temporarily  during  transients.  This  limit  philosophy  is 
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incorporated  into  the  control  logic  to  provide  the  best  overall 
system  performance. 

2.5.2  Steady  State  Accuracy  Requirements 

Control  accuracy  is  defined  as  the  difference  between  a 
desired  operating  point  and  the  actual  running  condition  in  the 
presence  of  actuator  uncertainties,  disturbances,  engine  varia¬ 
tions  and  degradation.  Many  control  system  designs  use  lag  net¬ 
works  or  integral  logic  to  achieve  a  high  d.c.  gain  to  meet  the 
specified  accuracy  levels.  Thus,  quantified  error  tolerances  can 
be  used  during  the  design  process  to  evaluate  the  regulator 
design  and  synthesize  appropriate  compensation.  Table  2.5  lists 
the  nominal  error  tolerances  allowed  at  all  flight  and  power 
conditions  for  the  GE-23. 


Table  2.5 

Error  Tolerances  on  Set  Points 


SCHEDULED  VARIABLE 

TOLERANCE 

Fan  RPM  (Military  to  Maximum) 

+0.45 

Fan  RPM  (Submilitary) 

+1.05 

Core  RPM  (Military  to  Maximum) 

♦0.25 

Core  RPM  (Submilitary) 

^0. 55 

iP/P 

♦  4.Q* 

Turbine  Blade  Temperature 

+30°R 

Fan  and  Compressor  Geometry 

+3*  to  +5X  F.S. 

Jet  and  Turbine  Nottles 

♦IS  to  +2*  F.S. 

Variable  8ypass  Geometries 

♦is  to  ♦as  F.S. 

Augmentor  Fuel  Flow 

♦4* 

2.5.3  Failure  Tolerance 

The  control  hardware  described  in  Section  2.3  is  designed 
to  produce  a  high  level  of  failure  tolerance  by  integrating 
hardware  and  software  capabilities.  The  control  requirements 
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for  failure  tolerance  can  be  formalized  for  sensor,  actuator, 
and  electronics  hardware.  Table  2.6  summarizes  the  detection, 
isolation,  and  accommodation  criteria  specified.  These  criteria 
represent  major  considerations  in  the  development  of  the 
controller  structure  described  in  Chapter  IV. 

Table  2.6 

Failure  Tolerance  Specification 


FAILURE  NODE 

PRIMARY  REQUIREMENTS 

ADDITIONAL  GOODNESS  CRITERIA 

•  Single  Sensor  Channel 
(all  types) 

a  Detect,  isolate,  and  indicate 
a  No  effect  on  performance 
a  Real  time,  no  hardware 
•  Allow  interaittants 

a  Satisfy  all  requirements 

•  Multiple  Sensor 
Channels 

•  Detect,  isolate,  indicate 
a  Allow  intermittants 

a  Maximum  insensitivity 
•  Conservative  control  action 

a  Actuator  Failures 
(single/multiple) 

a  Protect  Engine 
a  Smooth  transfer  to  backup 

•  Maximum  insensitivity 

•  Detect,  isolate,  indicate 

•  Exceed  backup  control  capability 

•  Processor  Hardware 

a  Detect 

a  Isolation,  Indication 

a  Transfer  to  backup 

•  Reinitialization 

2.6  SUMMARY 

The  requirements  governing  the  development  of  a  multivariable 
controller  for  the  GE-23  engine  are  summarized  in  the  following: 

The  general  function  of  the  digital  electronic  controller 
is  to  provide  acceptable  modulation  of  the  actuators  to  operate 
the  engine  safely  and  efficiently.  Specific  functional  require¬ 
ments  include: 

(1)  Sensor  signal  conditioning  to  eliminate  noise  and 
provide  reconstruction  of  failed  channels. 
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(2)  Command  input  interfacing  to  provide  external 
access  to  power  level  inputs,  inlet  commands, 
and  direct  digital  access  for  executive  inte¬ 
grated  airframe  commands. 

(3)  Regulation  of  outputs  to  operate  servo  actuators 
smoothly  and  stably. 

(4)  Failure  detection  and  accommodation  logic  to 
provide  fail  operate  and  fail  safe  outputs. 

(5)  Condition  indication  inputs  to  maintenance, 
diagnostics,  and  display  processors  via  direct 
digital  links. 

All  of  these  must  support  the  primary  controller  functions  providing 
smooth  engine  operation  in  transient  and  steady  state  regimes 
without  exceeding  engine  operating  restrictions. 

This  section  has  quantified  the  hardware  specifications, 
mission  definition  and  performance  requirements  for  the  con¬ 
troller.  In  Chapter  III,  design  procedures  for  multivariable 
control  synthesis  are  presented.  Chapter  IV  then  integrates 
the  design  problem  and  theoretical  techniques  to  form  a  control 
structure  definition  for  the  GE-23  MVC. 


SECTION  III 

MULTIVARIABLE  CONTROL  DESIGN  METHODOLOGY 

3.1  INTRODUCTION 

In  this  chapter,  the  methodology  required  for  the  design  of 
a  controller  for  the  General  Electric  GE23-JTDE  variable  cycle 
engine  is  reviewed.  Emphasis  is  on  the  design  techniques  which 
are  available.  This  chapter  forms  the  foundation  for  the  control 
structure  which  is  proposed  in  Chapter  IV. 

3.2  THE  CONTROL  DESIGN  PROBLEM 

The  primary  function  of  any  controller  is  to  maintain  stable 
and  smooth  operation  of  the  plant  at  acceptable  levels  of  perfor¬ 
mance  in  both  steady-state  and  transient  regimes  without  exceeding 
specified  operating  restrictions.  Acceptability  is  defined  in 
terms  of  the  specific  application.  Typical  measures  are: 

(1)  Steady- state  hangoff  errors 

(2)  Bandwidth 

(3)  Rise  time  and  overshoot 

(4)  Mean  square  levels  of  control  and  state  responses 

In  addition  to  the  classical  regulation  criteria,  modern  con¬ 
trollers  for  complex  dynamic  systems  must  also  include  logic  to 
provide  sensor  signal  conditioning,  command  input  processing,  and 
failure  detection  and  accommodation.  These  functional  requirements 
are  addressed  by  the  synthesis  procedures  in  the  sections  that 
follow. 

The  structure  and  type  of  the  controller  depends  on  the 
nature  and  complexity  of  the  problem  and  on  the  design  require¬ 
ments.  The  controller  may  be  open  loop  (feedforward)  or  closed 


loop  (feedback)  or  a  combination  of  both.  It  may  be  based  on  the 
so-called  classical  and/or  modern  control  design  techniques. 

Historically,  feedback  control  (i.e.,  using  plant  outputs  to 
effect  compensating  inputs)  has  been  the  foundation  of  control 
structures.  Feedback  control  tends  to  require  less  detailed 
response  models  and  is  therefore  particularly  suited  to  simple 
implementations.  In  modern,  digital  control  applications,  the 
overhead  associated  with  more  detailed  response  models  is  not  an 
important  design  factor.  Thus,  modern  multivariable  control 
designs  are  using  more  detailed  open  loop  strategies  in  conjunc¬ 
tion  with  feedback  regulation.  In  this  way,  superior  transient 
behavior  is  produced  without  many  of  the  unfavorable  stability 
and  noise  problems  often  associated  with  feedback  computation 
alone.  Feedforward  design  concepts  are  discussed  below. 

3.3  FEEDFORWARD  AND  TRANSITION  CONTROL 

If  a  perfect  plant  model  was  known,  and  if  there  were  no 
unknown  disturbances,  then  it  would  be  possible  to  control  the 
plant  using  only  feedforward  (open-loop)  inputs.  Furthermore, 
assuming  the  plant  was  stable,  the  feedforward  could  simply  be  a 
change  in  the  control  set  point  which  would  produce  the  desired 
steady-state  condition.  However,  using  this  strategy,  transition 
from  the  current  state  to  that  commanded  would  occur  with  the 
dynamics  characteristic  of  the  open-loop  plant.  There  would  be 
no  direct  control  over  the  maximum  output  excursions. 

More  control  over  the  response  dynamics  and  maximum  excur¬ 
sions  of  plant  variables  is  possible  if  the  input  is  more  sophis¬ 
ticated.  One  example  is  an  input  time  history  that  requires  the 
state  variables  to  ramp  (at  specified  rates)  from  their  current 
values  to  those  requested  at  the  new  steady  state.  Specifically, 
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if  the  system  is  linear  and  the  vector  of  ramp  rates  is  x  *  R, 
then  from  the  equations  of  motion 

x  =  Fx  +  Gu  (3.1) 

one  can  find  the  control  and  state  time  histories  during  the  tran¬ 
sient.  Defining  an  acceleration  steady  state  as  x  =  0  yields 

x  =  0  =  Fx  +  Gu  (3.2) 

or  0  =  FR  +  Gu  (3.3) 

If  the  number  of  controls  (m)  equals  the  number  of  states  (n) ,  then 
G'1  exists  and 

u  =  -G_1FR.  (3.4) 

If  m>n,  then  m-n  additional  constraints  must  be  imposed  on  the 
controls  to  yield  a  unique  solution  for  u.  If  m<n,  then  m-n  con¬ 
straints  must  be  removed  from  the  state  vector.  Using  the  example 
with  m=n  yields 


u(t)  »  u(o)  + 

w 

/ 

0 

udt 

®  u(o)  -G  ^FR  t 

(3.5) 

x(t)  °  x(o)  + 

t 

/ 

xdt 

»  x(o)  +  Rt. 

(3.6) 

o 


Another  more  sophisticated  approach  calculates  an  input  time 
history  that  is  optimal  with  respect  to  a  specified  performance 
index.  Current  techniques  for  trajectory  optimization  [7] 
utilize  an  engine  simulation  of  varying  complexity.  A  cost 
functional  is  derived  which  represents  desirable  trajectory 
attributes  (e.g.  maximum  thrust  reponse  or  minimum  turbine  inlet 
temperature  rates),  and  a  group  of  constraints  in  the  trajectory 
are  formulated  to  respect  stability  limits  and  other  physical 
operating  constraints.  A  trial  trajectory  is  calculated 
which  may  or  may  not  satisfy  the  constraints.  Standard  function 


optimization  procedures  can  be  used  to  derive  feasible  and  optimal 
paths. 

The  problems  with  this  type  of  procedure  are  threefold.  First, 
engine  transient  simulations,  especially  for  gross  transitions,  may 
poorly  represent  engine  behavior;  second,  implementation  of  exact 
trajectories  is  inevitably  impractical;  and  third,  there  is  no 
convenient  control  law  formulation  of  the  optimal  solution. 

An  interesting  technique  [8]  has  been  proposed  which 
could  yield  substantially  improved  implementation  properties  which 
make  the  scheme  compatible  with  real-time,  on-board  application. 

A  fairly  general  dynamic  optimization  problem  can  be  formulated 
where  the  performance  index,  system  equations,  and  path  constaints 
are  modeled  as  linear  functions  of  the  control.  The  solution  is  a 
bang-bang  control  law.  The  vector  control  history  u(t)  which 
minimizes  the  performance  functional 

J  =  *(x(tf),tf)  ♦  /tf  (a(x,t)  *  bT(x,t)u(t))dt  (3.7) 

*o 

subject  to  the  vector  system  differential  equations 

x  =»  f ( x , t )  +  g(x,t)u(t)  (3.8) 

and  of  path  inequality  constraints 

e^x.t)  ♦  d]  (x,t)u(t)*0  i  ■  1,2, ...,q  (5.9) 

must  be  calculated. 

The  initial  state  and  time  are  assumed  to  be  specified,  i.c., 
x(tQ)  -  x0  (3.10) 

while  the  terminal  state  and  time  are  subject  to  the  p  terminal 
constraints  (p  <  n  ♦  1) 

*i(x(tf) ,tf)  =  0  i  -  1,2 . p  (5.11) 
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The  functions  <J»,  a  and  ci  are  scalar  ftmctions  of  x  and  t 
while  b  and  d^  are  vector  functions  of  dimension  (r  x  1).  The 
terminal  time  t^  may  be  either  fixed  or  free. 

The  necessary  conditions  [8]  for  the  optimal  solution 
can  be  derived  as  follows: 


=  ft  ’  §|  u(t)  -  x1  |£  -  XT  !&  u(t) 


q  f  3ci  ad’*'" 

•  *  "i  TS  u'l) 


(3.12) 


bT  *  *Tg  * 


u-d!  .  0 

11 


(3.13) 


t  Oj  B  0  if  ci  ♦  dTj  <  o 


(3.14) 


The  switching  function  in  liq.  (3.13)  can  be  calculated  and 
perhaps  stored  as  a  feedback  control  law  in  a  modified  bang*bang 
form.  Another  approach  which  mimics  the  character  of  the  above 
optimal  trajectory,  while  avoiding  the  required  computation  is 
discussed  in  Section  4.1.5. 

As  indicated  in  the  beginning  of  this  section,  a  feedforward 
strategy  alone  (even  "optimal*')  may  not  be  adequate  to  control  a 
plant  due  to  modeling  errors  and  unknown  disturbances.  Regulator 
logic  is  usually  required.  In  general  terms,  the  objective  of 
regulator  design  is  to  determine  a  set  of  feedback  gains  which 
provide  a  system  with  desired  response  characteristics.  Such 
desired  characteristics  include  stability,  frequency,  damping, 
decoupling,  and  minimum  error  in  following  a  specific  command. 
Classical  and  modern  control  techniques  for  designing  regulators 
are  discussed  below, 
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3.4  CLASSICAL  CONTROL  SYNTHESIS 


Classical  design  techniques  are  we!l  known  and  will  not  he 
reviewed  in  detail  here.  When  applicable,  these  techniques  are 
very  powerful  and  can  lead  to  easily  designed  and  effective 
regulators.  In  fact,  over  40  years  of  successful  control  system 
designs  have  been  achieved  by  what  is  now  called  classical  synthe¬ 
sis.  This  approach  is  based  on  analysis  of  the  transfer  function 
representation  of  the  system  (e.g,,  the  frequency  domain  relation¬ 
ship  between  input  and  output). 

Examples  of  design  techniques  which  have  been  used  are  the 
root  locus  method,  the  Nyquist  method,  and  the  Bode  method  (a 
derivative  of  the  Nyquist  method).  The  basic  design  approach  is 
to  analyze  the  transfer  function  with  respect  to  desired  system 
characteristics  (e.g.,  gain  and  phase  margin,  transient  response) 
and  introduce  compensating  lead  or  lag  filters  to  modify  the  re¬ 
sponse  to  that  desired. 

The  compensators  can  provide  lead  or  lag  and  may  introduce  a 
pure  integral  effeet  to  eliminate  steady-state  errors  from  step 
inputs.  Sufficient  stability  margin  is  provided  such  that  small 
errors  in  the  design  parameter  values  do  not  afieet  the  stability 
of  the  system.  All  modes  can  be  made  adequately  fast  and  well 
damped  within  design  requirements.  These  methods  were  originally 
devised  for  single- input,  single-output  systems  and  are  difficult 
to  use  with  multi-input,  multi-output  systems.  A  summary  of  the 
advantages  and  disadvantages  of  these  techniques  is  given  in 
Table  5.1. 

5.S  MODERN  CONTROL  SYNTHESIS  METHODS 

Modern  control  design  methods  have  their  basis  in  the  state 
space  representation  of  dynamic  systems.  The  state  space  approach 
leads  naturally  to  the  description  of  both  linear  and  nonlinear 
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Table  3.1 

Classical  Control:  Advantages  and  Disadvantages 


ADVANTAGES 

DISADVANTAGES 

CLASSICAL 

CONTROL 

(1)  Provides  basis  for  engineering 
design  decisions  and  evaluation 
(for* low-order  systems). 

(2)  Inherent  use  of  stability  margin 
(e.g.,  gain  and  phase  margins) 
to  reduce  controlled  sensitivity 
to  uncertainties  in  system 
parameters. 

(3)  Best  for  single-input/single¬ 
output  systems. 

(1)  Much  trial-and-error  in  design, 
particularly  for  large  systems 
such  as  engines. 

(2)  Frequently  leads  to  controls 
fighting  each  other. 

(3)  Difficult  to  apply  to  multi- 
input/multi-output  systems. 

systems  in  time  domain.  The  control  inputs  may  be  chosen  to  mini¬ 
mize  a  certain  criterion  function  which  is  a  quantitative  measure 
of  the  penalty  for  undesired  state  response  and  the  control  input 
effort  required  (quadratic  synthesis  methods).  Alternately,  the 
control  may  be  selecied  by  specifying  the  desired  transient  response 
characteristics  of  a  system,  and  using  algebraic  synthesis  tech¬ 
niques  to  solve  for  the  required  gains  (e.g.,  pole  placement  meth¬ 
ods,  modal  control  methods,  etc.). 

A  major  control  theory  development  occurred  during  the  late 
1950's  which  was  fostered  by  high-speed  computers  and  an  advanced 
technology  space  program.  Now  known  under  various  names  (quad¬ 
ratic  synthesis,  linear-quadratic  Gaussian  synthesis,  modern  opti¬ 
mal  control  theory,  linear  quadratic  regulator  theory),  this  devel¬ 
opment  has  been  the  subject  of  intensive  theoretical  research.  It 
is  known  that  optimal  control  has  great  potential  since  it  is 
inherently  a  multi- input/multi -output  linear  system  control  design 
method  also  applicable  to  nonlinear  systems.  Much  of  the  demonstra¬ 
tion  of  the  versatility  of  the  method,  however,  has  been  limited  to 
simulation  studies,  and  the  applications  through  actual  implementa¬ 
tion  are  not  nearly  as  extensive.  Characteristics  of  quadratic 
synthesis  are: 
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(1)  Quadratic  synthesis  techniques  design  on  a  different 
objective  than  classical  techniques)  Instead  of 
attempting  to  obtain  a  specific  transient  response  fas 
do  classical  methods) ,  LQR  methods  minimize  the  control 
energy  required  to  keep  the  mean  square  response  of  the 
system  as  small  as  possible.  The  design  parameters  are 
weightings  or  penalties  on  deviations  of  states  and 
controls  as  well  as  the  choice  of  augmented  states  and 
controls).  Desirable  transient  response  is  obtained 
indirectly. 

(2)  Quadratic  synthesis  techniques  are  highly  automatic 
once  the  performance  index  and  the  design  model  have 
been  selected. 

(3)  Quadratic  controls  can  produce  simple  control  systems  if 
the  design  process  is  conducted  with  an  integrated  under- 
standing  of  the  system  physics.  This  is  because  many  of 
the  simplifications  to  an  "optimal  controller"  can  be 
based  on  certain  types  of  analysis  of  the  feedback  struc¬ 
ture  and  relating  this  structure  to  its  actual  effect  on 
the  system. 

Modern  optimal  control  theory  can  accommodate  nonlinear  systems 
and  criterion  functions.  However,  the  control  laws  thus  generated 
may  be  very  complicated  (e.g.,  non-state- variable  feedback).  An 
alternate  method  is  to  schedule  a  set  of  linear  (state  variable 
feedback)  control  laws  throughout  the  operating  region  of  the 
system.  Each  control  law  is  designed  using  a  linearized  model  of 
the  system  in  its  region  of  effectiveness.  The  piecewise  linear 
control  law  formed  in  this  fashion  is  generally  simpler  to  calculate 
and  easier  to  implement  than  the  corresponding  nonlinear  control 
law;  yet  it  sacrifices  little  performance. 

The  advantages  and  disadvantages  of  modern  control  synthesis 
methods  are  summarized  in  Table  3.2. 


3.5.1  State  Variable  Feedback 


Most  common  applications  of  LQR  theory  involve  calculation  of 
a  feedback  gain  matrix  on  the  system  state  variables.  These  gains 
minimize  a  criterie  function  (or  performance  index),  J,  of  the 
form 


J  *  <M*(tf),  ♦  /  l(x(t)  ,u(t)  ,t]  dt 

'o 


(3.15) 
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Table  3.2 

Modern  Control:  Advantages  and  Disadvantages 


ADVANTAGES 

DISADVANTAGES 

MODERN 

CONTROL 

(1)  Leads  to  automatic  synthesis  of 
controls  once  performance  index 
and  model  have  been  determined. 

(2)  Designs  are  inherently  minimally 
sensitive  if  proper  P.l.  is 
chosen. 

(3)  Automatically  includes  optimum 
control  and  state  cross¬ 
coupling. 

(1)  Places  more  emphasis  on 
correctness  of  model  used 
for  design  (e.g.,  if  model 
is  in  error,  controlled 
system  could  be  unstable) 
and  relevance  of  performance 
index. 

(2)  Most  applications  tend  to  be 
naive  because  too  much  confi¬ 
dence  placed  on  automatic  LQR 
methods  to  replace  knowledge 

of  system  and  modeling  accuracy. 

where  <j>  represents  a  penalty  on  the  terminal  error  and  L 
represents  a  penalty  on  trajectory  errors.  Several  choices  of  <f> 
and  L  are  possible.  A  common  choice  is  state  and  control  weight¬ 
ing  of  the  form 

J  =  /°f  (xT(t)A(t)x(t)  +  uT(t)B(t)u(t) ]  dt  (3.16) 

*o 

where  Aft)  is  a  weighting  matrix  on 
a  weighting  o"  control  effort. 

Assuming  the  state  equations  are 

x  *  Px  ♦  Gu 
y  -  Hx  ♦  Du 

the  control  law  which  minimizes  Equation  (3,16)  is  given  by 

u(t)  »  *C(t) x(t)  (3.19) 

where 

C(t)  *  IB(t)]‘l  CT(t)S(t) 


C3.17) 

(3.18) 


state  errors  and  B(t)  is 

in  matrix  form 


(3.20) 


« 

and  S(t)  is  the  solution  of  a  matrix  Riccati  equation 

S  =  -SF  -  FTS  +  SGB*1GTS  -A  (3.21) 

A  constant  feedback  gain  matrix  that  yields  "nearly  optimal" 
control  is  given  by  the  steady  state  solution  to  Equation  (3.21) 
(setting  S  =  0).  Computationally  efficient  algorithms  exist 
for  solving  this  problem.  An  eigenvector  decomposition  tech¬ 
nique  is  described  in  Ref.  9.  The  result  is  a  control  law  of 
the  form 

u(t)  *  -  C  x (t )  (3.22) 

3.5.2  State  Estimation  and  Sensor  Compensation 

In  the  preceding  discussion  it  was  assumed  that  the  controller 
had  perfect  knowledge  of  the  state.  A  model  reduction  procedure 
can  be  used  to  allow  the  control  designer  some  freedom  to  choose 
conveniently  measurable  quantities  as  state  variables.  However, 
this  is  not  always  possible.  Models  containing  unmeasured  states 
must  sometimes  be  formulated  to  describe  system  response  accurately. 
Measurement  errors,  model  errors  and  sensor  dynamics  all  tend  to 
make  the  knowledge  of  the  state  only  approximate.  The  dependence 
of  the  regulator  control  law  on  measurement  uncertainties  must  be 
analyzed  so  that  the  impact  on  performance  can  be  evaluated. 

The  regulator  gains  are  assumed  independent  of  the  uncertainty 
in  the  knowledge  of  the  state.  This  assumption  allows  the  synthe¬ 
sis  procedure  to  maintain  its  systematic  simplicity  and  separates 
the  compensation  of  the  transduced  signals  from  the  control  de¬ 
sign.  The  justification  of  this  assumption  is  discussed  below. 

For  the  general,  1 inear  stochastic  control  problem,  the  system 
follows: 


x  «  Fx  ♦  Gu  ♦  w 


(3.23) 
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y  *  Hx  +  v 


(3.24) 


In  this  case,  w  and  v  are  zero  mean,  Gaussian  random  processes 
with  covariances  defined  below: 

E[w(T)wT(t)]  =  Q6(t-x)  (3.25) 

E[v(x)vT(t)]  =  R6 (t -x)  (3.26) 

The  performance  index  for  the  problem  remains  unchanged  ex¬ 
cept  that  the  mean  value  of  the  cost  is  minimized. 

The  solution  to  the  problem  follows  directly  from  a  change  of 
coordinates.  The  maximum  likelihood  estimate  of  the  state  given 
all  the  measurements  to  time  t  is  the  function  x(t)  that  makes 
the  residual  process  white,  i.e., 

v ( t )  =  y(t)  -  Hx(t)  (3.27) 

E[v(t) vT(t) ]  =  M6(t-x)  (3.28) 


It  can  be  shown  that  there  is  an  (nxp)  matrix  function  of  time, 
K(t) ,  which  will  produce  the  maximum  likelihood  estimate  as 
follows: 


x  B  Fx  +Gu  +  K(t) v(t) 


(3.29) 


Now,  the  problem  can  be  rewritten  with  new  variables  as  follows: 


e  «  x  -  x 


J  -  ifi  |[xT|uT)A>N  i|+E(eTAeJ 

2  I  »  (n1;  bJ  LuJ  j 


(3.30) 


(3.31) 


The  error  term,  e,  follows  the  linear  system  relationship, 


e  ■  (F  -  KH)e  ♦  Ky  ♦  w 


(3.32) 
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and  is  independent  of ~ the  control.  The  regulator  problem  minimizes 
the  first  term  of  Eq.  (3.31).  The  solution  is  derived  [10] 
as  follows: 

u(t)  =  -  Cx(t)  (3.33) 

where  the  matrix  C  is  determined  under  the  assumptions  of  per¬ 
fect  state  knowledge. 

For  nonlinear  systems  and  for  systems  without  white  Gaussian 
disturbances,  this  separation  principle  may  not  hold.  However,  it 
is  often  approximately  true  and  utilization  of  separate  design 
procedures  for  control  and  estimation  yields  satisfactory  closed- 
loop  systems. 

This  philosophy  has  been  used  for  many  years  in  classical 
design.  The  goal  of  sensor  compensation  is  to  remove  the  random 
error  from  the  signal  and  produce  the  best  estimate  of  the  mea¬ 
sured  quantity  for  control. 

There  are  two  criteria  used  in  sensor  compensation.  If  a 
sensor  has  a  very  high  bandwidth,  spurious,  high  frequency  input 
may  enter  the  feedback  loop  and  cause  unacceptable  behavior.  The 
standard  remedy  is  a  low  pass  filter  at  a  bandwidth  higher  than 
the  closed-loop  frequency. 

For  most  digital  control  systems,  it  is  standard  practice  to 
filter  all_  input  signals  to  eliminate  aliasing.  Each  transduced 
input  will  have  an  ample  amount  of  noise  at  all  frequencies  due 
to  line  pickup,  power  supply  fluctuation,  EMI,  etc.  Sampling  into 
this  noise  can  often  lead  to  difficult  problems.  Higher  order 
filters  with  steeper  attenuation  characteristics  are  sometimes 
necessary. 

Sensor  compensation  must  also  improve  poor  dynamic  response 
relative  to  the  closed- loop  dynamics.  Many  transducers  can  be 
modeled  by  first  order  lags.  If  the  bandwidth  of  these  sensors  is 
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significantly  lower  than  the  desired  closed-loop  bandwidth,  poor 
transient  behavior  can  be  anticipated  due  to  the  phase  lag  intro¬ 
duced  by  the  instrument. 

Classical  compensation  philosophy  can  also  be  applied  [2] 
to  multivariable  systems.  With  many  measurements,  it  is  pos¬ 
sible  to  choose  the  high  bandwidth  variables  to  regulate  the  tran¬ 
sient  response  of  the  system.  Slower,  high  accuracy  sensors  can 
be  utilized  for  trim  or  bias  control.  The  design  proceeds  by  the 
appropriate  choice  of  model,  quadratic  weights  and,  finally,  fil¬ 
ter  compensation. 

It  is  possible  to  implement  a  full  order  filter  directly. 

This  increases  the  complexity  of  the  control  law  from  simple  pro¬ 
portional  feedback  to  an  nth  order  dynamic  compensation.  An  exam¬ 
ination  of  the  optimal  filter  roots  for  a  group  of  measurements 
typically  shows  that  some  filter  time  constants  will  be  well  above 
the  closed-loop  bandwidth.  Single-loop  prefiltering  can  be  used 
to  limit  noise  transmittance  in  these  inputs.  The  measured  vari¬ 
ables  will  naturally  decompose  into  a  set  of  high  bandwidth  signals 
which  are  necessary  to  the  control  function  but  which  do  not  have 
sufficient  response  for  use  directly  in  the  control  law  [2]. 

3.5.3  Output  Feedback 

In  many  applications,  the  quantities  that  are  to  be  regulated 
are  system  output  variables  and  not  system  state  variables.  Using 
state  variable  feedback  in  these  applications  (outputs  weighted  in 
the  performance  index)  generally  involves  estimation  of  unmeasured 
states  (e.g.,  the  Kalman  filter)  and  math. metical  modeling  of  the 
system’s  output  equations  (e.g.,  y  -  Hx  ♦  Du).  If  the  mathematical 
model  of  the  plant  were  known  exactly,  this  procedure  would  produce 
an  acceptable  control  law.  However,  modeling  errors  are  generally 
present.  Consequently,  additional  integral  or  trim  control  is 
often  necessary  to  insure  satisfaction  of  performance  requirements. 
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These  difficulties  may  often  be  avoided  by  designing  a  con¬ 
troller  that  feeds  back  output  errors  directly  [11].  One  dis¬ 
advantage  to  this,  however,  is  that  the  resulting  gain  matrices 
cannot  be  found  as  solutions  to  a  matrix  Riccati  equation.  They 
must  be  found  by  solving  a  set  of  nonlinear  equations. 

The  output  feedback  problem  can  be  formulated  as  follows. 

Given  the  set  of  plant  equations 

x  =  Fx  +  Gu 

y  =  Hx  +  Du 

find  the  matrix  C  of  the  form 

u  =  Cy  (3.36) 

which  minimizes  the  performance  index 

J  =  1/2  f*  [*TAxx  +  uTBuu  +  yTAy  y]  dt  (3.37) 

or  equivalently 

J  =  1/2  f"  [xTiuT]  dt  (3.38) 

0  N  ;B  u 

This  performance  index  is  clearly  a  function  of  the  gain  matrix, 
C.  One  technique  of  finding  its  minimum  (with  respect  to  C)  in¬ 
volves  the  use  of  a  gradient  search  procedure.  One  efficient  proce¬ 
dure  is  the  modified  Newton  Raphson  with  derivative  technique  [12]. 

It  requires  the  gradient  of  J  with  respect  to  C. 

The  gradient  of  J  with  respect  to  C  is  developed  in  part  in 
Ref.  13.  Summarizing  briefly,  let 

y*  ■  Hx  ■  y  -Du  (3.39) 


(3.34) 

(3.35) 
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and 


with 


u  =  C  y* 


C  =  C(I-DC) 


Then,  from  Refs.  13  and  14 


where 


^  =  (BCH  +  NT)  LHT  +  GTK  LHT 
dC 


F*  =  F  +  GCH 
0  =  LF*T  +  F*L  +  X 


C  3 . 40) 


(3.41) 


(3.42) 


(3.43) 

(3.44) 


0  =  KF*  +  F*TK  +  A  +  HTCTB£H  +  NG'H  +  HTCTNT  (3.45) 


XQ  *  E(x(0)  xT(0)  ] 


(3.46) 


Note  that,  in  this  formulation 


J  =  1/2  trace  (K  X) 


(3.47) 


Thus,  the  gradient  of  J  with  respect  to  C  is  provided. 

The  gradient  search  may  therefore  be  used  to  find  a  matrix 
C  (u  «  C  y*)  which  minimizes  J.  The  resulting  C  may  be  transformed 
to  C  (u  «  C  y)  using  Equation  (3.41). 

One  difficulty  remains.  The  resulting  feedback  gain  matrix, 

C,  is,  in  general,  full.  That  is,  there  is  a  feedback  path  from 
every  output  to  every  control.  Often,  unimportant  or  ineffective 
error  to  actuator  paths  can  be  identified  and  eliminated  -  yielding 
significant  reduction  in  the  number  of  feedback  paths  which  must 
implemented.  Optimization  of  this  resulting  fixed-structure 
gain  matrix  requires  an  extension  to  existing  theory,  however. 
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Specifically,  the  gradient  of  J  with  respect  to  C  (instead  of  C) 
must  be  found. 

The  derivation  of  ~  is  provided  in  Appendix  A.  The  result 
is 

~  »  (I  +  C8D)T  ~  £T  (3.48) 

dC  dC 

where 

S  -  (I-DC)'1  (3.49) 


With  —  available,  the  gradient  search  procedure  can  be  used 

to  optimize  specific  elements  in  the  C  matrix  with  the  remaining 
elements  fixed  (e.g.,  at  zero).  A  block  diagram  of  the  procedure 
is  presented  in  Figure  3.1. 

3.6  MODEL  GENERATION 

The  design  procedures  discussed  above  rely  heavily  on  an 
accurate  mathematical  model  of  the  system.  Consequently,  a  criti¬ 
cal  step  in  the  control  design  procedure  is  the  generation  of 
tractable  design  models.  Linear  models  are  one  class  of  such 
design  models,  and  as  will  be  shown,  can  be  effectively  integrated 
to  provide  a  nonlinear  control  function. 

A  dynamic  description  of  the  response  of  important  variables 
near  a  trim  condition  it  the  goal  of  linear  model  generation.  From 
these,  closed-loop  controls  can  be  developed.  Typically,  complex 
simulations  are  developed  before  the  actual  system  is  built.  The 
linear  models  used  in  control  design  can  be  produced  from  these 
large  simulations.  Alternately,  data  from  prototype  tests  can  be 
processed  to  develop  models  which  are  accurate  representations  of 
the  behavior  useful  in  control  design  and  validation  of  the  non¬ 
linear  simulation. 
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Figure  3.1  Block  Diagram  of  Output  Regulator  Design 
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An  equilibrium  point  is  determined  for  a  given  set  of  control 
inputs,  which  can  be  written  as  follows: 

0  »  f(x,  u,  h  (3.50) 

where  f(x,  u,  8)  are  the  detailed  nonlinear  equations  governing 
system  operation.  The  states,  6x,  are  a  group  of  independent 
dynamic  variables.  The  controls,  6u,  determine  the  inputs.  The 
linearization  point  is  determined  by  the  ambient  conditions,  0. 

Expanding  Eq.  (3.50)  about  trim  yields  the  following: 

6x  *  fx<$x  +  fu6u  +  g(6x,6u)  (3.51) 

where  6x  are  the  perturbation  state  rates,  and  g  is  at  least  second- 
order  small  in  the  perturbations. 

Standard  offset  derivative  generation  procedures  use  Eq.  (3.51) 
directly.  The  selected  control  and  ambient  inputs  are  applied  and 
an  equilibrium  point  is  reached.  All  dynamic  integrations  are  then 
held.  States  and  controls  are  perturbed  in  order  and  the  inputs  to 
the  integrators  calculated.  The  calculation  is  written  as  follows: 


(fx) ij  *  (fiXj/taj) 

i,j,k-l,.. . ,n 

(3.52) 

kj»j 

(Vh  *  C«ii/«u4) 

t»l,. . . ,m 

(3.53) 

dXjjO 


A  modification  to  this  method  can  be  used  to  improve  the  low 
frequency  accuracy  of  the  models  [15} .  First,  the  simulation 
is  trimmed  at  the  selected  flight  condition.  States  are 
perturbed  in  order  and  the  dynamics  matrix  calculated  as  in  Eq. 
(3.52).  Controls  are  perturbed  in  order  and  the  resulting  steady- 
state  trim  point  calculated.  The  models  are  formed  as  follows: 


<f*hj  ’  (W 


4V° 


i,j ,k-l, . . . ,n 
ki»j 


(3.54) 
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In  the  standard  method,  the  steady-state  response  gain  from 
the  control  to  the  state  perturbations  is  given  by  Eq.  (3. 56) 

[derived  from  Eqs.  (3.S2)  and  (3.53)].  This  gain  matrix  depends  on 
the  independently  calculated  quantities  lfv)  and  (f  ).  In  practice, 
these  independent  calculations  are  susceptible  to  numerical  inac¬ 
curacies  which  degrade  the  low  frequency  response  of  the  associated 
linear  models.  However,  the  accuracy  at  all  frequencies  in  model¬ 
ing  control  inputs  is  presumably  good.  The  forced  match  procedure 
uses  the  steady  response  as  the  defining  relationship  (Eq.  (3,57)] 
in  the  calculation  of  the  control  distribution  matrix.  This  forces 
the  steady- state  response  of  the  nonlinear  simulation  and  the 
linear  model  to  coincide  for  at  least  one  set  of  inputs.  In  general, 
the  low  frequency  response  of  these  models  is  quite  accurate. 

Both  procedures  are  computationally  efficient  and  yield  models 
suitable  for  closed-loop  control  design.  Also,  both  yield  linear 
models  of  the  fora 

x  *  Fx  ♦  Gu  (3.58) 


y  <*  Hx  ♦  Du 

directly.  Specifically,  the  elements  in  Ff  G,  H  and  D  are 

r(ij)  - 


(3.59) 

(3.60) 
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G(M) 


(3.61) 


(Vii 

■  or 


and,  though  not  discussed,  but  derived  by  similar  means 


H('k,i)  = 

^8xHi 

(3.62) 

D(k,n  = 

^8uHi 

(3.63) 

where  there  are  p  measurements  given  by 


and 


y  =  g(x,  u,  0) 


(3.64) 


(3.6S) 


An  alternate  strategy  is  possible  in  which  P,  G»  H  and  D  are 
determined  indirectly.  This  technique  requires  the  solution  of 
a  set  of  linear  equations  (and  the  associated  numerical  difficul¬ 
ties)*  but  may  be  more  convenient  for  some  forms  of  simulation 
programs.  Specifically,  the  incremental  changes  in  the  states  and 
the  outputs  are  found  due  to  perturbations  in  the  state  derivatives 
and  controls.  Hence  the  following  matrix  is  identified; 


•  a 

X 

* 

An  :  Ai  z " 

s  •  a*  a  m*  >4k  4*  «* 

x  ' 

K  m  m 

y 

L  -j 

.  A21  !  A22  . 

U 

(3.66) 


where 


An(ij: 

-  ■  (Sf) 

6V° 

i=l , . . .  ,n 
j -1 , . . .  ,n 

(3.67) 

A1 2  ^ 1  ^ - 

>  ■  Q) 

s*r° 

(3.68) 

^2i  (i > j . 

1  •  (5j> 

«ik=° 

(3.69) 

A22  ^  ■ 

>■© 

6V° 

(3.70) 

F,  G,  H  and  D  may  then  be  found  as 

F  =  (3.71) 

G  =  '  A11  A12  (3.72) 

H  =  A21  AxJ  (3.73) 

D  *  A22  '  A21  A11  A12*  (3.74) 

Generating  perturbation  size  is  important  in  all  techniques 
because  of  inherent  nonlinearities  in  the  behavior  and  numerical 
accuracy  in  both  the  hybrid  and  digital  environment. 

Often,  linear  models  generated  numerically  do  not  contain  the 
most  convenient  parameterization  of  the  dynamics  and,  indeed,  con¬ 
tain  far  too  complex  a  description  to  be  practically  utilized  for 
design.  Some  method  is  required  to  analyze  the  dynamical  models  and 
establish  simpler  systems  which  include  only  elements  important  to 
the  desired  control  function.  Without  such  simplification,  appli¬ 
cation  of  design  procedures  can  result  in  highly  complex  and  param¬ 
eter  sensitive  controlled  systems.  Utilizing  the  reduced  order 
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models,  regulator  synthesis  procedures  become  physicaly  intuitive 
and  far  less  sensitive  to  parameter  variation. 

3.6.1  Modal  Decomposition 

A  modal  decomposition  provides  the  framework  for  reducing 
arbitrary  linear  models  to  design  models  containing  the  appro¬ 
priate  parameterization  for  the  control  function.  The  procedures 
are  well  known  [2]  and  their  application  to  model  reduction  is 
described  below. 

The  linear  Eqs .  (3.58)  and  (3.59)  can  be  transformed  to  block 
diagonal  form  assuming  the  n*n  dynamics  matrix,  F,  has  no  repeated 
eigenvalues : 


x  =  Tz  (3.75) 

z  »  Az  +  =  u  (3.76) 

y  *  HTz  ♦  Du  (3.77) 

where  A  is  an  n*n  block. diagonal  matrix,  T  is  an  nxn  matrix  com¬ 
posed  of  the  column  eigenvectors  of  F,  Z  is  an  nxl  modal  coordinate 
vector,  and  2  is  the  n*ro  modal  control  distribution  matrix.  Also, 

FT  -  TA  (3.78) 


The  system  of  Eqs.  (3. 7S)  through  (3.77)  can  be  partitioned 
into  a  set  of  q  states  and  q  eigenvalues  (time  constants)  and  n-q 
states  and  eigenvalues  as  follows: 


~*l~m 

;Tu- 

- 

Y22  '  T22 
m*  ni 

'  rq 

1  M 

! _ 1 

(3.79) 


(3.80) 
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where  x^  and  z^  are  qxl  vectors  partitioning  the  states  and  modes 
and  x2  and  z2  are  (n-q)xl  vectors  partitioning  the  remaining  states 
and  modes . 

If  the  following  equilibrium  relationship  is  approximately 
true  (within  the  time  frame  of  control  interest) , 

z2  %  0  (3.81) 

then  the  following  reduction  can  be  made 

*1  “  Frxl  +  Gru  (3.82) 

where  x^  is  now  the  qxl  state  vector,  Fr  is  the  qxq 
dynamics  matrix,  and  Gr  is  the  qxm  control  distribution 
matrix.  Also, 


[*?: 

5 

X.  + 

~D*~ 

v: 

i _ 

1 

J  r_ 

(3.83) 


where  x2  is  treated  as  an  additional  (n-q)  x  1  output  vector  with 
a  (n-q)  xq  state  distribution  matrix  H*  and  a  (n-q)  xm  control 
distribution  matrix  D*.  The  original  output  distribution 
matrices,  H  and  D,  are  modified  to  Hr  and  Df,  respectively. 
The  formulas  for  these  matrices  are  shown  below  in  terms  of  the 
partitioned  modal  decomposition: 


Fr  ■  TnAi  Tn 


Gr  •  Tll»lTiiT12»ilW 


«*  -  T2iTii 

°*  ■  tT21TliT12"122}AZls2 


(3.84) 

(3.85) 

(3.86) 


(3.87) 


(3.88) 


H  *  H,  +  H,H* 

r  l  ^ 

Dr  =  D  +  H2D*  (3.89) 

Care  must  be  exercised  in  the  partitioning  of  the  system  to 
assure  that  states  retained  span  the  space  of  the  reduced  modal 
system.  In  other  words,  must  be  invertable. 

3.6.2  Second  Stage  Reduction 

The  control  structure  proposed  in  Chapter  IV  requires  an  ac¬ 
curate  mathematical  model  of  the  engine  throughout  its  operating 
flight  envelope.  It  is  highly  desirable  to  provide  this  model  in  a 
simple  form  to  minimize  computer  storage  and  computation  require¬ 
ments.  A  first  step  in  a  model  reduction  procedure  has  been  de¬ 
scribed  above.  It  involves  the  reduction  in  order  of  the  system 
equations  to  include  only  those  modes  which  are  important  for 
control  purposes.  Outlined  below  is  a  procedure  for  further 
simplifying  the  reduced  order  system  (system  order  remains  un¬ 
changed).  The  procedure  is  based  on  eliminating  (setting  to  zero) 
those  elements  in  the  system  dynamics  matrix,  F,  which  have  "little1' 
effect  on  the  system's  state  dynamic  response.  The  result  is  a 

a 

"simplified"  system  dynamics  matrix,  F.  Note  that  concentrating 
on  system  dynamics  is  justifiable  if  the  equations  are  implemented 
in  the  form 

* 

x  «  F  (x  -  g(x,u,0)) 

or  •  «  a 

x  ■  F  (x  -  g(x,u,0)) 

since  the  correct  steady  state  is  guaranteed  by  the  steady  state 
schedules. 

The  procedure  is  comprised  of  three  steps: 

Step  1:  Identify  those  elements  in  F  which  have  little  effect  on 
the  engine's  dynamic  response.  Specifically,  let  JQ  represent 
a  measure  of  the  total  system  state  response. 

^  J  x*A  X  dt  (3.90) 
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with 


Step  3:  Modify  the  non- zero  elements  in  F  to  minimize  the  mean 
square  difference  between  the  initial  condition  responses  of  the 
reduced  system  (x  ®  Fx)  and  the  original  (x  ■  Fx) .  Specifically,  find 

F  »  min  jj  "  \  f  (x-x)1  A2(x-x)  dtj  (3.P4) 

F  0 


where 


x  -  Fx 
l  -  Fx 
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and 


with 


E  x(o)  x(o)^  =  X. 

J 


x(o)  =  x(o) . 


F  is  the  optimum  reduced  F  matrix. 

A  procedure  for  solving  Step  3  is  presented  in  Appendix  B.  A 
block  diagram  of  the  overall  model  reduction  procedure  is  presented 
in  Figure  3.2. 


Figure  3.2  Generation  of  Design  Models 
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so 


The  solution  is 


F  = 


-6.495  0 

0  -20 


Note  that  due  to  the  one-way  coupling  in  the  system,  the  F22  element 
was  not  affected;  only  the  element  was  altered. 

These  results  are  plotted  in  Figure  3.3;  x^(t)  is  the  original 
response,  x^(t)  is  the  response  of  the  reduced  system,  and  x^ft)  is 
the  response  of  the  system  obtained  by  setting  F12  «  0  without  alter¬ 
ing  Fjj,  (i.e.,  =  -10  x^).  These  results  show  that  by  increasing 

the  system  time  constant,  the  reduced  system  can  be  made  to  dupli¬ 
cate  more  closely  the  original. 


Figure  3.3  Comparison  of  Time  Histories  for  Second  Stage 
Reduction  example 


3.7  NONLINEAR  MODEL  DEVELOPMENT 


The  procedures  used  to  derive  linearizations  of  a  complex  non¬ 
linear  simulation  at  selected  flight  points  is  discussed  in  Section 
3.6.  These  linearizations  are  the  basis  for  point  design  of  control 
regulation  logic  discussed  above.  The  simplified  dynamic  models  can 
also  be  used  in  the  controller  to  implement  transition  logic,  filter 
update  equations,  and  fault  tolerant  strategies.  The  primary  require¬ 
ment  for  a  simple  model  equation  arises  from  the  real  time,  micro¬ 
processor  implementation.  For  example,  a  typical  nonlinear  digital 
design  simulation  of  an  advanced  engine  requires  64-128  thousand 
words  of  memory  and  uses  several  thousand  multiplication  operations 
for  each  integration  step.  These  programs  describe  each  component 
using  empirical  operating  maps  and  detailed  aerothermodynamic  equa¬ 
tions  to  derive  the  cycle  balances  existing  during  engine  operation. 
Clearly,  this  type  of  implementation  is  far  too  complex  to  use  in  a 
real  time  microprocessor  environment.  However,  the  program  derived 
linear  models  can  be  used  as  a  starting  point  for  development  of 
real-time,  simplified  nonlinear  equations  which  closely  match  the 
full  engine  simulation  program.  This  method  is  described  below. 

Engine  dynamics  can  be  represented  by  the  nth  order  set  of  nonlinear 
equations: 


X  »  f(x,U,0) 

(3.95) 

y  ■  h(x,u,9) 

(3.96) 

where  m  quantities,  u,  (controls)  and  q  inputs,  0 ,  (ambients)  determine 
the  n  states  and  state  rates,  x  and  x  respectively.  These  variables 
appear  as  nonlinear  algebraic  functions  in  the  output  equations  for 
p  quantities,  y  of  interest.  Much  is  known  about  the  steady  state 
engine  response.  This  can  be  written: 

0  “  f(xSS’  u*  0) 

(3.97) 

y  “  h(x$s>  a,  0) 

(3.98) 
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Most  of  the  modeling  effort  should  go  into  matching  steady  state 
response  since  it  is  at  these  points  that  engine  performance  is 
usually  evaluated.  Dynamical  response  is  also  important  to  the 
overall  description  of  the  behavior  since  these  dynamical  proper¬ 
ties  determine  stability. 

A  linear  approximation  to  equation  (3.95)  can  be  evaluated  if 
derivative  matrices  of  f  are  available  as  described  in  Section  3.6. 
The  function,  f,  in  Eq.  (3.95)  can  be  expanded  about  an  equilibrium 
point  (x0,uQ)  as  follows: 

6x  *  fx(VU0)4x  +  fu<xo'uo){u 

Sx  =  x-x  (3.99) 

o 

6u  *  U-UQ 

If  it  is  assumed  that  u(t)  is  piecewise  constant  on  an  interval 
[nTj(n+l)T],  equation  (3.95)  can  be  rewritten  relative  to  the  con¬ 
stant  value  of  u  and  the  value  of  x  which  would  be  reached  if  u 
remained  at  that  value,  or, 

x  «  Fx(*ss»  u(n))  [x(t)-xss]  (3.100) 

where  xcc  must  satisfy  the  nonlinear  equilibrium  relationship 

0  -  f(xss,  u(n),  8)  (3.101) 

or,  equivalently, 

xss  “  s(u(n)»  Q)  (3.102) 

where  g(u(n),  8)  is  defined  as  the  reference  schedule  of  the  states 
given  the  control  input  levels. 

The  time  interval,  T,  can  be  made  small  relative  to  the  dy¬ 
namics  so  that  the  constant  control  assumption  is  valid  for  any 
u(t). 
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Equation  (3.100)  is  the  linear  simulation  of  the  nonlinear 
dynamics  in  Eq.  (3. 95).  It  is  an  accurate  model  only  when  6x 
and  6u  remain  within  a  suitable  region  of  linearity.  Unfortunately, 
the  controls  in  complex  engines  often  change  quickly  and  the  engine 
operates  away  from  the  equilibrium  points  during  most  large  power 
level  changes.  The  concept  of  the  linear  simulation  of  the  form 
in  Eq.  (3.99)  can  be  expanded  to  produce  a  more  accurate  represen¬ 
tation  of  the  true  nonlinear  dynamics. 

The  model  reduction  procedure  described  in  Section  3.6  pro¬ 
duces  low  order  linear  models  parameterized  with  terms  which  are 
important  to  the  response.  These  models  represent  the  important 
dynamic  behavior  of  the  engine  throughout  its  flight  envelope.  The 
first  step  in  the  nonlinear  model  development  procedure  produces 
polynomial  functions  which  fit  the  dynamic  parameters  with  engine 
operating  variables.  This  is  a  curve- fitting  problem  which  can 
be  stated  as  follows.  For  a  set  of  operating  variables,  z,  e.g., 
rotor  speeds,  ambient  conditions,  and  control  inputs,  terms  of 
the  following  form  are  created  at  each  model  generation  point: 

t.(z)  *  n  Z.  (3.103) 

J  i*l 

The  set  of  coefficients,  ,  which  minimize  the  following 
cost  function  are  found: 

J*  •  min  !£  fp(k)  -  £  a j  t(z(k))l  J  (3.104) 

a^.r.q  (fo  L  j-1  3  /J  l 

where  the  non-zero  elements  of  the  reduced  dynamics  matrix,  p (k ) 

at  each  of  N  flight  points  are  fit  with  a  set  of  polynomial  func¬ 
tions  of  the  operating  variables  t(z(k)j.  (Note  that  the  set  of 
variables,  z,  may  contain  other  p's).  The  minimization  is  carried 
over  the  constant  coefficients,  a^,  the  number  of  terms  used,  r, 
and  the  highest  powers  represented  in  each  term,  q.  This 


procedure  is  known  as  optimal  subset  selection  in  a  ridge  regres¬ 
sion  problem  [18] .  In  Section  5.2,  this  procedure  is  illustrated 
for  the  GE  23. 

The  resulting  dynamic  matrix  can  be  written  as  follows: 

|f  *  F(x)  (3.105) 

x=0 

i.e.,  the  matrix  F  represents  the  gradient  of  the  nonlinear  func¬ 
tion,  f (x,u)  along  the  steady  state  operating  line,  A  Taylor 
series  expansion  of  f(x,u)  matches  the  function  at  an  arbitrary 
(x,u)  to  the  gradient  at  the  static  operating  line,  F(x).  This 
development  is  shown  in  Appendix  C.  The  results  indicate  that  to 
match  the  nonlinear  dynamics  at  arbitrary  point  to  order  ||6x||5, 
the  following  form  can  be  used: 


where 


f(x,u)  •  ifFfXj)  ♦  F(X2)]  *  o((x-x0)S) 

(3.106a) 

£  F  (x,xss) 

(3.106b) 

x-x 

X.  ■  X  ♦  - - 

1  0  <*1 

(3.107) 

x*x 

X,  •  xrt  ♦  - - 

2  0  12 

(3.108) 

-  3*/!T 

(3.109) 

^2  *  3-/31 

(3.110) 

*0  *  g(u) 

(3,111) 
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The  implication  of  Eqs.  (3.106)  through  (3.111)  is  that  the  linear 
simulation  in  Eq.  (3.100)  can  be  replaced  by  a  far  more  accurate 
( j  1 6x | {  vs.  j  1 5 x J | 5)  nonlinear  simulation  as  follows: 


x  *  F(x,x  ) (x-x  ) 
ss  ss 

xss  =  g(u(n),  8) 


(3.112) 


This  method  can  be  illustrated  with  the  following  simple 
example.  Suppose  the  nonlinear  model  is  the  following  scalar  equa¬ 
tion  (which  is  assumed  unknown) : 


k  »  X4  -  2x3  u2  *  u3  (3. 113) 

Gradients  of  this  function  were  available  and  the  linearized  models 
were  fit  to  form  fitted  dynamics  matrix: 

F(x,u)  =  (4x3  -  6x2  u2)  (3.114) 


The  reference  point  schedule  indicates  that  usg  *  1,  xRS  *  1,  is 
a  static  equilibrium.  The  “engine  simulation"  for  a  step  input  to 


2 '  3+/S1 


(3,115) 


The  model  is  (3.94)  can  be  compared  to  the  actual  equation  in  (3.92) 
as  follows: 


•  1 1 
x  • 

2 1 


x-l  3  x-1  31 

(1  *  w?  * “  *  j 


«f»  *  *  in 


]»)  ]  |(x-i) 


(3.116) 


-  |x5  ♦  x2  ♦  x  ♦  1  -  2(x2*x*l) j  (X‘l) 


*  x4  -  1  -  2(x^-l) 


*  x4  -  2x5  ♦  1 
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i.e.,  the  procedure  creating  Eq.  (3.115)  from  the  gradients  pro¬ 
duces  the  exact  representation  of  Eq.  (3.113)  since  the  original 
dynamics  function  is  a  polynomial  of  order  less  than  x^. 

3.8  SUMMARY 

Several  procedures  for  developing  important  components  of 
the  multivariable  controller  have  been  presented.  Trajectory 
generation  methods  have  been  discussed  which  produce  feedforward 
commands  to  the  plant.  Output  regulator  design  methods  have  been 
derived  which  result  in  simplified  regulator  gains  for  the  pertur¬ 
bation  logic.  Finally,  a  model  analysis  methodology  has  been 
explained  which  uses  a  complex  nonlinear  engine  simulation  to  pro¬ 
duce  low  order,  linear  design  models,  reduced- parameter  linearisa¬ 
tions  and,  finally,  simplified  nonlinear  equations  for  the  plant. 

In  Chapter  IV,  the  structure  of  a  multivariable  controller  for 
the  GE  23  is  presented  which  addresses  the  control  requirements 
discussed  in  Chapter  III  and  uses  procedures  presented  in  Chapter  III 
to  produce  a  reliable  and  efficient  full  envelope  control  law  for 
the  GE  23  variable  cycle  engine. 
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SECTION  IV 

CONTROL  STRUCTURE  DEFINITION 


4.1  INTRODUCTION 

The  design  requirements  for  a  controller  for  the  General 
Electric  GE23-JTDE  variable  cycle  engine  were  developed  in  Chapter 
II.  The  purpose  of  this  chapter  is  to  describe  a  control  structure 
and  design  methodology  proposed  by  SCI  (Vt)  which  satisfies  these 
requirements . 

The  proposed  controller  is  diagramed  in  Figure  4.1.  It  uti¬ 
lizes  a  model  following  structure  for  the  control  law,  and  incor¬ 
porates  a  nonlinear  engine  model  that  is  common  to  several  of  the 
functional  blocks.  Transient  regulation,  steady-state  performance 
and  failure  accommodation  criteria  can  be  met  with  this  technique. 
The  primary  purpose  of  this  controller  is  to  operate  the  engine 
at  desired  levels  of  thrust  and  airflow  in  the  presence  of  dis¬ 
turbances,  failures,  and  plant  variations  without  violating  con¬ 
straints  on  component  stability  and  physical  limits.  Any  flexi¬ 
bility  available  after  these  primary  goals  have  been  attained  is 
used  to  maximize  secondary  goodness  criteria  such  as  response 
time,  fuel  consumption,  and  life. 

The  hardware  environment  in  which  the  control  logic  will  be 
implemented  is  defined  in  terms  of  the  FADEC  microprocessor  con¬ 
figuration  (Figure  4.2).  The  signals  entering  and  leaving  the 
controller  are  defined  in  the  figure.  They  are  listed  and  cate¬ 
gorized  in  Table  4.1.  These  signals  can  be  grouped  into  gas  path 
sensor  outputs,  actuator  stroke  feedbacks,  logic  inputs,  operating 
point  inputs,  actuator  commands  and  logic  outputs.  A  precise 
definition  of  these  quantities  is  important  to  initial  structure 

r  definition. 
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IlGURE  4.1  CONTROL  LOG l C/COMPUTER 


Table  4.1 

Input/Output  Signal  Categories 


SIGNAL 

SOURCE 

CONVERSION 

TYPE 

ng 

F/0 

:  Sensor  Input 

NF 

F/0 

Sensor  Input 

T2 

Sensor 

A/D 

Sensor  Input 

t2.5 

Sensor 

A/0 

Sensor  Input 

tab 

Sensor 

A/D 

Sensor  Input 

P53 

Sensor 

F/D 

Sensor  Input 

*P21 

Sensor 

F/D 

Sensor  Input 

ps21 

Sensor 

F/0 

Sensor  Input 

*P23 

Sensor 

F/0 

Sensor  Input 

ps23 

Sensor 

F/D 

Sensor  Input 

LOO 

Sensor 

Switch 

Logic  Input 

3F  («) 

Sensor 

P/0 

Sensor  Input 

flc  (FB) 

Sensor 

P/O 

Sensor  Input 

WFN  (FB) 

Sensor 

P/0 

Sensor  Input 

FBTV  (FB) 

Sensor 

P/O 

Sensor  Input 

VABIF  (FB) 

Sensor 

P/0 

Sensor  Input 

VASIA  (FB) 

Sensor 

P/O 

Sensor  Input 

ASA  (FB) 

Sensor 

P/0 

Sensor  input 

«FR  («) 

Sensor 

p/0 

Sensor  Input 

AS  (FB) 

Sensor 

P/O 

Sensor  Input 

W.A 

Date  Bus 

Thrust  Request 

HOBS 

Oita  Bui 

Mode  Select 

POWER 

Oata  Bus 

On/Off 

INLET 

Data  Bus 

FAIL 

Oata  Bus 

Logic  Output 

STATUS 

Oata  Bus 

Binary  Output 

H 

Modulator 

PMM 

Actuator  Comand 

«c 

Modulator 

PMH 

Actuator  Comand 

WFN 

Modulator 

PMN 

Actuator  CoHwand 

FBTV 

Modulator 

PUN 

Actuator  Cotaund 

VABIF 

Modulator 

rWT 

Actuator  CcMwnd 

VASIA 

Modulator 

PUN 

Actuator  Comand 

ASA 

Modulator 

rKW 

Actuator  Coeaund 

WFR 

Modulator 

PtM 

Actuator  Coaaaand 

AB 

Modulator 

PUN 

Actuator  Comand 

MAIN  IGNITION 

Driver 

Suite* 

Output  Logic 

AUGNENTOR  tCNITION 

Driver 

Suite* 

Output  logic 

A  functional  description  of  the  controller  is  presented 
in  Table  4.2.  Described  in  the  sections  below  are  the  functional 
blocks  of  logic  which  accomplish  these  indicated  tasks  as  well  as 
the  signal  flow  paths  between  the  blocks. 

Note  that  the  control  logic  is  discussed  from  a  continuous 
time  viewpoint  since  this  approach  tends  to  be  more  clear.  Actual 
implementation  of  the  logic,  however,  will  be  through  full  dis¬ 
cretization  of  the  equations.  Program  structure  will  be  a  multirate/ 
multiloop  design.  Depending  on  the  availability  of  loop  clocks,  the 
controller  may  operate  in  a  free -running  mode  where  actual  loop 
times  are  variable.  Synthesis  techniques  are  available  to  accommo¬ 
date  this  mode  of  operation. 


4.2  OVERVIEW  OF  SYSTEM 

The  operation  of  the  interconnected  system  is  described  below* 
Figure  4.3  shows  a  detailed  diagram  of  the  signal  connections 
defined  at  this  time.  Details  of  each  block  are  presented  in  the 
following  section. 

Pilot  and  inlet  commands  are  the  primary  performance  setting 
variables.  The  pilot  input  is  treated  as  a  thrust  request  and  an 
inlet  input  as  an  airflow  request.  In  addition  to  these  two  con¬ 
tinuous  inputs,  a  digital  mode  input  is  provided  to  define  the 
engine  performance  and  operating  line  behavior  at  various  mission- 
oriented  flight  conditions.  The  modes  are  discussed  in  the  des¬ 
cription  of  the  reference  point  schedules.  These  input  quantities 
are  used  to  set  the  desired  operating  point.  Ambient  pressure  and 
temperature  signals  are  synthesized  in  the  fault  tolerant  filter 
for  the  remaining  degrees  of  freedom  in  specifying  the  desired 
performance. 

Reference  point  schedules  produce  an  estimate  of  the  equili¬ 
brium  point  quantities  at  the  requested  operating  points.  These 
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schedules  represent  the  static  portion  of  the  nonlinear  engine 
model. 

The  estimated  equilibrium  values  generated  by  the  reference 
point  schedules  are  passed  to  the  trajectory  generator.  This  block 
uses  the  dynamic  engine  model  to  form  compatible,  nominal  trajec¬ 
tories  in  the  engine  variables  and  controls  to  move  from  the  present 
point  to  the  desired  point  in  an  optimal  fashion  without  limit 
exceedances. 

The  nominal  path  specifications  for  the  regulator  are  the 
trajectory  generator  outputs.  The  regulator  uses  sensed  engine 
state  information  compared  to  these  nominal  values  to  form  an 
actuator  command.  The  regulator  has  no  time  varying  components. 

The  d.c.  gains  from  set  point  quantities  to  controls  will  be 
designed  so  that  acceptable  static  performance  can  be  achieved 
without  the  integral  trim  or  d.c.  gain  enhancing  compensation. 

Outputs  from  the  regulator  are  actuator  commands.  These  are 
processed  through  output  failure  detection  and  accommodation  logic 
to  prevent  abnormal  actuator  inputs  to  the  engine.  Limit  protection 
overrides  appear  in  this  block.  Also,  sensed  gas  path  a'td  aetuator 
stroke  inputs  are  processed  through  a  fault-tolerant  filter  which  is 
designed  to  produce  the  best  estimate  of  the  required  regulator 
inputs,  given  a  set  of  noisy  and  failed  sensors.  Passive  adaptation 
of  logic  to  ambient  and  power  conditions  is  accomplished  with  a 
centralized  gain  schedule. 

Failure  accommodation  is  provided  throughout  the  logic.  Thrust, 
airflow,  and  mode  inputs  are  monitored  according  to  slew  rate  and 
allowable  range*,  sensor  inputs  have  full  analytical  redundancy  in 
the  fault-tolerant  filter;  control  commands  are  monitored  for 
abnormal  conditions;  and,  servo  error  signal  response  is  monitored 
to  detect  failed  output  servos.  The  procedures  proposed  are  discussed 
in  a  later  section. 

Section  4,3  describes  the  engine  model  setup.  The  function 
of  the  reference  point  schedules  is  developed  in  Section  4.4, 


Several  approaches  to  the  trajectory  generator  are  described  in 
Section  4.5.  Section  4.6  describes  the  format  of  the  regulator 
and  the  procedure  for  deriving  feedback  gains.  Fault  tolerant 
filtering  and  the  engine  protection  blocks  are  described  in  Sec¬ 
tions  4.7  and  4.8.  Finally,  failure  detection  and  accommoda¬ 
tion  protocols  are  discussed  in  Section  4.9. 

4.3  NONLINEAR  MODEL 

A  simple,  accurate  description  of  the  dynamic  and  static 
behavior  of  the  engine  is  required.  The  model  structure  should 
allow  separate  programming  of  steady  state  and  dynamic  response. 

The  implementation  of  the  model  will  be  in  a  subroutine  format  so 
that  various  control  blocks  can  utilize  the  same  code  with  slightly 
different  inputs  and  outputs. 

The  initial  approach  is  described  in  Chapter  III  and  is  sum¬ 
marized  below.  The  engine  operating  point  is  derived  as  a  func¬ 
tion  of  the  inputs.  Thus,  PLA,  mode,  TT2  and  PT2  (derived)  repre¬ 
sent  a  complete  set.  Set  point  schedules  currently  used  will  be 
incorporated  so  that  engine  set  point  performance  will  meet 
specification.  Peripheral  variables  such  as  sensor  output  and 
non-scheduled  actuator  deflections  are  also  predicted  in  this 
portion  of  the  model.  Nonmeasured  outputs  such  as  thrust,  stability 
margin  and  airflows  can  also  be  approximated.  The  general  form  is: 

y  •  h(x,u)  (4.1) 

where  u  is  a  group  of  terms  that  specify  the  set  point  uniquely 
and  v  is  the  estimate  of  all  other  states,  outputs,  measurements, 
and  controls  of  interest  at  this  set  point  The  functions,  h(x,u), 
will  consist  of  tables  or  polynomial  fits  depending  on  the  most 
efficient  implementation. 

The  dynamics  of  the  model  will  be  constructed  from  reduced 
order  linear  models  produced  at  various  points  in  the  flight 
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envelope  and  mapped  by  a  group  of  low  order  functions  of  auxiliary 
variables.  These  functions  will  most  likely  be  implemented  as 
polynomials  of  low  order.  Sensitivity  methods  will  be  used  to 
reduce  the  problem  to  a  tractable  for:s.  The  model  for  the  linear 
dynamics  is: 

F(x,u)  «  ”(x,u)  (4.2) 

oX 

where  f(x,u)  are  the  nonlinear  plant  equations. 

x  =*•  f(x,u)  (4.3) 

The  advantage  of  using  the  gradients  rather  than  f  directly 
lies  in  the  independent  calculation  of  the  reference  point.  The 
dynamic  model  has  the  following  form: 

x  *  #(F(x(u)}  ( ,x*x1$)  (4.4) 

where  4^(.)  is  a  linear  functional  described  in  Appendix  C, 

4.4  REFERENCE  POINT  SCHEDULES 

The  reference  point  schedules  are  basically  given  in  Equation 
(4.1)  with  x*xe  Figure  4.4  shows  the  reference  point  schedules 
in  more  detail.  A  procedure  for  generating  logic  is  discussed 
below. 

The  inputs  to  the  schedules  are  FLA*  mode,  inlet  commands  and 
derived  ambient  variables  TT2*  PT2.  PLA  requests  thrust  continuous¬ 
ly  and  freely  from  idle  to  maximum  power,  inlet  inputs  are  also 
available.  Ambient  variables  determine  the  gas  state  at  the  engine 
face.  The  mode  command  is  a  means  of  moving  the  operating  point 
to  a  region  near  the  nominal  point  in  order  to  satisfy  specific 
operating  constraints.  These  mode  criteria  are  shown  in  Table  4.3, 
where  the  aircraft  envelope  is  divided  as  shown  in  Figure  4.5. 
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Figure  4.4  Reference  Point  Schedule  Generator 


Table  4.3 

Mission/Aircraft  Definition 
(Control  Performance  Criteria) 


l 
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igure  4.5  Control  Criteria  Within  Flight  Envelope 


The  data  have  been  used  to  define  three  modes  of  operation 
which  will  be  developed  to  satisfy  a  group  of  flight-dependent 
constraints.  Each  control  criterion  in  Table  4.3  can  be  grouped 
within  one  of  four  constraints  shown  in  Table  4.4.  The  first 
constraint  is  satisfied  by  the  engine  manufacturer's  operating 
point  design  as  supplied  by  the  control  schedules.  Further  opti¬ 
mization  in  this  mode  could  be  achieved  by  outer  loop  adjustment 
of  the  operating  point  to  minimize  TSFC.  This  portion  of  the 
logic  will  be  optional  until  later  in  the  development  effort.  The 
second  constraint  is  based  upon  inlet  limits  and  matching.  The 
details  of  these  interactions  have  not  been  identified  at  this  point. 
The  third  constraint  involves  increasing  component  stability  margins 
to  accommodate  distortion,  weapons’  discharge  and  low  ambient  pressure 
effects.  The  increased  stability  will  be  produced  at  constant 
thrust  and  increased  fuel  consumption.  The  final  constraint  is  a 
high  response  match  point  which  will  provide  continuous  thrust 
modulation  at  increased  response  rates. 


Table  4.4 

Operating  Constraints 


Constraint 

Definition 

ci 

Within  engine  limits,  run  to  maximum  aug¬ 
mented  and  dry  thrust  points  with  mini¬ 
mum  TSFC  modulation  continuously  between 
them. 

C2 

Inlet  airflow  compatibility 

C3 

Increase  stability  margins  at  constant 

thrust. 

C4 

Maximize  thrust  response  rates  w/o 

*T 

discontinuities. 

Three  modes  are  defined  and  regions  of  the  envelope  where 
these  modes  are  operative  are  shown  (Table  4.5).  The  normal  operat 
ing  mode  is  shown  to  be  active  in  all  parts  of  the  envelope  except 
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Table  4.5 

Modes  of  Operation 


MODE 

SYMBOL 

DESIGN 

CONSTRAINT 

ACTIVE 

REGION* 

Normal 

Increased  Stability 

High  Response 

Mo 

M2 

C2*  C1 
c3,  C2,  Cj 

C4*  C1 

I,  II.  IV 

I,  II,  III 

IV 

*REGI0N$:  I  Supersonic  III  Left-Hand  Corner 

II  Cruise  IV  Approach 


Region  III  (left-hand  corner).  In  this  region,  selection  of  mode  0 
or  1  or  2  will  have  the  same  result.  Similarly,  M3  is  not  distin¬ 
guishable  from  mode  anywhere  except  approach. 

The  initial  calculation  of  Mq  (baseline)  reference  schedules 
will  be  accomplished  as  follows.  Steady  state  control  and  outputs 
will  be  derived  from  the  linear  model  base  points  and  the  nonlinear 
simulation  when  available.  The  engine  set  point  quantities  will 
be  determined  from  the  control  schedules  and  output  schedules 
provided  by  the  manufacturer.  It  is  anticipated  that  the  output 
schedules  will  have  to  be  modified  to  incorporate  control  satura¬ 
tions  which  are  only  implicitly  handled  in  the  schedules.  This 
effect  is  shown  in  Figure  4.6  for  the  Np  schedule.  At  low  PLA  set¬ 
tings,  the  variable  low  turbine  area  saturates  and  the  rotor  speed 
falls  off.  For  the  reference  point  logic,  a  more  feasible  schedule 
is  shown  in  the  figure.  Steady  state  data  is  schematically  repre¬ 
sented.  Predicted  values  of  Np  will  be  closer  to  actual  running 
values  and  control  action  will  not  tend  to  compensate  inappropriately 
for  large  offsets.  A  flow  chart  of  the  reference  point  calculation 
is  shown  in  Figure  4.7.  The  PLA  and  mode  inputs  form  the  initial 
set  point  command.  Mode  biasing  will  be  discussed  below.  The  first 
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Figure  4.6  Control  Saturation  Effects  on  N,,  Schedule 
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Figure  4.7  Reference  Point  Calculation  Flowchart 
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calculation  uses  the  set  points  as  inputs.  Control  actuator  satu¬ 
rations  on  output  are  reset  to  their  maximum  or  minimum  value  and 
the  routine  is  entered  with  these  new  inputs.  The  resulting 
reference  point  is  an  estimate  of  the  desired  equilibrium  within 
operating  limits  and  without  actuator  limit  saturation.  Set  point 
schedules  and  actuator  limits  will  be  constructed  with  a  small 
'’overlap"  so  that  operating  points  with  predicted  control  satura¬ 
tions  will  actually  cause  the  engine  to  operate  with  its  controls 
limited.  This  outer  bound  is  illustrated  in  Figure  4.6. 

Mode  switching  will  be  handled  by  alteration  of  the  reference 
point.  Transition  and  regulation  will  be  accommodated  as  in  any 
other  operating  point  shift.  This  procedure  simplifies  the  logic 
significantly. 

The  procedure  for  generating  control  mode  shifts  involves 
generating  a  group  of  trial  biases  from  the  linear  models.  These 
biases  will  be  evaluated  on  the  nonlinear  simulation  to  determine 
actual  performance  level  shifts.  Final  schedules  will  be  deter¬ 
mined  from  these  results.  In  general,  the  simplest  structure  will 
be  sought. 

This  procedure  is  described  for  the  increased  stability  mode. 
In  this  case,  the  fan  and  compressor  stabilities  will  be  increased 
at  constant  thrust  and  airflow.  A  direction  in  state  space  for 
this  effect  can  be  deduced  from  the  linear  models  in  steady  state. 
The  general  motion  of  the  operating  point  for  changes  in  control 
can  be  written  as  follows: 


6y  »  H6u 

where  H  is  steady  state  performance  sensitivity  matrix  derived  from 
the  models.  By  choosing  appropriate  values  of  6y  to  be  zero  (e.g., 
6Fn  and  6Wfan)  and  regarding  others  to  be  positive  (6SMF,  6SMC) , 


appropriate  linear  combinations  of  the  6u's  can  be  calculated. 
Reasonable  values  for  these  can  be  used  to  generate  a  perturbed 
operating  point.  Actual  performance  shifts  must  then  be  evaluated 
in  the  nonlinear  deck. 

4 . 5  TRAJECTORY  GENERATOR 

The  outputs  from  the  reference  point  schedules  represent  a 
group  of  engine  variables  and  controls  which  satisfy  equilibrium 
conditions  and  which  are  within  physical  and  operating  limits. 

These  quantities  are  based  on  the  values  of  PLA,  mode,  and  inlet 
conditions.  These  may  instantaneously  change  from  sample  to 
sample  because  of  pilot  inputs.  The  actual  reference  outputs 
could  be  discontinuous  over  a  large  range.  If  these  reference 
values  were  linked  directly  to  the  regulator,  moderate  engine 
transitions  would  saturate  the  actuators.  To  rate  limit  the  refer¬ 
ence  points  or  the  PLA  commands  would  seriously  degrade  small 
signal  response.  Also,  since  the  system  response  to  very  large 
inputs  is  nonlinear,  the  response  achieved  without  some  input  com¬ 
pensation  would  be  suboptimal  (and  most  likely,  catastrophic).  The 
transition  generator  is  designed  to  produce  an  ideal  reference 
between  the  present  engine  state  and  the  state  most  recently 
requested  by  the  reference  values. 

The  nominal  reference  trajectory  should  have  the  following 
attributes: 

(1)  it  should  be  nearly  compatible  with  actual  engine 
response,  i.e.,  the  reference  input  trajectory  should 
nearly  produce  the  reference  output  trajectories, 

(2)  it  should  nearly  track  all  engine  and  actuator  limits, 
and 

(3)  it  should  exhibit  optimized  response  for  both  large  and 
small  inputs. 

The  approach  chosen  for  this  function  is  to  utilize  the  dynamic 
model  of  the  engine  compensated  with  a  nonlinear  feedback  law  and 
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driven  by  the  reference  point  schedule  as  the  generator  of  this 
nominal  path. 

Several  issues  can  be  addressed  at  this  point.  A  common  ques¬ 
tion  that  arises  concerns  the  requirement  for  a  trajectory  generator. 
Since  the  model  is  compensated  by  a  feedback  law,  a  similar  control 
strategy  could  be  used  on  the  actual  engine  without  the  need  for  a 
complex  software  block.  As  the  control  strategy  is  simplified,  this 
argument  becomes  more  viable.  However,  the  internal  model  can 
exhibit  the  response  of  the  engine  without  actuator  and  sensor 
dynamics.  The  frequency  response  of  the  model  can  be  made  quite 
high  in  order  to  produce  responsive,  limit  accommodating  paths. 

These  paths  can  be  used  at  the  regulator  inputs.  Stiff  regulator 
compensation  can  be  designed  to  assure  accuracy  of  the  actual  engine 
in  the  presence  of  plant  variation  and  disturbances.  Small  and 
large  signal  responses  are  optimized  in  this  process.  In  a  more 
complex  strategy,  the  model  can  be  used  iteratively  to  calculate 
on-line  optimized  trajectories  relative  to  nonlinear  cost  functions 
including  minimum  time.  Thus,  advanced  control  strategies  can  be 
applied  to  the  system  with  the  model  without  a  large  impact  on  the 
control  structure.  Since  this  is  a  viable  option  in  the  design, 

SCI  has  chosen  to  maintain  an  accurate  dynamic  model  in  the  control 
and  use  the  model  reference  structure  as  the  design  basis. 

The  preliminary  trajectory  generation  logic  is  shown  in 
Figure  4.8.  A  nonlinear,  proportional  override  logic  is  used  as 
the  compensation.  The  critical  element  in  the  design  is  that  the 
model  compensation  gains  are  generated  directly  from  the  regulator 
gain  schedules.  This  procedure  uses  an  asymptotic  assumption  on 
the  optimal  response.  The  procedure  results  in  a  SO  percent  reduc¬ 
tion  of  the  gain  schedule  storage  requirement.  Operation  of  the 
trajectory  generator  and  the  initial  design  approach  is  discussed 
below. 

For  small  transients,  reference  inputs  are  used  as  direct  com¬ 
mands  to  the  engine  model.  The  mode  response  is  further  compensated 
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with  output  gains  which  are  proportional  to  the  optimal  regulator 
gains.  When  no  limits  are  exceeded  by  the  model,  the  response  of 
the  model  is  controlled  by  the  locally  linear  gains.  When  a  model 
output  approaches  a  limit,  the  servo  rate  limits  are  proportionally 
reduced.  This  has  the  effect  of  transferring  the  control  law  from 
the  unlimited  regulator  to  one  of  several  specifically  designed 
multivariable  limit  loops.  The  rate  limit  is  proportional  to  the 
error  signal  from  the  limiting  blocks.  This  feedback  in  one  or 
more  of  the  control  channels  will  tend  to  cause  the  model  to  move 
smoothly  onto  a  limit.  If  the  proportional  error  signal  indicates 
that  the  system  can  move  away  from  the  boundary,  the  limit  is 
smoothly  removed.  The  feedback  gains  for  the  limit  loops  are 
designed  using  output  weightings  on  the  specific  engine  constraints. 
This  yields  a  feedback  vector  for  the  limit  loop.  This  vector  will 
be  simplified  using  sensitivity  calculation  and  scheduled  as  a 
function  of  ambient  conditions. 

The  large  input  performance  of  the  system  will  tend  to  behave 
in  a  time  optimal  fashion.  This  assertion  is  justified  from  optimal 
solution  to  the  minimum  time  problem  for  a  linear  system.  Here, 
optimal  trajectories  consist  of  a  minimum  time  (corresponding  to  a 
bang-bang  control)  trajectory  to  a  limit,  tracking  the  limit,  and 
moving  off  the  limit  to  the  final  point.  The  character  of  the  tra¬ 
jectories  generated  in  this  method  will  be  very  similar  to  this 
type  of  motion  without  the  requirement  of  explicit  solution  of  non¬ 
linear  optimization  problems. 

A  second,  higher  risk  option,  will  be  investigated  during  the 
second  phase  of  the  program.  Using  the  system  model  and  auxiliary 
differential  equations,  solutions  to  optimal  trajectories  can  be 
calculated  in  real  time.  Trajectories  of  this  type  will  be  generated 
off-line  and  used  to  tune  the  trajectory  generation  implemented 
above.  Specific  performance  tradeoffs  will  be  investigated.  If 
warranted,  further  development  of  on-line  optimization  will  then  be 
incorporated  into  the  transition  logic. 
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It  is  also  possible  to  modify  the  trajectory  generator  logic 
to  include  a  variable  transient  response  mode  for  increasing  engine 
life. 

Two  possible  strategies  are  proposed.  Both  require  specifica¬ 
tion  of  a  "transient  life  function"  (TLF)  in  the  form 

fT 

TLF  =  /  f(T4,  TB4,  N2S...)dt 

Technique  1: 

The  current  philosophy  behind  the  trajectory  generator  is  to 
produce  a  feedforward  control  which  takes  an  engine  to,  and  holds 
it  on,  its  operating  limits  (T4,  fB4,  N25»  etc.).  These  limit  sched¬ 
ules  are  the  same  as  those  employed  in  the  engine  protec¬ 
tion  block  of  the  controller  logic.  By  modifying  the  limit  sched¬ 
ules,  the  engine  could  be  made  to  track  a  less  "severe"  limit.  The 

transient  response  would  be  degraded,  but  the  engine  life  would  be 

* 

increased.  A  simple  example  would  bo  to  require  a  lower  TB4  at 
higher  values  of  1‘4  as  indicated  in  Figure  4.9.  Specification  of 
the  actual  schedules  requires  knowledge  of  the  TLF  function. 

Technique  2: 

A  full  nonlinear,  optimal  life,  trajectory  generation  procedure 
is  possible.  This  is  a  higher  risk  approach  than  that  proposed  as 
Technique  1,  but  it  may  provide  a  corresponding  payoff.  The  calcu¬ 
lation  could  be  performed  on-line  replacing  the  present  Trajectory 
Generator,  or  it  could  be  performed  off-line  and  the  results  (trends, 
etc.)  incorporated  into  the  current  trajectory  generator.  These 
calculations  would  be  impossible  without  the  TLF. 

4.6  MULTI  VARIABLE;  REGULATOR  LOGIC 

The  outputs  of  the  trajectory  generator  are  continuous  control 
commands  and  engine  variable  references  which  predict  the  transition 
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Figure  4.9a  T,  Limit  Schedule 


H 


* 

Figure  4.9b  TB^  Limit  Schedule  for 
Increased  Engine  Life 


of  the  engine  from  one  point  to  another  without  exceeding  physical 
or  operational  limits.  If  this  prediction  were  exact  or,  alternately, 
if  the  mathematical  model  matched  the  engine  exactly,  there  would 
be  no  reason  to  include  the  regulator.  The  regulator  is  designed  to 
cause  the  engine  to  track  small  perturbations  from  a  nominal 
trajectory.  Further,  it  attempts  to  cause  the  engine  to  track 
the  reference  in  steady  state.  As  the  model  accuracy  gets 
better,  a  higher  bandwidth  can  be  sustained  by  the  regulator  logic 
and  the  d.c.  accuracy  will  be  proportionally  improved.  The  regulator 
model  is  proportional-only  control  with  the  set  point  accuracy  and 
stability  characteristics  designed  using  linear  optimal  theory. 


Two  procedures  will  be  used  to  calculate  the  optimal  regulator 
gains.  The  first  procedure  will  use  a  reduced  order  state  model 
produced  from  the  linearized  engine  model.  Linear  optimal  control 
laws  will  be  designed  using  output  weighting.  The  state  feedback 
laws  will  be  converted  to  output  feedback  laws  using  the  invorti- 
bility  of  the  output  distribution  matrix  partitioned  as  follows: 


Model : 

Control  Law: 

Output  Partition: 

Output  Control  Law: 


x  “  Fx  ♦  Gu 
u  ■  Cx 
y  *  Hx  ♦  Du 

u  •  (I+CH^DrWV 


An  alternate  procedure  is  preferred  over  this  method.  The 
models  are  reduced  to  include  the  important  engine  dynamics  and 
then  augmented  by  important  sensor  and  actuator  time  constants. 

The  output  feedback  matrix  is  then  directly  calculated  by  minimiz¬ 
ing  quadratic  cost  function  weighting  the  response.  This  procedure 
is  more  complex  than  the  LQG  approach;  however,  it  offers  several 
advantages.  In  particular,  sensitivity  calculations  can  be  used  to 
reduce  the  number  of  feedback  gains  and  the  system  can  be  reopti¬ 
mized  for  the  fixed  structure  response.  Thus,  a  single  synthesis 
algorithm  will  yield  the  optimized  control  law.  The  procedure 
offers  a  large  flexibility  in  model  and  feedback  components  which 
is  more  difficult  to  realize  in  the  LQG  approach. 
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The  accommodation  of  variation  in  flight  conditions  and  the 
development  of  the  gain  schedule  algorithm  will  be  accomplished 
using  several  procedures.  In  addition  to  the  ad  hoc  methods  of 
gain  fitting,  several  systematic  approaches  will  be  explored.  An 
attractive  approach  involves  using  the  output  regulator  design 
routine  to  optimize  a  fixed  structure  control  law  for  a  group  of 
operating  points.  Additional  operating  points  are  incorporated 
into  the  design  group  until  a  performance  degradation  threshold  is 
exceeded.  Gains  are  then  held  at  the  optimal  value  within  the 
design  envelope  and  scheduled  in  regions  between  the  envelopes.  A 
variation  in  this  procedure  is  to  assume  a  functional  relationship 
in  the  gain  law  which  is  derived  from  point  designs.  The  param¬ 
eters  of  gain  relationship  are  then  optimized  to  produce  a  point- 
wise  optimal /nonlinear  control  law  or  schedule.  The  procedures 
will  be  investigated* 

The  gain  schedule  algorithm  will  be  used  by  the  regulator  and 
the  transition  logic  for  gain  inputs.  Transition  generator  will 
use  a  proportional  gain  constant  to  vary  the  gain  values  to  achieve 
the  desired  model  response. 

4.7  FAULT  TOLERANT  FILTER 

The  signals  entering  the  controller  represent  various  types 
of  discretized  information.  These  signals  provide  many  orders  of 
redundancy  concerning  the  actual  measured  quantities.  The  filter 
block  functions  to  provide  the  best  available  estimate  of  required 
information  to  the  control  law.  The  fauit  tolerant  filter  must 
operate  on  all  analog  sensor  channels  except  PLA  and  inlet  commands 
to  provide: 

(1)  noise  attenuation, 

(2)  dynamic  compensation, 

(3)  error  correction,  and 

(4)  fault  insensitivity. 
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The  filter  uses  the  dynamic  model  of  the  engine  as  the  primary 
means  for  providing  analytical  redundancy  necessary. 

The  most  important  function  of  the  filter  is  to  provide  atten¬ 
uation  of  errors  due  to  various  sensor  noise  sources  and  compensate 
sensor  dynamics  which  would  otherwise  compromise  control  performance. 
Accurate  d.c.  response  must  not  be  degraded  by  the  filter.  The  form 
chosen  will  be  a  series  of  decoupled  extended  or  nonlinear  Kalman 
filters  with  reduced  gain  matrices.  The  gains  may  be  a  function  of 
the  flight  condition.  The  filter  update  will  be  provided  by  the 
nonlinear  dynamical  engine  model.  The  design  goals  of  this  block 
will  be  high  frequency  noise  rejection  and  d.c.  accuracy.  Foldover 
and  aliasing  due  to  sampling  high  frequency  noise  will  be  addressed 
by  the  design.  Analog/digital  prefiltering  requirements,  random 
sampling  time  algorithms  and  adaptive  noi.e  rejection  techniques 
will  be  investigated. 

A  sensor  failure  detection  algorithm  will  be  associated  with 
the  filter.  The  purpose  of  this  block  will  be  to  correlate  filter 
residuals,  sensed  levels  and  model  outputs  into  inferences  on  soft 
and  hard  channel  failures.  Multiple  failures  will  be  detected. 

This  block  will  examine  all  gas  path  variables  (including  TT2) . 

In  addition  to  the  detection  function,  sensor  failure  accommodation 
will  be  provided.  This  will  consist  of  manipulation  of  filter  in¬ 
puts  to  reduce  the  effects  of  bad  sensor  channels  and  optimize  the 
integrity  of  the  filter  outputs.  Also,  logic  indications  will  be 
passed  to  the  central  failure  accommodation  block  discussed  below 
for  GO/NO-GO  and  secondary  diagnostic  processing. 

4.8  ENGINE  LIMIT  LOGIC 

The  outputs  from  the  trajectory  generator  are  measured  against 
their  limits  and  are  used  to  produce  nominal  trajectories  which 
track  predicted  engine  response.  Because  of  model  inaccuracies, 
sensor  and  actuator  lags  and  nonlinearities,  build  differences, 
etc.,  it  is  possible  that  actual  engine  limits  may  be  exceeded 
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transiently  or  in  steady  state.  It  is  currently  planned  that  no 
integral  control  will  be  used  to  compensate  these  types  of  errors. 

The  engine  limit  logic  provides  additional  d.c.  gain  to  the  system 
when  limits  are  approached  or  slightly  exceeded.  Also  incorporated 
with  this  logic  is  a  failure  indication  concerning  the  overall 
control  operation. 

The  purpose  of  this  logic  is  to  produce  control  action  which 
tends  to  keep  the  engine  state  under  its  limits.  As  thresholds  are 
exceeded,  this  control  action  will  have  increasing  authority.  This 
initial  design  philosophy  is  sketched  in  Figure  4.10. 

As  a  limit  is  approached,  a  nonlinear  factor  is  provided  to 
the  output  regulator  gains  specifically  designed  for  this  3 iir.it 
(These  gains  are  also  used  by  the  trajectory  generator.)  As  the 
exceedance  increases,  larger  control  effort  is  commanded.  Above 
a  certain  value,  it  is  concluded  that  the  control  logic  is  unable 
to  produce  safe  regulation.  At  this  point,  an  integrated  over-limit- 
flag  is  passed  to  the  error  accommodation  logic.  The  failed  condi¬ 
tion  represents  a  non-specific  failure  in  the  overall  control  func¬ 
tion  due  to  input  failure,  pilot  error,  actuator  failure,  engine 
component  deterioration,  engine  auxiliary  failure,  or  related 
causes.  The  control  regards  this  as  a  NO-GO  situation  and  appro¬ 
priate  backup  sequences  are  initiated 

A  preliminary  diagram  of  the  detailed  logic  is  presented  in 
Figure  4.11. 

4.9  FAILURE  ACCOMMODATION  LOGIC 

Failure  indications  are  monitored  in  a  central  location.  Logic 
is  provided  to  produce  coordinated  transfer  to  backup  options.  This 
logic  Integrates  with  the  hardware  BIT  capabilities  but  provides  a 
far  more  powerful  tool  in  assessing  the  function  of  the  overall  con¬ 
trol  in  running  the  engine.  Emphasis  is  placed  on  fault  detection. 
The  important  decision  point  is  the  transfer  to  the  backup  mode. 
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gure 


FIOWE  «.))  OVERVIEW  OR  ENGINE  LIMIT  LOGIC 


This  decision  is  not  made  when  a  failure  is  detected,  but  rather 
when  that  failure  has  caused  the  engine  to  operate  away  from  pre¬ 
dicted  nominal  performance  levels.  In  this  situation,  the  controller 
has  failed  since  it  no  longer  can  be  used  to  satisfy  its  primary 
requirement. 

Peripheral  failure  diagnostics  are  also  generated  for  cockpit 
and  flight  line  utilization.  Failure  accommodation  is  not  stressed 
except  in  removing  the  effect  of  noncritical  instrumentation  and 
bringing  the  engine  to  a  fail  safe  backup  transition  point  after 
a  critical  failure  has  been  detected. 

Inputs  are  monitored  by  the  logic.  Redundant  PLA  channels  may 
be  monitored.  Failed  channel  selection  will  be  based  on  fail-safe 
operation.  In  addition,  maximum  slew  rate  limits  will  be  detected. 

The  action  for  PLA  failures  will  be  a  NO-GO  state  and  transition  to 
the  backup  control  mode. 

Actuators  are  monitored  for  correct  servo  following.  This 
system  monitors  error  signal  magnitudes  compared  to  a  threshold 
which  is  a  function  of  the  command  level  and  changes  in  command 
levels.  Threshold  exceedance  is  considered  as  an  actuator  loop 
failure.  The  failure  could  be  either  the  actuator  or  the  LVPT, 
however.  If  the  failure  is  in  the  LVPT,  it  is  possible  that  the 
actuator  could  still  be  used  to  control  the  engine.  Logic  for 
accomplishing  this  is  proposed  in  Figure  4.12  (detection  of  the 
failure  is  addressed  below). 

Normal  operation  is  described  in  Figure  4.12a.  Here  the 
actuator  is  driven  by  an  error  signal  (compensated)  which  is  formed 
by  comparing  the  actuator  command  to  the  LVPT  output.  In  the  event  of 
a  failed  LVPT,  the  actuator  can  be  driven  by  an  estimate  of  the 
error  signal  as  indicated  in  Figure  4.12b.  Here  the  error  signal 
estimate  is  provided  by  an  internal  mathematical  model  of  the 
actuator  in  the  form  of  a  subroutine.  It  is  further  suggested  that 
one  common  subroutine  could  be  used  to  model  all  the  servoloops  in 
a  computer  storage  efficient  structure. 
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Figure  4.12a  Normal  Operation  of  Actuator  Servo  Loop 
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Figure  4.12b  Operation  of  Actuator  with  Failed  LVPT 


Once  the  switch  to  this  backup  mode  has  been  made,  the  band¬ 
width  of  the  actuator  "loop"  will  be  reduced.  Modeling  errors 
combined  with  the  loss  of  the  feedback  will  result  in  hangoff 
and  drift  in  the  actuator's  output.  Outer  loops  in  the  regu¬ 
lator  logic  and  engine  protection  logic  will,  however,  act  to 
reduce  these  effects. 

Using  this  strategy,  actuator  effort  is  estimated  without 
relying  on  information  contained  in  the  engine  output  variables 
(slow  dynamic  response)  as  would  be  the  case  if  the  actuator  loop 
was  included  as  part  of  the  fault  tolerant  filter.  Also,  a  switch 
to  this  backup  mode  will  have  negligible  effect  on  any  other  compo¬ 
nent  of  the  control  logic. 

It  is  anticipated  that  most  actuators  will  be  "non-critical ." 
That  is,  the  controller  will  be  capable  of  running  the  engine  within 
its  constraints  with  one  of  these  actuators  failed.  Consequently, 
reconfiguration  of  the  control  law  will  be  unnecessary.  The  "dis¬ 
turbance"  produced  by  a  failed  actuator  will  be  compensated  by 
aetion  in  the  remaining  healthy  control  loops  -  and  an  allowable 
(but  degraded)  engine  responso  will  be  maintained.  A  failure  in  an 
actuator  servoloop  can  now  be  accommodated  by  proceeding  as  if  the 
feedback  sensor  (LVPT)  was  at  fault,  i.e.  if  corrective  action 
is  taken  for  a  failed  LVPT  and  the  failure  is  in  fact  in  the 
actuator,  no  accommodation  incompatibility  exists. 

Some  actuators  may  be  "critical,"  That  is,  loss  of  one  of 
these  will  result  in  the  control  structure  being  incapable  of 
running  the  engine  within  its  constraints.  Consequently,  a  failure 
will  require  a  switch  to  either  hydromechanical  backup  or  a  recon¬ 
figuration  of  the  control  structure,  tfhieh  actuators  are  critical 
is  dependent  on  the  logic  that  is  implemented.  Therefore*  identi¬ 
fication  of  critical  actuators  cannot  be  addressed  until  the  control 
law  is  actually  designed. 
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The  failure  diagnostic  logic  consists  of  a  series  of  decision 
points  to  determine  the  ability  of  the  control  to  safely  operate 
the  engine.  When  failure  indications  in  a  group  of  monitored  points 
occur,  transition  to  the  backup  mode  is  initiated.  This  procedure 
involves  a  pilot  indication  of  control  failure  and  simultaneous 
movement  of  the  engine  operating  point  to  the  nearest  "fail  safe 
region. "  No  action  is  taken  if  the  engine  is  operating  in  a  safe 
region.  The  fail  safe  region  will  be  defined  in  detail  later. 
However,  initially  an  operating  point  in  the  non- augmented  power 
and  away  from  temperature,  pressure  and  rotor  speed  limits  is 
desired.  This  transition  is  treated  as  an  overriding  mode  switch 
to  the  reference  point  logic  and  normal  control  regulation  to  this 
point  is  attempted. 

4.10  SUMMARY 

The  passively  adaptive  control  structure  is  defined  for  the 
VC11,  Each  block  will  be  further  refined.  The  results  of  the 
Phase  I  activity  will  be  a  more  detailed  definition  of  each 
block  and  a  group  of  proposed  design  specifications  and  pro¬ 
cedures  for  one  or  more  synthesis  method  which  will  be  under¬ 
taken  in  Phase  II. 
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SECTION  V 

PRELIMINARY  DESIGN  RESULTS 


This  chapter  is  a  compendium  of  design  results  generated  dur¬ 
ing  the  Phase  I  study.  These  results  serve  two  purposes.  First, 
they  form  a  base  from  which  to  build  an  innovative,  practical 
controller  in  Phase  II.  Second,  they  both  demonstrate  and  verify 
the  design  procedures  presented  in  Chapters  III  and  IV.  This  is 
important,  since  many  of  the  proposed  procedures  have  never  been 
implemented  in  digital  propulsion  controls  and  represent  extensions 
to  the  present  levels  of  complexity  in  control  law  development. 

Section  5.1  summarises  the  results  of  the  linear  model  analy¬ 
sis.  A  set  of  reduced  order  and  reduced  parameter  models  is 
generated.  In  Section  S.2,  an  accurate  nonlinear  model  is  developed 
using  these  reduced  models.  This  nonlinear  model  is  computer 
storage  efficient  and  will  become  the  base  of  the  implemented  con¬ 
troller.  Section  5.3  discusses  the  design  and  simulation  of  com¬ 
pensation  for  the  actuator  servo  valves.  Finally,  Section  5.4 
presents  preliminary  simulation  results  of  the  trajectory  genera¬ 
tion  logic  of  Section  4.1.5. 

S.l  LINEAR  MODEL  ANALYSIS 

A  mathematical  model  which  characterises  accurately  both 
static  and  dynamic  operation  of  the  6E23  throughout  its  flight 
envelope  is  highly  nonlinear.  However,  there  exist  two  major 
advantages  to  linearising  Jiis  model  before  attempting  to  design 
a  controller.  The  first  is  that  control  logic  design  procedures 
for  linear  mathematical  models  are  well  developed  and  easy  to  use. 
The  second  is  that,  with  a  linear  model ,  important  dynamic  inter¬ 
actions  may  be  identified  and  analysed  using  well  known  time 
domain,  frequency  domain  and  eigensystem  analysis  procedures;  thus, 
physical  insight  may  be  gained  which  is  necessary  for  the  design 
of  a  practical  controller. 
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There  are  dangers  which  must  be  avoided  in  the  linearization 
process,  however.  First,  the  linear  models  may  be  inaccurate  due 
to  procedural  errors.  Choosing  too  large  a  linearization  interval 
and  ignoring  numerical  truncation  are  examples.  Second,  the  set 
of  linearization  points  chosen  may  ignore  important  couplings  and/ 
or  dynamics  which  arise  as  a  result  of  some  nonlinearity.  Conse¬ 
quently,  an  insuff icient  set  of  models  may  be  in  use.  Care  must 
therefore  be  exercised  to  validate  the  linear  models  before  and 
during  the  design  phase. 

Described  below  are  the  results  of  an  analysis  of  linear 
models  supplied  by  the  engine  manufacturer.  Section  $.1.1  des¬ 
cribes  how  the  linearization  points  were  specified;  Section  5.1, 2 
describes  the  form  of  the  models;  Section  5.1.3  justifies  the  re¬ 
duction  in  order  from  eleventh  to  fifth;  and  Section  5.1.4  demon¬ 
strates  the  use  of  second- stage  model  reduction  (see  Section  3.6.2) 
to  eliminate  terms  in  the  system  dynamics  matrix  (F  in  x»Px*Gu}. 

5.1.1  Model  Linearization  Points 

Model  linearization  points  must  be  chosen  to  explore  fully 
the  envelope  of  engine  operation  (Section  2.9).  Specifically, 
points  must  be  chosen  which  span  the  inlet  pressure,  density  and 
temperature  ranges  (e.g. ,  altitude  and  Mach  number)  at  ail  power 
levels.  Sixty  sueh  points  have  been  specified  during  Phases  I  * 

They  are  detailed  in  Figures  5,1  and  $*2  and  in  Tables  $, i  and  S.2. 
The  flight  points  are  displayed  against  altitude  and  Mach  number 
in  Figure  5. la  and  against  inlet  pressure  and  temperature  in 
Figure  5.1b.  The  power  levels  examined  at  these  flight  points  are 
pictured  in  Figure  S,2  and  summarized  in  Table  S, J.» 

Table  S. 2  presents  a  summary  of  the  sixty  combinations  chosen. 
Included  In  this  table  are. comments  explaining  why  each  flight 
point  was  selected. 


ALTITUDE 
MAX  ALTITUDE 


Figure  5.1a 
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Figure  S.lb 


Definition  of  Flight  Points  Vs.  Altitude  and 
Mach  Number 


INLET  TEMPERATURE 
MAX  "'INLET'  TWE’RATDW 

Definition  of  Flight  Points  Vs.  PT2  and  TT2 
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Table  5.1 

Definition  of  Power  Points 


-POINT 

DEFINITION 

A  and  2 

Intermediate  Power  (Maximum 

Dry  Power) 

B 

Bypass  Transition 
e  A14B  Open 

•  A16B  Closed 

$  FBTV  Scheduled 

C  and  4 

Power  Break 

•  A54  Saturation  Point 

•  Begins  to  Fall  Off 

D 

Off- Idle 

•  Point  at  20%  of  Distance 
between  Idle  and  C 

Measured  in  Thrust  Change 

1 

Maximum  Power 

2 

Above  Intermediate 

•  A16B  Open 

•  A14B  Closed 

•  FBTV  Closed 

5 

Idle 

5.1.2  Form  of  the  Linear  Models 


The  linear  models  generated  by  the  engine  manufacturer  are 
11th  order  with  9  control  inputs  and  22  outputs.  The  state  vari¬ 
ables,  controls  and  outputs  used  are  listed  in  Tables  5.3,  5.4  and 
5.5  respectively.  In  Figure  5.3,  a  schematic  of  the  engine  is  pre 
sented  showing  the  physical  locations  of  the  controls  and  output 
quantities. 
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Table  5.3 


Engine  State  Variables 


NUMBER 

VARIABLE 

SYMBOL 

1 

Front  Fan  Percent  Corrected  Speed 

PCN2 

2 

HP  Compressor  Percent  Corrected  Speed 

PCN25 

3 

HP  Compressor  Metal  Temperature 

TM3 

4 

HP  Turbine  Metal  Temperature 

TM41 

5 

HP  Turbine  Blade  Temperature 

TB41 

6 

Bypass  Duct  Gas  Weight  Per 

Unit  Volume 

RH15B 

7 

Bypass  Duct  Gas  Entropy  Per 

Unit  Volume 

SRH15B 

8 

Main  Combustor  Gas  Weight  Per 

Unit  Volume 

RH31 

9 

Main  Combustor  Gas  Entropy 

Per  Unit  Volume 

SRH31 

10 

Tailpipe  Gas  Weight  Per  Unit 

Volume 

RH6 

11 

Tailpipe  Gas  Entropy  Per 

Unit  Volume 

SRH6 

Table  5.4 
Engine  Inputs 


NO. 

VARIABLE 

SYMBOL 

1 

Second  Fan  IGV  Stator  Angle 

STP22 

2 

Exhaust  Nozzle  Throat  Physical  Area 

A8 

3 

Main  Combustor  Effective  Fuel  Flow* 

WF36 

4 

IP  Turbine  Nozzle  Position 

STP49 

5 

Inner  Ouct  Mixer  Exit  Physical  Area 

A14B 

6 

Bypass  Ouct  Mixer  Effective  Area 

AE16 

7 

Augmentor  Effective  Fuel  Flow 

DWF6 

8 

High  Compressor  Stator  Angle 

STP25 

9 

Forward  Blocker  Door 

AE94 
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Table  5.5 

Engine  Output  Variables 


NUMBER 

VARIABLE 

SYMBOL 

1 

Front  Fan  Percent  Corrected  Speed 

PCN2 

2 

HP  Compressor  Percent  Corrected  Speed 

PCN25 

3 

Front  Fan  Exit  Mach  Number 

XM93 

4 

Second  Fan  Exit  Mach  Number 

XM23A 

5 

Front  Fan  Exit  Total  Pressure 

P21 

6 

HP  Compressor  Inlet  Total  Pressure 

P25 

7 

HP  Compressor  Discharge  Static  Pressure 

PS3C 

8 

Front  Fan  Discharge  Tip  Total  Temperature 

T93 

9 

HP  Compressor  Inlet  Total  Temperature 

T25 

10 

HP  Compressor  Discharge  Total  Temperature 

T3 

11 

HP  Turbine  Rotor  Inlet  Total  Temperature 

T41 

12 

HP  Turbine  Blade  Temperature 

TB41 

13 

HP  Turbine  Discharge  Total  Temperature 

T48 

14 

LP  Turbine  Discharge  Total  Temperature 

T5 

15 

Exhaust  Nozzle  Inlet  Total  Pressure 

P7 

16 

Front  Fan  Stall  Margin  at  Constant  Flow 

SM2 

17 

Second  Fan  Stall  Margin  at  Constant  Flow 

SM22 

18 

HP  Compressor  Stall  Margin  at  Constant  Flow 

SM25 

19 

Front  Fan  Inlet  Air  Flow 

WA2 

20 

Bypass  Duct  Total  Air  Flow 

Ml  53 

21 

Net  Thrust 

FN 

22 

Uninstalled  Specific  Fuel  Consumption 

SFC 

Figure  5.3  Variable  Cycle  Engine  -  Controls  and  Outputs 
5*1.3  Model  Reduction 

the  AS.deSCribed  in  Section  S'l-2.  the  linear  models  generated  bi 
the  engine  manufacturer  are  nth  ordor.  Using  lincar  analysis 

“it  :  U,  P°SSibU  t0  id0"tify  th0Se  ”odos  »f  these  models 

controller  t  "h”"’  th°  d°Slr0d  bandwldth  of  a  ‘los°d  lo»P 
...  *  S  hon  Possible  to  generate  reduced  order  models 

law'd  T  °nl>’  the  dynMical  Ihtoractions  relevant  to  control 
w  development.  By  eliminating  "unimportant"  dynamics,  this 

fimnortlVr  (1)  yield  ”°d0lS  that  ar®  in  their  si»Pl»st  form 

tational  \  8aini"8  P'’>'SiCal  inSi8ht)-  ond  <2>  reduce  the  compu 
tational  complexity  required  in  control  law  calculations. 

One  point  must  be  remembered,  however.  Any  control  law 
evelopod  using  these  reduced  order  models  must  be  verified  when 

r;dr  l"\tUU  ^  «*  U  possible  the  cent", I”, 

affect  the  higher  modes  -  perhaps  leading  to  an  instability.  A 

check  with  the  full  nonlinear  model  (e.g.,  hybrid  simulation)  is 
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also  required  to  assure  that  no  adverse  effects  arise  due  to 
unmodeled  nonlinearities. 

5.1. 3.1  Eigensystem  Analysis 

The  objectives  of  an  eigensystem  analysis  are;  (1)  to  iden¬ 
tify  the  modes  of  the  linear  model  that  are  within  the  bandwidth 
of  the  desired  closed  loop  system  (0  to  20  sec'*),  and  (2)  to 
identify  the  state  variables  which  dominate  those  modes.  This 
information  forms  the  basis  for  the  model  reduction  procedure 
employed  in  Section  5. 1.3. 3. 

Eigenvalues  for  the  eleventh  order  full  power,  sea-level- 
static  linear  model  are  presented  in  Table  5.6.  (These  eigen¬ 
values  have  been  normalized  by  the  time  constant  of  the  collective 
response  mode  of  the  spools).  A  clear  frequency  separation  is 
apparent.  Only  five  of  the  eleven  modes  are  in  the  frequency 
range  of  interest  of  the  controller.  These  modes  are  isolated 
in  Table  5.7. 

The  rotor  response  is  dominated  by  a  second  order  system  with 
two  real  roots.  Examination  of  the  eigenvectors  (Figure  5.4) 
shows  that  the  faster  response  can  be  associated  with  the  rotor 
speeds  differentially  rematching,  i.e.,  the  fan  spool  decelerating 
while  the  compressor  accelerates.  The  slower  root  involves  the 
collective  or  rigid  body  response  of  the  two  shafts.  The  differ¬ 
ential  and  collective  response  is  characteristic  in  the  motion  of 
two  heavily  damped  and  heavily  coupled  inertia  elements.  (The 
two  spools  are  aerodynamical ly  coupled  in  the  turbines.)  The 
characteristic  motion  is  present  in  all  regions  of  the  operating 
envelope.  The  three  temperature  response  modes  shown  in  Table  5.7 
are  essentially  first  order  lags  resulting  from  the  thermal 
capacitance  and  resistance  of  the  engine  components. 
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Eigenvalue  locations  are  functions  of  flight  condition. 

However,  the  frequency  separation  demonstrated  above  is  maintained 
throughout  the  flight  envelope.  Figures  5.4  through  5.7  demon¬ 
strate  this.  Figure  5.4  is  a  display  of  the  collective  and  differ¬ 
ential  response  modes  of  the  spools  for  the  first  six  flight  con¬ 
ditions.  Figure  5.5  displays  the  same  information  for  the  three 
temperature  roots  associated  with  TM3,  TM41  and  TB41.  Figures  5.6 
and  5.7  display  the  eigenvalues  associated  with  the  remaining  six 
gas  dynamic  modes.  These  are  clearly  much  faster. 

Note  that  all  the  eigenvalues  have  been  normalized  by  the 
collective  response  mode  time  constant  at  sea-level-static,  full-power. 

The  state  variables  which  dominate  each  mode  can  be  found  by  an 
examination  of  its  associated  eigenvectors.  In  order  to  reduce  the 
system  order,  an  equal  number  of  states  and  modes  must  be  chosen. 
Furthermore,  the  states  chosen  must  be  able  to  describe  the  re¬ 
tained  modes  if  the  model  reduction  procedure  of  Section  3.6.1  is 
to  be  possible. 
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Table  5.6 

Eigenvalues  of  VCE  at  SLS  -  Maximum  Power 


NUMBER 

VALUE  * 

DESCRIPTION 

1 

-0.13 

Compressor  Metal  Temperature  Response 

2 

-0.39 

Turbine  Metal  Temperature  Response 

3 

-0.47 

Turbine  81ade  Temperature  Response 

4 

-1.00 

Collective  Response  Mode  of  Spools 

5 

-3.24 

Differential  Response  Mode  of  Spools 

6.7 

-22.7i;2.83 

Has  Oyn 

antics 

3.9 

-60.75U9.30 

10 

-65.35 

11 

-131.22 

' 

•Norm*] Ued  by  Collective  Response  Node  Eigenvalue 


Table  5.7 

Principal  Dynamic  Modes  of  VCE  at  SLS  -  Maximum  Power 


MODE 

TIME 

CONSTANT* 

DESCRIPTION 

1 

0.271 

Differential  Response  Mode  of  the  Spools 

2 

1.000 

Collective  Response  Mode  of  the  Spools 

3 

6.368 

Compressor  Metal  Temperature  Response 

4 

2.206 

Turbine  Metal  Temperature  Response 

5 

1.975 

Turbine  Blade  Temperature  Response 

‘Normalized  by  collective  response  mode  time  constant. 
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Figure  5.4  Differential  and  Collective  Spool  Response  Time 
Constants  at  Flight  Points  1  Through  6 


Figure  5.S  Temperature  Response  Time  Constants  at 

Flight  Points  1  Through  6 
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Figure  5.6  Eigenvalues  Corresponding  to  Engine  Gas  Dynamic 

at  Caeh  of  Six  Flight  Conditions 
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Figure  5.7  Eigenvalues  (Gas  Dynamics)  Vs.  Flight  Conditions 
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Results  are  presented  in  Figure  5.8  and  Table  5.8.  The 
figures  picture  the  eigenvectors,  and  the  table  lists  the  states 
chosen. 


(a)  eigenvector  of  Made  1  (Differentia!  Spc*!  response) 


(b)  Eigenvector  of  Mode  2  (Collection  bocal  Pesocnte) 


-1.0  1*0 
(c)  Eigenvector  of  Mode  3  (Cowpmior  Metal  Teeoeroture) 


td)  Eigenvector  of  Mode  *  (Turoine  Metal  ’*K*rnur*> 


(e)  Eigenvector  of  Node  $  (Turoine  Slade  ftwperatore) 


Figure  5.8  Eigenvector  Ootapos  i  t  ion  of  the  Five  Principal  Modes 
at  SLS/Maxitatusi  Power  (Table  5.7) 


Table  S.8 

State  Variables  for  Reduced  Order  System 


NUMBER 

VARIABLE 

SYMBOL 

1 

Front  Fan  Percent  Corrected  Speed 

PCN2 

2 

HP  Compressor  Percent  Corrected  Speed 

PCN25 

3 

HP  Compressor  Metal  Temperature 

TM3 

4 

HP  Turbine  Metal  Temperature 

TM41 

5 

HP  Turbine  Blade  Temperature 

TB41 

5.1. 3.2  Transfer  Function  Analysis 


The  eigensystem  analysis  above  is  most  useful  in  determining 
the  system  modes  which  might  be  important  in  the  controller  design. 
It  is  entirely  possible,  however,  that  some  of  these  modes  may 
be  ignored,  'lhat  is,  they  may  not  contribute  significantly  to 
the  output  quantities  which  are  to  be  controlled.  This  subject 
is  most  easily  addressed  using  transfer  function  analysis  tech¬ 
niques  . 

Transfer  function  techniques  examine  the  structure  of  the 
input-output  relations.  They  rely  on  an  investigation  of  the 
system  zaroes,  residues  and  dc  gains. 

The  transfer  function  from  an  input,  u(s),  to  an  output,  y(s), 
may  be  written  as 

y (s)  (s+tj)  (s  +  z2)  » ♦ » 
u(sj^  Ts+PjUs+Pj)  -r* 


or,  equivalently,  as 


y(s)  _  v  ^  V  JPj  .  4.  k'  -El-  4.  4.  If - 5 

Sfit  "  KDF  K1  s+Pl  k2  s<p2  *•*  KN  s+i 


Here.  KDp  is  the  direct  feedthrough  component,  and  the  define 
the  various  modal  components  of  the  dc  response  (K^p^  is  the  ■ 
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residue  at  the  i**1  pole).  An  examination  of  the  relative  magni¬ 
tudes  of  these  and  K.  will  determine  those  modes  which  must  be 
retained  for  an  accurate  description  of  the  engine’s  output 
response. 

Figure  5.9  summarizes  an  analysis  of  transfer  functions  at 
flight  conditions  one  through  six.  Shown  are  K^p  and  the  foi 
PCN2  PCN2  FN  FN 

~WF~*  A8  aiU*  WF‘  T*‘e  m0<*es  indicated  were  summarized  in 

Table  5.7.  In  the  figure  refers  to  the  differential  response 
mode  of  the  spools;  Ac  refers  to  the  collective  response  mode;  A^ 
refers  to  the  comprecsor  metal  temperature  mode;  and  A41  refers  to 
the  turbine  metal  temperature  mode. 


(•)□□□□□ 


Figure  5.9  Residue  Spectra  of  Selected  Transfer  Functions 
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The  results  shown  in  Figure  5.9  are  characteristic.  Several 

of  the  transfer  functions  are  dominantly  low  order,  but  every  mode 

FN 

is  important.  For  example,  —  is  zeroth  order  (direct  feedthrough 

Ao 

nnv| 

dominates)  at  all  flight  points;  and  ■  is  first  order  at  flight 

A8 

conditions  1  and  2,  but  second  order  at  flight  condition  4.  The 
modal  composition  of  each  transfer  function  is  different;  how¬ 
ever,  every  mode  contributes  significantly  to  at  least  one  trans¬ 
fer  function. 

The  turbine  blade  temperature  mode  is  excited  by  the  con¬ 
trols,  but  is  observable  in  TB41  only.  The  reason  this  mode 

must  be  retained  is  that  turbine  blade  temperature  is  important  to 
engine  protection  constraints. 

The  conclusion  which  must  be  drawn  is  that  none  of  the  five 
modes  of  Table  5.7  can  be  ignored  in  an  accurate  description  of 
the  engine. 

5. 1.3. 3  Modal  Reduction 


Based  on  the  above  analysis,  a  fifth  order  model  of  the  en¬ 
gine  which  incorporates  the  five  modes  listed  in  Table  5.7  and  the 
five  states  listed  in  Table  5.8  will  characterize  the  engine’s 
dynamic  behavior  sufficiently  well  for  control  purposes.  The 
procedure  for  generating  this  reduced  order  model  is  highly  auto¬ 
matic.  It  was  explained  in  Section  3.6.1. 

During  Phase  I  of  this  program  28  linear  models  were  generated 
by  the  engine  manufacturer  and  processed  at  SCI  (Vt). 


5.1.4  Second  Stage  Reduction 


As  explained  in  Section  3.6.2,  not  all  of  the  elements  in  the 
sy.Sem  dynamics  matrix,  F,  are  important.  Specifically,  the  sensi¬ 
tivity  of  the  dynamic  response  (initial  conditions)  to  a  particular 
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element  may  be  small.  If  so,  that  element  could  be  fixed  at  zero 
with  little  effect.  Furthermore,  the  remaining  elements  can  be 
adjusted  to  minimize  any  error  that  does  occur.  A  procedure  for 
performing  these  calculations  is  presented  in  Appendix  B. 

Before  elements  in  the  F  matrix  can  be  eliminated,  however, 
it  is  first  required  to  determine  the  response  sensitivities. 

This  can  be  done  as  follows.  Let  the  integral  square  response  be 
defined  as 


J  =  i  /  xT  Av  x  dt. 

0  Z  t  x 


where 


and 


x  =  Fx 
x(o)  =  XQ 


with 

e[xq  xj]  -  X  (X  »  I  for  random  initial  conditions) 


3J0 

The  object  is  then  to  find  -gy  . 

This  is  easily  accomplished.  Consider  the  following  substi 
tions: 


Geq  ■  1 

v  ■ 1 


c„„  -  4 r. 

eq 
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3J. 


Then  finding 


3C 


eq 


for  the  system 


x  =  Fx  +  G  u 

eq  eq 


y  =  H  x 
1  eq  eq 


u  =  C  y 
eq  eq  'eq 


is  the  same  as  finding 


3J_  3J. 


3AF  3F 


aF=0 


in  the  system 


x  =  Fx  +  AFx 
x  =  (F+AF)x. 


In  other  words,  the  output  regulator  algorithm  discussed  in 
Section  3.5.3  may  be  used. 


The  values  of  have  been  calculated  throughout  the  flight 
dr 

envelope,  and  the, results  are  presented  in  Figure  5.10.  The  figure 
is  arranged  in  a  matrix  format  so  that  the  diagram  in  the  (1,1) 

9J 

block  of  the  figure  corresponds  to  ©tc.  Those  elements  which 

dtll 

may  be  eliminated  are  indicated  by  shading  in  the  figure. 
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MAXIMUM  SENSITIVITY \QVERAU 
LINEAR 

MINIMUM  SENSITIVITY  I*0051-5 


Figure  5,10  Sensitivity  of  Mean- Square  Response  to  Elements  in 
F-Matrix 
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Example: 


Here  the  second  stage  reduction  procedure  is  applied  to  a 
fifth  order  linear  model  of  the  GE23-JTDE.  The  model  chosen  has 
the  following  F  matrix: 


F  = 


-  5 . 6600E+00 
2 . 6860E+00 
6. 8520E-02 
6. 4180E- 03 
1 . 6360E- 01 


4 . 6230E+00 
-  8 . 8090E+00 
2 . 6050E-01 
-1 . 0110E+00 
-1 . 6780E-01 


1 . 1860E-01 
1.8910E-01 
-4. 6350E-01 
1 . 0520E- 01 
3. 2950E-01 


2 . 3240E- 01  0. 

4. 2300E-01  0. 

9. 0940E-03  0. 

- 1 . 2150E+00  0. 

9 . 7680E-02  -1.6220E+00 


It  is  desired  to  eliminate  as  many  elements  as  possible, 
yet  still  retain  an  acceptable  fit  to  the  original  response.  Step  1 
of  the  procedure  is  to  identify  the  sensitivity  of  the  system 
response,  JQ,  to  each  element  in  the  F  matrix.  For 

J0  =  *Tx  dt 


and 


I  (random  initial  conditions) 


the  result  is 


dJ, 


dF 


lii 


ij 


l.lxlO'1  -5.3xl0*2 
»-  2 


-6.5  xio 


-3 


3.7x10  6  9.7x10  6  -1.1  xlO 


-S.5xl0*3  -1.5xl0'2 

,-S 


.-1 


-6.5x10' 

4.2x10 


6.0  xio  A  -7.1x10 


-3 

-4 


5.2x10  3  5.8x10  J  -7.81x10 


.-3 


.-4 


-2 


2.0x10 


1.4x10  J  9.1x10  H  -3.68x10  6  -1.4x10 


0 

0 

0 

0 

1.9x10' 


Clearly,  the  diagonal  elements  must  not  be  eliminated.  For  example, 
setting  ^33*0  would  cause  a  601  change  in  JQ.  However,  should 
be  eliminated  since  this  would  produce  only  a  0.07%  change  in  JQ. 
Using  this  logic,  the  elements  chosen  for  elimination  are 


F13*  F14*  F24»  F31 *  F34»  P41»  F42*  F43»  F51* 


F52*  FS4  * 
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This  leaves  ten  non-zero  elements  which  must  be  optimized. 
The  results  of  this  optimization  are 


F  = 


-5.7207 

2.7746 

0 

0 

0 


4.7626 
-9.1784 
0. 3294 
0 
0 


0 

0.3494 

-0.4651 

0 

0.3357 


0 

0 

0 

-1.2885 

0 


0 

0 

0 

0 

-1.6212 


Note  that  the  original  JQ  was 

JQ  =  1.03 


Setting  the  indicated  elements  to  zero  and  performing  no  optimiza¬ 
tion  yielded  only  a  1.0%  change  in  JQ.  Performing  the  optimization 
reduced  this  already  small  error  to  0.9%.  This  confirms  the  "unim¬ 
portance"  of  those  terms  which  were  eliminated  and  thus  provides  a 
much  simplified  F  matrix  for  use  in  the  controller. 

The  eigenvalues  of  the  original  and  reduced  systems  are 
presented  for  comparison  in  Table  5.9. 


Table  5.9 

Comparison  of  Eigenvalues*  between  F  and  F 


ORIGINAL 
EIGENVALUES  (F) 

REDUCED  A 

EIGENVALUES  (F) 

-0.49 

-0.49 

-3.38 

-3.50 

-1.00 

-1.05 

-0.41 

-0.39 

-0.14 

-0.14 

♦Eigenvalues  normalized  by  collective  response  mode 
eigenvalue. 
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5.2  NONLINEAR  MODEL  GENERATION 


The  reduced  parameter  models  developed  from  the  nonlinear 
simulation  have  been  described.  The  procedure  for  developing  a 
nonlinear  model  of  the  engine  from  these  matrices  was  described  in 
Section  3.7.  Twenty-eight  linear  matrices  have  been  used  to 
generate  a  low  order  nonlinear  model  of  the  VCE  valid  throughout 
the  operating  envelope.  These  preliminary  results  are  presented 
below  and  a  comparison  of  the  linear  and  nonlinear  response  is 
shown. 


5.2.1  Regression  for  Dynamic  Parameters 


The  model  generation  procedure  is  shown  in  Figure  5.11. 


Figure  5.11  Nonlinear  Model  Generation  Procedure 
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The  reduced  parameter  models  discussed  in  Section  5.1  repre¬ 
sented  fifth  order  engine  dynamics  and  were  parameterized  with 
10  nonzero  elements.  Curves  were  derived  for  each  dynamic  element 
using  subset  selection  methods.  Independent  variables  were  chosen 
as  transgenerated  groupings  of  engine  running  variables.  A  typical 
set  of  operating  variables  was  initially  established  by  using  the 
stepwise  regression  procedure  on  a  large  set  of  candidate  terms 
(e.g.,  80-100).  The  variables  selected  during  this  procedure  are 
then  used  as  the  generators  for  a  regression  using  the  optimal 
subset  selection  method.  This  approach  limits  the  number  of 
25  regressed  terms  because  the  computational  overhead  increases 
exponentially  with  the  number  of  terms. 

The  models  derived  for  the  elements  of  the  dynamics  matrix 
form  the  basis  for  the  nonlinear  model.  Table  5.10  shows  the 
preliminary  model  selected  to  represent  engine  dynamics  throughout 
the  operating  envelope.  Thirty-three  (33)  terms  are  sufficient  to 
match  the  linearized  dynamics  at  all  flight  points.  The  model's 


Table  5.10 

Preliminary  Nonlinear  Model  Structure 


MATRIX 

ELEMENT 

MODEL  FORM 

RMS 

ERROR* 

Pll 

alPs3c  +  a2P2  +  a3 

1.0 

f33 

blps3c  +  b2N2T2  +  b3 

0.3 

f12 

clp2  +  c2t2  +  C3N25P2  +  c4 

1.7 

f21 

<*lp2  +  d2p21  +  d3ps3c  +  d4 

1.5 

f22 

elFll  +  e2p|  +  ®3n2p2  +  e4 

1.2 

f23 

flp2  +  f2p21  +  F3F33  +  f4 

1.5 

f32 

9lp2  +  92p21  +  93ps3c  +  94 

1.2 

f44 

hiF33  +  hg 

0.3 

f53 

klPS3c  +  k2F33  *  k3 

0.4 

f55 

ijF 33  +  *2 

0.1 

*  Percentage  of  mean 
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accuracy  is  summarized  in  Figure  5.12,  which  shows  a  comparison  of 
the  actual  engine  time  constants  and  those  calculated  from  the 
curve  fit  model  of  the  dynamics  matrix.  This  procedure  represents 
a  significant  simplification  over  other  linear  modeling  methods 
and  allows  a  tractable  simulation  of  nonlinear  dynamics. 


(RAD/SCC) 


Figure  5.12 


0  I  i  3  4  b 


AC1UAI  COUECTIVE  RUOl 
(RAH/ SEC) 


Comparison  of  Linear  Eigenvalue  Locations  and  Root 
Locations  Found  by  Linearizing  the  Nonlinear  Model  at 
28  Flight  Points  Spanning  the  Operating  Envelope. 

The  Three  Remaining  Temperature  Modes  Had  Less  Error 
Than  the  Collect iye/Differential  Response  Models 
Shown  Above 


The  nonlinear  dynamics  matrix  developed  from  the  linear  models 
was  used  in  a  nonlinear  simulation  of  the  engine  response  for  a 
large  acceleration  from  idle  to  military  power.  A  step  input  of 
the  control  variables  was  used  in  this  test  because  the  definitions 
of  the  control  logic  governing  such  a  maneuver  has  not  been 
developed.  The  response  of  the  nonlinear  simulation  is  compared 
to  a  linear  model  derived  at  intermediate  power  and  a  model  which 
uses  a  linearization  at  the  trajectory  point.  (Figure  5.13)  It 
is  observed  that  the  response  of  the  nonlinear  equations  falls 
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CflHPfiESSOR  SPEED  FAN  SPEED 

{X  DESIGN)  DESIGN) 


Figure  S.13  Comparison  of  Idle  to  Military  Power  Step  Response 
Simulated  Using  Linear  Model  Derived  at  Intermediate 
Power,  the  Nonlinear  Model  Derived  from  linear  Matrices 
and  a  Model  Using  Linearizations  at  Each  Point  along 
the  Trajectory 
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between  the  fixed  linear  and  pointwise  linear  simulation.  Further 
comparison  of  the  approximate  nonlinear  response  to  the  detailed 
digital  simulation  will  be  made  during  the  next  phase  of  the 
program. 

5.3  ACTUATOR  SERVO  VALVE  COMPENSATION 

As  part  of  Phase  I,  compensation  for  the  actuator  servo  valves 
has  been  designed  and  a  detailed  simulation  performed.  The  purpose 
of  the  simulation  was  to  determine:  (1)  the  sensitivity  of  the 
proposed  compensation  to  known  errors  (e.g.,  bias  and  gain),  and 

(2)  the  effects  of  variable  sample  and  computation  intervals  (i.e., 
the  impact  of  not  having  a  clock). 

Two  systems  have  been  simulated.  They  are  based  on  the 
original  and  revised  versions  of  the  servo  valve  characteristic 
respectively.  These  characteristics  are  reproduced  in  Figures  S.14 
and  S.1S. 

The  control  design  specifications  for  both  systems  were 
assumed  to  be  the  same.  Specifically  required  were: 

(1)  bandwidth  of  10  rad/sec  (linear  operating  region) 

(2)  linear  response  (no  saturation)  for  stlOl  of  full  scale 
step  changes  in  commanded  fuel  flow 

(3)  dc  hangoff  error  of  0.11  stroke 
5*3.1  Simulation 

The  simulation  performed  involves  a  detailed  description  of 
the  loop  error  sources  and  nonlinearities.  Included  are  the 

(1)  valve  characteristic  (nonlinear) 

(2)  servo  valve  time  constant 

(3)  servo  valve  gain  error 

(4)  variable  bias  (s2  ma) 
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VALVE  FLOW 


Figure  5.14  Failed- Fixed  Servo  Valve  Characteristic 


Figure  S.1S  Revised  Failed-Pixed  Servo  Valve  Characteristic 
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(5)  2  ma  hysteresis 

(6)  additive  white  measurement  noise 

(?)  8  bit  quantizer 

(8)  12  bit  sampler 

(9)  variable  loop  sample  time 

(10)  variable  control  computation  time 

(11)  427.25  Hz  pulsewidth  modulator 

A  block  diagram  showing  these  terms  is  presented  in  Figure  5.16. 

Note  that  the  variable  times  and  modulation  interval  introduce 
effective  time  delays  (variable)  to  the  system.  These  delays  and 
the  hysteresis  are  sources  of  phase  lag  which  adversely  affect 
loop  stability. 

The  simulation  was  designed  specifically  to  mimic  operation 
of  an  actual  digital  control  loop.  Implementation  of  the  variable 
sample  interval,  modulation  time  (23.5  ms),  and  control  computa¬ 
tion  time  is  described  with  the  aid  of  Figures  5.17  and  5,18. 

Figure  5.17  describes  the  sequence  of  events  in  time.  Note  that 
the  sample  intervals  shown  are  not  constant  (this  is  intended  to 
describe  operation  without  a  clock).  The  sample  interval  is  assumed 
to  be  randomly  distributed  between  5  and  IS  ms.  At  the  beginning 
of  each  sample  interval  the  computer  memory  location  containing 
the  LVPT  output  is  sampled.  However,  this  memory  location  is  up¬ 
dated  only  every  2  ms.  Consequently,  since  the  sample  interval  is 
not  synchronized  with  this  update  interval,  the  information  con¬ 
tained  in  memory  will  lag  the  actual  LVPT  output.  This  delay  is 
evenly  distributed  between  0  and  2  ms. 

After  sampling  the  LVPT  memory  location,  the  computer  must 
compute  the  commanded  control  input.  This  takes  time.  Currently, 
this  time  interval  is  assumed  to  be  randomly  distributed  between 
5  and  15  ms. 
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:d/fixed 

LTERISTIC 


SL0PE=343pph/% 


INITIALIZE 
TIME  »  0 
TSAMP  »  T SAMP I 
TMOO  «  0 
TCOMP  *  0 


TSAKP'TSAMPI 


THOD-TMOD1 


TCOHP>TCOMPI 


(1!  SAMPLE  LVPT  MEMORY  LOCATION 

(2)  UOLD-OID  COMMANDED  CONTROL 

(3)  COMPUTE  NEW  COmANDED  CONTROL 

(4)  UNEW*NEW  COMAANDED  CONTROL 

(5)  CHOOSE  NEXT  SAMPLE  INTERVAL  (TSAMP1) 

(6)  TSAMP*0 

(7)  COMPUTE  TIME  UNTIL  END  OF  MODULATION 
INTERVAL  (TLEFT) 

(8)  TCOMP  •  0 


.YES  (1)  UPDATE  APPLIED  CONTROL:  U.«U„ 

> - *  (2)  TMOO'O  * 

(3)  COMPUTE  TIME  UNTIL  NEXT  SAMPLE 
INTERVAL  (TTGOJ  _ _ 


COMMANDED  CONTROL 
EQUALS  UNEW: 
U,*UNEW 


COMPUTE  NEW 
INTEGRATION  TIME  (DT) 
BASED  ON  TLEFT 


COMPUTE  NEW  INTEGRATION 
TIME  (DT)  BASED  ON  TTGO 


COMMANDED  CONTROL 
EOUAIS  UOLD: 

U  «UOLD 


INTEGRATE  (ONE  TIME  STEO 
i ■ f(x,uA,») 
y  •hU,uA,t.) 


TIME  •  TIKE  ♦  CT 
TSAMP  •  TSAMP  ♦  DT 
TMOO  •  TMOO  ♦  DT 
TCOMP  ■  TCOMP  ♦  DT 


Figure  S.18 


Flowchart  of  Actuator  Valve  Servo  Loop 
Simulation 
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After  the  commanded  control  input  has  been  calculated,  it  is 
applied  to  the  pulsewidth  modulator.  This  modulator,  however,  only 
outputs  a  new  pulse  at  23.5  ms  intervals.  Consequently,  there  is 
a  further  delay  (up  to  23.5  ms)  before  the  actuator  sees  the  new 
control.  Note  that  it  is  possible  for  a  second  commanded  control 
to  be  computed  before  the  first  has  been  applied  to  the  torque 
motor. 


Figure  5.18  is  a  flow  chart  of  the  computer  code  which  simu¬ 
lates 
are: 

Sample  Interval 

Running  Sum  for  Comparison  with  TSAMPI 

Modulation  Interval  (23.5  ms) 

Running  Sum  for  Comparison  with  TMODI 

Computation  Interval 

Running  Sum  for  Comparison  with  TCOMPI 

je  Kutta  integration  is  used.  The  intej 
time  step  is  varied  during  execution  to  assure  that  the  end  of  each 
sample  interval  and  the  end  of  a  modulation  interval  both  coincide 
with  an  integration  time  (variables  TLEFT  and  TTGO) . 


the 

process 

(1) 

TSAMPI: 

(2) 

TSAMP: 

(3) 

TMODI : 

(4) 

TMOD: 

(5) 

TCOMPI: 

(6) 

TCOMP: 

Fourth  order 

5.3.2  FADEC  Valve  I 


The  preliminary  compensation  designed  for  the  failed-fixed 
servo  valve  (Figure  5.1  0  is  presented  in  Figure  5.19.  This  com¬ 
pensation  is  essentially  a  stabilized  inverse  dead  zone.  With  this 
compensation,  the  10  rad/sec  bandwidth  requirement  is  easily  satis¬ 
fied.  However,  in  a  worst-case  error  configuration  (hysteresis  and 
bias),  the  dr  hangoff  specification  is  not  met. 

I 

The  design  proceeds  as  follows.  In  the  linear  operating 

region,  two  approximations, are  first  made:  (1)  the  failed/fixed 

*  ■  » 
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Figure  5.19  Nonlinear  Compensation  for  Actuator  Servovalve 
(FADEC  I) 

characteristic  is  roughly  a  constant  gain  (101/sec/ma) ,  and  (2)  the 
compensation,  Hc>  is  a  constant  gain  (KK^) .  The  open-loop  transfer 
function  is,  therefore, 

KK^IO) 

G  =  s  (0 . 025s+l )" 

and  the  closed-loop  transfer  function  is 

piston  position  400  KK^ 

-  =*  -j - 

commanded  position  s  +  40s  +  400  KK^ 

The  closed-loop  pole  locations  are,  therefore, 

sx  -  -  20  -  20  /r"-'TK1 

S2  »  -20+20  A  -  KKX 

The  bandwidth  requirement  is  met  if 

S2  »  -  10  -  -  20  +  20  /T^THTj 

or 

KKj  -  0.75 

The  dc  hangoff  error  specification  (c50. 1 %)  requires  (no  bias  or 
hysteresis) 

ij  i  Ke 
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or 


6X  i  0. IK 


A  set  of  parameters  which  satisfies  all  the  requirements  is 

Kx  =  1.00 

K  »  0.75 

«  0.075 

This  compensation  has  been  incorporated  in  the  detailed  Simula 
tion  discussed  above.  The  results  are  presented  in  Figure  5.20. 
Shown  are  two  responses.  In  the  first,  there  are  no  time  delays. 

In  the  second,  all  the  error  terms  in  the  simulation  were  included. 


Figure  5.20  Step  Response  of  FADBC  Valve  *  1  to  a  101  of  Full 
Scale  Step  in  Commanded  Stroke  (Measurement  Noise 
Equals  II  of  Full  Scale  Stroke) 
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The  destabilizing  effect  of  the  phase  lag  associated  with  the 
hysteresis  and  time  delays  is  evident.  Note  that  measurement  noise 
signal  to  II  of  the  full  scale  stroke  has  been  assumed. 

5.3.3  FADEC  Valve  II 


During  the  Phase  I  effort,  the  engine  manufacturer  revised  the 
servovalve  characteristic  to  that  presented  in  Figure  5.15.  This 
new  characteristic  does  not  have  the  large  deadband  of  the  FADEC 
Valve  I;  consequently,  it  may  be  compensated  using  linear  techniques 

For  this  valve,  both  the  dc  hangoff  error  and  bandwidth  require 
ments  are  easily  met.  The  dc  hangoff  error  requirement  is  satisfied 
by  specifying  a  lower  limit  on  the  dc  gain  of  the  compensation 
Hc(s).  Specifically,  a  dc  current  error  of  4  ma  is  possible  (2  ma 
in  bias  and  2  ma  in  hysteresis).  This  requires 

,  4  [ma]  „  ,  ..  , 

xss[%]  "  Hc(o)[ma/tf  =  0A[i] 

or 

H.(o)  =  40  ma/t. 

The  bandwidth  requirement  of  10  rad/sec  is  satisfied  using  a 
compensation  of  the  form 


Hc(s) 


HcC°)(*s  s+1) 


The  time  constants,  t  and  x  ,  are  easily  found  using  Bode  analysis 

s  p 

One  set  which  satisfies  the  constraints  is 


Hc(o)  -  40.0 

x  »  1.S9  sec 

t„  -  20.0  sec 
P 
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This  compensation  loop  has  been  simulated.  The  results  are 
presented  in  Figure  5.21. 


Figure  5.21  Step  Response  of  FADEC  Valve  #2  to  a  101  of  Full 
Scale  Step  in  Commanded  Stroke  (Measurement  Noise 
Equals  II  of  Full  Scale  Stroke) 

The  destabilizing  phase  resulting  from  the  time  delays,  hysteresis 
and  sampling  is  apparent.  Furthermore,  the  effect  of  the  noise* 
passing  through  the  linear  compensation  is  obvious  in  that  the 
final  valve  position  tracks  the  noise.  This  effect  was  not 
observed  in  the  FADEC  Valve  U  response  due  to  the  presence  of  the 
dead  zone  (see  Figure  5.20).  The  final  compensation  for  FADEC 
Valve  12  will  include  a  deadzone. 


5.4  TRANSITION  GENERATOR  DEMONSTRATION 


Trajectory  generator  logic  was  designed  and  demonstrated  at 
one  flight  condition  to  evaluate  the  concept  (see  Section  4.1.S). 

A  linear  model  of  the  engine  at  sea  levels,  static,  intermediate 
power  was  used  as  the  engine  simulator.  An  '’acceleration"  from 
901  intermediate  power  to  intermediate  power  was  performed  using 
fuel  flow  as  the  only  modulated  variable.  While  this  test  did  not 
exercise  the  full  flexibility  of  the  general  system,  it  served  to 
illustrate  feasibility  and  provide  a  prototype  software  implementa 
tion  for  the  final  system. 

The  response  of  the  engine  to  a  step  input  is  compared  to  the 
linear  servomechanism  response  in  Figure  5.22.  The  structure  of 
this  demonstration  system  is  shown  in  Figure  5.23.  The  compensa¬ 
tion  for  the  command  generator  will  be  designed  from  the  output 
regulator  processes  described  above.  It  is  observed  that  the 
compensated  command  generator  produces  a  faster  thrust  response  at 
the  cost  of  a  temperature  overshoot  and  unacceptable  surge  margin 
loss. 

The  nonlinear  variable  rate  limit  design  was  also  implemented 
using  the  T4  j  (turbine  entrance)  gas  temperature  limit  and  the 
compressor  surge  margin  limits  as  inputs  to  the  system.  The 
response  is  shown  in  Figure  5.24.  It  can  be  seen  that  the  trajec¬ 
tory  generated  with  this  system  provides  a  fast  response  without 
causing  a  predicted  overtemperature  or  unacceptable  surge  margin 
loss  in  the  model.  The  implementation  of  the  final  form  of  the 
system  will  include  several  limit  schedules  driving  variable  rates 
in  each  of  the  compensated  command  servos  used  to  produce  th?r 
nominal  control  and  output  time  histories. 
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Figure  S .22  Comparison  of  Step  Response  and  Linear  Servomechanism 
Response  for  a  Ten  (10)  Percent  Acceleration  to  SLS, 
Intermediate  Power 


Densons  t  rat  ion  Trajectory  Generator  for  VCE  Linear  Simulation  Using 
a  Fuel  Flow  Actuator  and  Two  Nonlinear,  Variable  Rate  Feedbacks 


SECTION  VI 

SUMMARY  AND  PROGRAM  STATUS 

6.1  SUMMARY 

The  development  of  a  multivariable  control  structure  for 
the  GE23-JTDE  variable  cycle  engine  has  been  described.  The 
variable  cycle  engine  provides  the  necessary  performance  flex¬ 
ibility  to  meet  advanced  mission  requirements  specified  for  the 
Grumman  design  623  fighter.  A  digital  control  processor  is  used 
in  conjunction  with  a  hydromechanical  backup  control  and  a 
variety  of  actuating  and  sensing  devices  to  produce  this  cap¬ 
ability.  The  specified  hardware  represents  a  major  step  for¬ 
ward  in  complexity  and  versatility  in  propulsion  system  design, 

A  multivariable  control  design  methodology  has  been  presented 
which  addresses  each  functional  requirement.  Linear  dynamic  models 
form  the  basis  of  the  design  methods.  Output  regulator  synthesis, 
failure  detection  and  accommodation  and  nonlinear  model  development 
are  presented  as  they  address  the  major  specifications  of  the  digital 
control.  Preliminary  design  results  are  described  in  regulator 
design,  actuator  compensation,  model  reduction  and  nonlinear  model 
development. 

The  multivariable  control  structure  has  been  defined.  A  refer¬ 
ence  point  schedule  uses  power  commands  and  ambient  conditions  to 
estimate  the  operating  point  quantities.  A  trajectory  generator 
produces  optimized  paths  between  operating  points  which  respect 
engine  limits.  An  output  regulator  provides  compensation  for 
modeling  errors  and  disturbances.  Failure  detection  and  accommoda¬ 
tion  is  provided  by  a  fault  tolerant  filter  algorithm  to  compensate 
for  arbitrary  sensor  failures;  engine  protection  logic  is  provided 
to  coordinate  transfer  to  the  hydromechanical  backup;  and  a  failure 
accommodating  actuator  compensation  module  provides  a  high  performance 


interface  to  the  GH23  fail-fixed  actuators.  The  specified  control 
structure  and  design  methodology  in  conjunction  with  the  sophisti¬ 
cated  engine  hardware  has  the  potential  to  realize  substantially 
improved  performance,  reliability  and  flexibility  in  the  next 
generation  propulsion  systems. 

6.2  PROGRAM  STATUS 

The  multivariable  design  program  for  the  GE23  is  being  per¬ 
formed  in  two  phases.  The  Phase  I  study  results  have  been  docu¬ 
mented  in  this  report.  They  represent  the  control  structure 
definition  and  an  analysis  of  linear  dynamics  of  the  engine. 

The  results  of  this  study  indicate  that  the  proposed  design 
methods  have  the  potential  for  quickly  and  efficiently  de¬ 
veloping  the  digital  multivariable  control  software  in  Phase  II 
of  the  program. 

The  development,  validation  and  evaluation  of  the  control 
design  will  be  undertaken  in  the  second  phase  of  the  program. 

The  nonlinear  simulation  will  be  used  as  the  test  bed  for  the 
control  logic  design.  The  developed  software  will  be  transferred 
to  the  engine  manufacturer  for  evaluation  on  a  hybrid  computer 
model  of  the  engine.  A  careful  analysis  of  the  results  will 
validate  the  system  for  future  test  in  prototype  engine  programs. 


APPENDIX  A 


DERIVATION  OF 


PROBLEM: 

Given  J  =  J(C(C)) 
and 

dJ 

dC 

find 

dJ 

IC 

where 

Z  =  CCI-DC)'1 
or 

C  «  (I+CD)'1^  . 

J  is  scalar 
C  is  mxp 
Z  is  mxp 
D  is  pxm 

Presented  first  is  a  summary  of  notations  described  in 
detail  in  Refs.  17  and  18.  Specifically 

Ei,  «  matrix  whose  (ij)  element  equals  one; 
rest  of  matrix  is  zero 

E  ■  matrix  whose  (i,j)tlt  block  *  E^ 

rsA  ■  row  string  of  A  *  {a,*  '  a-**--'  a  *] 
pxq  (Ixpq) 

csA  *  column  string  of  A  « 
pxq 

(pqxi) 


a*2 

• 

i 

i 


ai* 


i*"*1  row  of  A 
it^1  column  of  A 


Derivative  of  one  matrix  with  another 

(d/dB)  {A}  =  <dA/dbu>  kj  =  1 ,2 , . . .  ,s 
sxt  Pxc*  spxtq  1+  -  1,2 , . . .  ,t 

Kronecker  product 


A  = 
pxq 


:  B 
sxt 


<aijB> 

psxqt 


i|  °  1,2,. ..,p 
j =  1,2,.  •  •  >q 


The  derivation  proceeds  as  follows. 

From  J  =  J(C(C)),  chain  rule  differentiation  yields: 


dJ  _  n 

K  (I 


dJ 


m  d(rsC] 


)  ( 


d[rsC] 
dC 


^I> 


or 


but 


dJ  «  fi 

K 


dJ 


m  d{rs(Jj 


.  rs  41 

d(rs£]  dC 


)  ( 


dC 


therefore 


dJ  _  ....  dJ  x  ,  dfrs£]'  ^ 

ac  (*xrs  5  c  ar  } 

mxp  lxmp  mmpxp 


(A.l) 
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1  j; 


kLl  :*V  V  -*r~ 


is  known.  The  problem  is  to  find 
dC  dc 

First,  find 

dC 

1C 

From 

C  =  C(I-DC)  1 


H?  =  1C  [CCI-DC)’1]. 


Let 

B  =  (I-DC)'1 
pxp 

and 

F  =  (I-DC)  . 
pxp 


Then,  the  product  rule  yields 


dC 

3C 


*z 


nunxpp 


(I  X B)  +  (I  -s—sr C) 
p  m 


To  find  ^  use 


Using  the  product  rule 


[0). 


(A. 2) 
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Therefore 


or 


Now 


a!  ■  -  <r  — =»  >  a?  d^-F)'1 

m  pxp  p 


d  B 
3C 


(I 

m 


^ B  )  ac 

pxp  p 


a?B  3U  (I'DC)  =  c°3  -  ^  (DC) 


Using  the  product  rule  again 
dl>* 

1 


dF 

ac 


r° 

(i 

p 


c)  -  (I 

m 


D}  ac 


or 


dF 

ac 


(I  s»-er  D) 
m 


dC 

3C 


Inserting  Eq.  (A. 4)  into  Eq.  (A. 3) 

all  "  (r  ^  B)  (I  •*-*!»  S£  (lx  B) 
n  in  p 

Inserting  Eq.  (A. 5)  into  Eq.  (A. 2) 
d?  dC 

au  *  atr  ♦  (i^c)(i^b)(i^d)  4£ 

P  w  in  in  ®'* 


(A.  3) 


(A. 4) 


(A. 5) 


(I^C.B) 

P 
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which  can  be  written 


dC 

HC 


[I  +(I^C)(I^B)(I: 
mm  m  m 


W1  3§  <■' 


or 


dC 

HC 


d: 

m 


:A)  $  (I- 

P 


where 


A  =  I  +  CBD 


Noting  that 


dC  c 
Hf  t 
mmxpp 

th 


■B) 


where  the  (i,j)th  block  of  E  is 


Then  Eq.  (A. 6)  may  be  written 


E^  and  E  i 
mxp 


dC 

HU 


(csA)  (rsB) 


By  expanding  Eq.  (A. 7)  it  may  be  shown  that 


aigfl:  ■  (csa)— -bt 


By  substitution  of  Eq.  (A. 8)  into  Eq.  (A.  I) 


^  -  (I  >-irs  “  )(csA^BT) 
dC  n  dC 


B) 


(A. 6) 


s  m  by  p  blocks . 


(A. 7) 


(A. 8) 
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which  equals 


dJ 

3C 


at  "bt 

dC 


The  solution  is,  therefore, 


=  (I+CBD)1  —  (I-DC)'T 
aL  dU 


(A.  9) 


APPENDIX  B 

SECOND  STAGE  MODEL  REDUCTION 


Given  the  system 

x  »  Fx  x(o)  ®  xQ, 


(B.l) 


the  problem  of  "second  stage  reduction"  may  be  stated  as:  Find 
the  elements  of  a  matrix,  F,  (which  has  a  different  structure 
than  F)  such  that  the  response  of  the  system 
* 

x  -  F  x  x(o)  -  xQ  (B. 2) 


most  nearly  matches  the  response  of  the  system  defined  by  Eq.  (B>1). 
Specifically,  if  the  error  in  the  response  is  defined  as 

A 

2  =  X  -  x,  (B. 3) 


then  the  object  is  to  minimise 


J  •  1/2  f*  sT  Azz  dt. 


(B.4) 


For  example,  it  might  be  desired  to  find  an  F  of  the  form 


(B.S) 


which  most  nearly  matches  the  dynamic  response  of  the  original 
system 


(B.6) 
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That  is,  the  object  is  to  find  the  optimum  system  with  the 
element  fixed  at  2ero. 

A  convenient  procedure  for  solving  this  problem  follows. 
First,  define  a  matrix,  AF,  as 

AF  =  F-F  (B.7) 

Then,  from  Eqs,  (B.1)-(B.3) 

z  a  x  -  x  (B.8) 

*  Fx  •  F  x 
a  Fx  -  (F-AF) (x- z) 

=  Fz  -  AF(x-z) 


where  z(o)  =  0.  An  equivalent  statement  of  the  problem  now 
becomes : 


Find  AF  in 


*  P]‘P  3P]-[» -IP] 


(B.9) 


which  minimizes 


J  -  1/2  /“  txT2T] 


p  yp] 


dt 


(B. 10) 


In  the  example  used  above,  AF  would  be  of  the  form 


(B. 11) 


I 


and  the  problem  becomes  finding  6fn,  6f21,  and  <Sf22- 

An  interesting  observation  can  be  made  at  this  point.  If 
the  following  equivalences  are  made 


*  [:] 

(B. 12) 

Fe«  •  [o  3 

(B. 13) 

Ge,  ■  [“] 

(B. 14) 

%  ■  [o  '!] 

(B. 15) 

De,  ■  tO] 

(B. 16) 

ceq  ■  [«;o] 

(B* 17) 

’  [» ? 

j  (B. 18) 

■-“>  •  d  a 

(B. 19) 

then  the  optimisation  problem  is  identical  to  the 

output  regulator  problem  with  a  fixed  structure  gain  matrix. 

Specifically,  the  problem  is  finding  in  the  system 


eq  eq  eq  eq 


CB.2C) 


u  •  C  Y 
eq  *eq 
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1  vi'SpS 


jgsasss’serr-.-v 


and 


y 


eq 


(B. 21) 


with 

J  =  1/2  C  ylq  Ayoq  y  dt  (B.22) 

and 

E!xeq  sIq>  =  Xeq  (0)  (B. 23) 

Consequently,  this  problem  may  be  solved  with  the  same  algorithm 
used  for  the  calculation  of  output  regulator  gains. 


APPENDIX  C 

MODELING  PROCEDURES  FOR  THE  VARIABLE  CYCLE  ENGINE 


C.l  INTRODUCTION 

A  critical  part  of  the  multivariable  control  law  is  a  dy¬ 
namical  model  of  the  engine  response.  While  this  model  does 
not  predict  engine  response  exactly,  the  closer  to  the  real 
behavior,  the  better.  The  VCE  is  a  complex  engine  and  an 
accurate  mathematical  model  could  be  far  too  complex  to 
implement  on  the  digital  processor.  An  approach  to  model 
synthesis  is  described  below  which  appears  to  balance  accuracy 
and  complexity  considerations  in  a  straightforward  fashion. 

The  modeling  problem  has  been  solved  in  several  other 
applications.  It  is  useful  to  review  these  before  discussing 
the  VCE  in  particular.  Simulation  modeling  in  a  hybrid  computer 

environment  has  been  around  a  long  time.  Thermodynamics  and 
rig  data  are  combined  to  give  the  best  match  of  engine  response 
as  measured  on  prototype  engines  or  previous  models.  Simplified 
phenomenological  models  can  be  written  and  solved  in  the 
processor.  This  approach  is  the  basis  of  the  Spang-Corley  (1§1 
fault  detection  algorithm  for  the  QCSBE  engine.  In  this  case, 
maps  are  simplified  and  the  thermodynamic  equations  are 
incorporated  directly.  This  procedure  has  been  used  by 
several  people  to  model  static  engine  response  for  fault  monitor¬ 
ing  purposes.  Control  of  the  accuracy  and  complexity  trade-offs 
is  not  apparent. 

Linear  behavior  of  the  engine  is  used  to  design  the  control 
logic.  Direct  application  of  the  linear  dynamics  matrices  is 
not  generally  used  for  simulation  because  of  the  poor  steady 
state  response  of  the  models  to  moderate  inputs.  It  has  been 


IS! 


suggested  that  the  linear  dynamic  matrices,  along  with  accurate 
set  point  information,  could  be  used  to  simulate  the  nonlinear 
response  of  the  engine  in  a  more  accurate  fashion  than  using 
constant  coefficient,  linear  dynamical  equations.  This  pro¬ 
cedure  is  attractive  to  MVC  design  since  there  is  already  a 
reference  point  generator  requirement  in  the  synthesis  procedure. 

C .  2  COMPARISON  OF  LINEAR  AND  NONLINEAR  RESPONSE 

To  clarify  the  following  sections,  the  utilization  of 
linear  dynamics  matrices  for  simulation  is  compared  to  in¬ 
tegration  of  the  actual  nonlinear  equations.  Engine  dynamics 
can  be  represented  by  the  nth  order  set  of  nonlinear  equations: 

x  =  f(x,u,0)  (C. 1) 

y  »  h(x,u,0)  (C. 2) 

where  m  quantities  u  (controls)  and  q  inputs  0  (ambients) 
determine  the  states  and  state  rates,  x  and  x  respectively. 

These  quantities  appear  as  nonlinear  algebraic  functions  in 
the  output  equations  for  p  quantities,  y  of  interest.  Much  is 
known  about  the  steady  state  engine  response.  This  can  be 
written: 

0  3  fC*SS’  u>  9)  (c*3) 

y  *  h(xss,  u,  8)  (C.4) 

Most  of  the  modeling  effort  should  go  into  matching  steady 
state  response  since  it  is  at  these  points  that  engine  per¬ 
formance  is  usually  evaluated.  Dynamical  response  is  also 
important  to  the  overall  description  of  the  behavior  since 
these  dynamical  properties  determine  stability. 

A  linear  approximation  to  Eq.  (C.l)  can  be  evaluated  if 
derivative  matrices  of  f  are  available.  The  function  in 


Eq.  (C.l)  can  be  expanded  to  first  order  about  an  equilibrium 
point  (xQ,uo)  as  follows: 

a*  ■  fx(Vuo)Sx  *  fu<xo>uo)aa  (C.  5a) 

<Sx  *  x-xQ  (C.5b) 

6u  *  u-u„ 
o 

If  it  is  assumed  that  u(t)  is  piecewise  constant  on  an  interval 
T 

[n  ,(n+l)T],  Eq.  (C.5)  can  be  rewritten  relative  to  the  constant 
value  of  u  and  the  value  of  x  which  could  be  attained  if  u 
remained  at  that  value,  or, 

x  *  Fx(xss,  u(n}) (x(t)-xss)  (C. 6) 

where  xgs  must  satisfy  the  nonlinear  equilibrium  relationship 

0  *  f(xs$,  u(n) ,  8)  (C, 7) 

or,  equivalently, 

xs$  •  gCu(n),  8)  (C.8) 

where  g(u(n),  8)  came  from  the  reference  schedules. 

The  time  interval,  T,  can  be  made  small  relative  to  the 
dynamics  and  the  equation  used  to  solve  for  tho  response  for 
an  arbitrary  function,  u(t). 

The  principal  drawbacks  of  this  approach  are  the  accuracy 
and  the  storage  requirements  necessary  to  implement  Eq,(C.6).  The 
accuracy  limitation  arises  because  it  is  typically  true  that 
||x-xss!l  may  be  quite  large  if  u(n)  changes  quickly.  This 
produces  large  errors  in  the  prediction  of  the  engine  response. 
Also,  the  gradient  matrices  must  be  written  as  functions  of  the 
operating  line  variables. 
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C. 3  A  NEW  APPROACH  TO  SIMPLE  ENGINE  MODELS 


Linear  models  are  simple  and  accurate  near  the  linearization 
point.  Extension  of  these  models  to  larger  operating  regions 
should  be  approached  systematically  to  achieve  the  most  accurate 
results.  Steady  state  accuracy  is  extremely  important  in  evaluating 
the  model  results.  A  procedure  for  allowing  accurate  dynamic 
and  static  models  to  be  combined  is  attractive. 


One  approach  to  the  problem  would  be  to  generate  engine  data 
(from  a  simulation)  of  the  form 


(x^  x i ,  u^  i»l,  N 

Then,  the  nonlinear  relationships  shown  in  Eq.  (C.l)  could  be 
fit  with  a  suitable  model  using  curve  fitting  techniques.  Data 
weighting  could  be  used  to  assure  that  points  whete  x»Q  were 
matched  accurately.  The  primary  problem  in  this  procedure  is  that 
weighting  the  static  data  points  degrades  the  dynamic  per'* 
forraance  significantly.  Direct  modeling  of  the  process  results 
in  these  types  of  equations  with  presumably  the  same  problems. 


An  alternate  approach  appears  more  attractive  and  is  based  on 
the  linearization  in  Eq.  (C.6).  The  linearized  matrices, 
fx(x,  u)  are  easily  generated  from  the  nonlinear  simulation. 
These  matrices  can  be  generated  along  the  static  operating  line. 
Sensitivity  calculations  and  multivariate  regression  techniques 
can  be  applied  to  develop  accurate,  and  simple  representations 
of  these  matrices  as  low  order  polynomial  forms.  The  models 
determine  only  the  dynamics  and  there  is  no  trade  off  between 
dynamic  and  static  performance.  These  curve  fitting  techniques 
are  not  pursued  here.  Consider  that  these  functions  have  been 
derived  and  can  be  represented  as  follows; 


ff(x,u) 


.  -  F(x,u) 

x*»0 


(C.9) 
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The  equilibrium  points  are  modeled  by  a  separate  equation. 


*ss  "  S<uss-e> 


(C.10) 


It  is  tempting  to  combine  Eqs.  (C.9)  and  (C - 10)  as  follows: 

X  *  F(x,  uss)(x-xss)  CC. 11) 

and  integrate  away.  A  related  equation  can  be  used  to  accurately 
match  the  true  nonlinear  dynamic  behavior  in  Eq.  (C.  1)- without' 
a  specific  representation  of  the  dynamics,  f(x,u,6). 

C. 4  MODELING  NONLINEAR  DYNAMICS  BY  GRADIENTS 

A  comparison  of  Eqs.  (C.l)  and  (C.ll)  can  be  made  by 
considering  a  Taylor  series  expansion  of  Eq.  (C.l)  around  the 
equilibrium  point  (xgs,  uss)*  First,  define  the  perturbation 


6x  -  x(t)  -  xgs 


Then, 


*  ’  FCxss-  “,*)«*:*  t  H 


fixfix  ♦  . 


(C. 12) 


•  itiEl  5x...4x 


X*X 


ss 


3  x 


u*u 

ss 


x-x 

u*u 


ss 

ss 


where  dyadic  notation  is  used.  This  can  be  more  conveniently 
written  as  follows: 

x  -  (F0  ♦  ^F^fix  ♦  ^F^^S*  4  ^F^xfixdxjdx  / 

♦  llo(s*5)||  cc.13) 

Note  that  Eq.  (C.ll)  represents  the  true  dynamics  tip  to  order, 

| | 6x | 1  .  Suppose  Eq.  (C.ll)  is  modified  so  that  it  matches  ' 
Eq.  (C.ll).  In  this  case,  Eq.  (C.ll)  would  match  the  true 
dynamics  to  order  ||6x||^.  As  it  is,  it  matches  only  to 
order  j|6x||.  Write  Eq.  (C.ll)  in  a  slightly  different  form 
as  follows:  ' 
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x  -  la]F(V^>  +  a2F(Viif»  CC,14) 

An  expansion  of  one  of  the  terms  in  Eq.  (C.14)  would  look 
as  follows: 


6x 


Ft'Vr'  uss>  -  r0TrlqJ  '  Tr257  qj 


F  +F,  r—  +  iF,—  —  + 


(C. 15) 

By  matching  terms  in  Eqs.  (C.14)  and  (C.13)  an  equivalence  can 
be  attained.  There  are  four  parameters,  a^,  ct2,  qj ,  and  q2< 
Thus,  Eq.  (C.14)  can  be  made  to  match  the  first  four  terms  in 
Eq.  (C.13)  -  resulting  in  the  desired  accuracy.  The  coefficients 
yield  the  following  equations: 


“l  *  “2  "  1 

(i_)ai  .  (l-,a2  -  1 


(— i-r)o,  *  (-i-r)a 


2q 


1 


2q. 


(  1  )a  ■,  (l)a 
6<»i  <V 


1 

6 

1 

2? 


The  solution  to  these  equations  is  as  follows: 
1 


?1  M 

a2 


1 

7 


qx  -  3  -  vT  «  1.268 


^2 


•  3  ♦  i/T  •  4.732 


(C. 16) 
(C. 17) 

(C. 18) 
(C. 19) 

(C. 20) 

(C.21) 

(C .22} 

(C. 23) 


the  resulting  form  of  equations  which  match  the  nonlinear  system 
to  order  | } <Sx 1 1  5  can  now  be  written  as  follows: 


.  1 

**  • 


ftPCx 


X*xee 

•» - eJL  q  )  ♦  p(x  ♦—§4  „ 

ss  s./r  ss  ss  wF  ss 


uss)l(x»xss)  (C. 24) 


xss  “  8<uss,0) 


(C. 25) 
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The  simulation  procedure  can  be  derived  by  assuming  u  con* 
stant  on  a  small  interval  and  solving  a  set  of  constant  coef¬ 
ficient  equations  at  each  step.  Higher  order  approximations 
can  be  derived  directly  from  matching  the  expansions. 

Example :  An  explicit  equation 

Suppose  the  engine  model  had  the  following  form  (which 
was  unknown) : 

x  «  x4-2x3u2+u3  (C. 26) 

Gradients  were  available  and  the  linearized  models  were  fit  to 
form  the  ‘’simulated  model”: 

F(x,u)  «(4x3*6x2u2)  (C. 27) 

The  reference  point  schedule  indicates  that  u  *1,  x  *1,  in 

0  9  5 

a  static  equilibrium.  The  "engine  simulation”  for  a  step  input 
to  this  equilibrium  for  any  initial  condition  would  be  as 
follows: 

X  .  l-tFO-’Cl,  1)  *  FO*-^.  D)(x-1)  CC. 28) 

*  3+/3  3-/3 

The  accuracy  of  this  model  can  be  verified  by  expanding 
Eq.  (C.28)  : 

x  -  *  ci*—)3] 

L  3+/j  3-/y 

-  6{(i*-*ll)2  ♦  (i+-5zi)2]  }(X-1) 

3«-/5  3-/5 

■  {x3*x2+x«-l  -  2(x2*x*l))(x-l) 

-  x4-l  -  20c3-l) 

•  x4-2x3+l  (C.29) 
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In  this  case,  the  resulting  equation  (C.29)  is  identical  to  the 
unknown  equation  (C.26)  at  the  equilibrium  because  Eq.  (C.26) 
does  not  have  terms  higher  than  0(6x5)  and  the  form  used  is 
exact.  The  fits  for  the  derivative  matrix,  Eq.  (C.27),  can 
be  implemented  in  subroutine  form  and  evaluation  of  Eq.  (C.28) 
is  straightforward  during  simulation. 

Example  2:  Matrix  Case 

A  complex  nonlinear  set  of  differential  equations  is  listed 
below: 

*1  '  •xi4V*23'xiVvz‘1*Vu2  (c,30) 

*2  "  ‘x2  *xlx2X3X4Ul’XlX4U2  ’U1  u24x3+2ul+^u24 

x3  *  "X3U2^*X1U1U2*X2U2^U1^*X4U1^+^U2^U1^  (C. 32) 

x4  •  *X44U1U2*X2X33U1-X1X2U1U2'X1X4-X4+S  (C. 33) 

These  equations,  for  example,  might  represent  an  approximation 
to  engine  dynamics  derived  from  a  nonlinear  dynamic  simulation. 
These  equations  could  be  replaced  by  the  trim  map  or  steady 
state  schedules  and  a  functional  representation  of  the  time 
derivatives.  The  equivalent  linear  dynamics  function  is  shown 
below: 

r  .  •  •  *1 


.  i 

0 

**j*jwW 

‘VsW1*:* 

wv*iV 

s 

•viwv:* 

..  * 

'V: 

0 

• 

•wy*4 

mm 

VvW t 

(C. 34) 


The  procedure  is  verified  for  the  equilibrium  point, 


X  *X«*X, »x  *U-  aU  *1 
1  2  3  4  1  2 

The  ’’simulation'*  procedure  integrates  the  following  equations: 

X  ^  1  X  ^  1 

i-  i))c*-y  (c-36) 

2  z*jt  3-/r 

for  any  xtoj-x^.  In  this  case,  the  nonlinear  equations  are 
of  lower  order  than  xS  so  the  solutions  are  identical.  The 
equations  can  be  integrated  numerically  to  verify  the  results. 
Notice  that  inaccuracies  in  F(x,u)  do  not  affect  the  equi¬ 
librium  point  of  the  simulation  equations. 
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